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Importance  of  the  cloud  position  and  the  atmospheric  temperature  profile  on  the 
brightness  temperature  values.  In  addition,  investigation  of  the  effects  of  pre¬ 
cipitation  on  the  temperature  retrieval  using  both  the  theoretically  simulated 
values  and  real  SSM/T  data  has  also  been  carried  out  on  the  basis  of  a  statisti¬ 
cal  method.  The  hypothetical  retrieval  experiments  reveal  that  the  temperatures 
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This  finding  is  further  supported  by  the  analysis  employing  the  real  SSM/T  data 
for  a  number  of  case  studies  in  which  temperature  profiles  from  radiosondes  are 
available  for  comparisons. 
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SECTION  1 


INTRODUCTORY  REMARKS 

The  first  microwave  sounder  intended  for  operational  use 
was  flown  aboard  the  Air  Force  Defense  Meteorological  Satellite 
Program  (DMSP)  Block  5D  satellite  system  launched  in  June,  1979. 

The  microwave  sensor  (SSM/T)  is  a  7-channel  scanning  radiometer 
in  the  50-60  GHz  oxygen  band  region.  The  prime  objective  of 
SSM/T  is  to  provide  data  for  deriving  global  basis  temperature 
profiles  in  the  troposphere  and  lower  stratosphere.  Channel  1 
(50.5  GHz)  is  a  window  channel  responding  strongly  to  the  earth's 
surface  characteristics,  dense  clouds,  and  precipitation.  Moreover, 
the  peaks  of  weighting  functions  for  channels  2,  3,  and  4  (53.2, 

54.35,  and  54.9  GHz)  are  below  about  10  km. 

It  Is  known  that  the  prime  advantage  of  microwave  over 
Infrared  temperature  sounders  is  that  the  longer  wavelength  micro- 
waves  are  much  less  affected  by  clouds  and  precipitation.  However, 
to  what  degree  and  quantitative  extent  will  clouds  and  precipita¬ 
tion  Influence  the  brightness  temperature  values  which  would  have  been 
observed  under  clear  conditions?  The  objectives  of  the  present 
project  under  the  sponsorship  of  the  Air  Force  Geophysics  Laboratory 
have  been  (1)  to  develop  a  program  for  the  transfer  of  microwave 
radiation  In  cloudy  and  precipitating  atmospheres,  (2)  to  compare 
the  theoretical  simulated  microwave  and  infrared  radiances  with 
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observed  values  utilizing  the  required  ground  truth  data  when  the 
DMSP  microwave  and  infrared  data  become  available,  and  (3)  to  in¬ 
vestigate  the  inference  of  temperature  profiles  under  all  weather 
conditions. 

Within  the  scope  of  these  objectives  and  subject  to  the 
availability  of  the  DMSP  data,  this  final  report  first  presents  the 
theoretical  foundation  for  the  transfer  of  microwave  radiation  in 
clear  and  cloudy  atmospheres.  Applications  of  the  microwave  pro¬ 
gram  to  DMSP  SSM/T  channels  are  then  described.  Sensitivity  analyses 
on  the  effects  of  clouds  and  precipitation  on  the  upwelling  brightness 
temperatures  for  channels  1  -  4  are  subsequently  carried  out.  In 
addition,  we  report  the  results  of  temperature  profile  retrieval 
exercises  using  the  simulated  brightness  temperatures  computed  from 
the  microwave  radiative  transfer  program  described  previously. 

Lastly,  a  number  of  selected  case  studies  utilizing  the  real  SSM/T 
data  recently  made  available  to  us  are  also  carried  out  for  the 
temperature  retrieval  exercises  by  means  of  a  statistical  method. 

These  case  studies  cover  the  conditions  for  clear,  cloudy  and  pre¬ 
cipitating  atmospheres. 
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SECTION  2 


THEORY  OF  MICROWAVE  RADIATIVE  TRANSFER 
IN  PLANE -PARALLEL  ATMOSPHERES 


2.1  Microwave  Radiative  Transfer  in  Clear  Atmospheres 

Consider  a  non-scattering,  plane-parallel  atmosphere  which 
is  in  local  thermodynamic  equilibrium  and  assume  that  thermal  radiation 
from  the  earth  atmosphere  is  independent  of  the  azimuthal  angle  <|>,  the 
general  equation  of  transfer  may  be  expressed  in  the  frequency  domain 
as 


d  Iv(z,u) 

w  k  pdz 

v 


*  -  Iv(z.n)  +  Bv[T(z)]  , 


(2.1) 


where  Iy  Is  the  radiance  at  frequency  v,  ky  is  the  absorption  coeffi¬ 
cient,  p  is  the  density,  p=cos0,  0  is  the  emergent  angle,  z  is  the 
height,  and  denotes  the  Planck  function.  For  satellite  sounding 
applications,  the  observations  are  normally  taken  close  to  the  up- 
welling  direction  (y~l).  Thus,  for  the  simplicity  of  discussions  on 
the  transfer  of  thermal  radiation  in  a  non-scattering  atmosphere, the 
dependence  on  y  contained  in  the  radiance  expression  is  omitted. 

The  radiance  solution  of  Eq.  (2.1)  at  the  top  of  the  atmosphere 
in  the  upwelling  direction  Is  given  by 

ijL-)  -  iv(o)tv(o,-)  +  jr  bv[t(Z)]  dz»  (2.2) 


where  Iv(0)  represents  the  radiance  contribution  from  the  surface, 


and  the  transmittance  rv(z,«)  is  expressed  with  respect  to  the  top  of 
the  atmosphere  in  the  form 

r  (z,«°)  =  exp[-J  k  (z‘ )p{z'  )dz‘ ]  .  (2.3) 

z 


In  the  microwave  region,  since  the  emissivity  ev  of  the  surface  is 
normally  less  than  unity,  there  will  be  a  reflection  contribution  from 
the  surface.  The  radiance  emitted  from  the  surface  having  a  tempera¬ 
ture  T$  would  therefore  be  given  by 


M°>  =  +  (1'ev}  /  Bv£T(z)] 


3r(0,z) 


9Z 


dz 


(2.4) 


The  first  terms  in  the  right-hand  side  of  Eq.  (2.4)  denotes  the  sur¬ 
face  emission  contribution,  whereas  the  second  term  represents  the 
emission  contribution  from  the  entire  atmosphere  to  the  surface,  which 
is  reflected  back  to  the  atmosphere  at  the  same  frequency.  The  trans¬ 
mittance  rv(0,z)  is  now  expressed  with  respect  to  the  surface. 

Inserting  the  lower  boundary  condition  into  Eq.  (2.2),  we 

find 


IJL»)  =  sJBv(T#)rv(0..)  +  (l-ev)rv(0,-)  /  By[T(z)] 


32V (0*z) 


9Z 


dz 


+ 


3rv( z,») 

3z 


dz  . 


(2.5) 


In  the  frequency  domain,  the  Planck  function  is  given  by 
Bv(T)  *  2h  v3/[c2(ehv/KT-l)]  , 


(2.6) 


where  h  is  the  Planck  constant,  K  is  the  Boltzmann  constant, and  c  is 
the  velocity  of  light.  Moreover,  in  the  microwave  region  hv/KT«l, 


the  Planck  function  may  be  approximated  by 


BV(T)  =  (2Kv2/c2)T  . 


(2.7) 


This  is  the  Rayleigh-Jeans  law,  which  states  that  the  Planck  radiance 
is  linearly  proportional  to  the  temperature.  Furthermore,  radiometers 
which  measure  the  thermal  emission  are  usually  calibrated  with  sources 
at  certain  reference  temperatures.  Thus,  we  may  define  an  equivalent 
brightness  temperature  Tg  such  that 

Iv  -  (2Kv2/c2)Tb(v).  (2.8) 


Substituting  Eqs.  (2.7)  and  (2.8)  into  Eq.  (2.5),  the  solution  of 
microwave  radiative  transfer  may  now  be  written  in  terms  of  tempera¬ 
ture  as  follows: 

0  32V  (0,z) 

Tb(v)  =  e  Tsf  (0,»)  +  (l-ev)ry(0,»)  /  T(z)  - dz 

<30 


+  /  T(z) 
0 


a^z.oo) 

3Z 


dz  . 


(2.9) 


The  transmittance  is  generally  available  with  respect  to 
the  top  of  the  atmosphere;  i.e.,  T^iz)  =  ^(z,®).  Thus,  for  computa¬ 
tional  purposes,  it  is  desirable  to  express  Tv(0,z)  in  terms  of 
2*v( z ,<») .  For  monochromatic  frequencies,  the  transmittance  is  an 
exponential  function  of  the  optical  depth  as  shown  in  Eq.  (2.3). 
Hence,  we  may  write 

Tv(0,z)  =  exp[-^  kv(z,)p(z')dz'] 

■  exp[-jf  kv(z')p(z')dz'  +  J  kv(z' )p(z' )dz' ] 
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(2.10) 


Moreover,  since  T  (0 ,«>) ,  the  transmittance  of  the  entire  atmosphere. 


is  a  constant  value,  we  also  find 

_ _ I  *Tjj Z, 

fm  /  _  ” Tt2”  9Z 


3Tv(0,z) 


[rv(z,»)r 


(2.11) 


Substituting  Eq.  (2.11)  into  Eq.  (2.9),  rearranging  terms,  and  letting 
^(zH  =  rv(z),  Eq.  (2.9)  may  be  rewritten  to  yield 


tb(v)  =  evVv(0)  +  £  Jv(z) 


33'v(z) 

9z 


dz  , 


(2.12) 


where  the  atmospheric  source  term  is  given  by 

Jv(z)  =  (1  +  (l-ev)[!Tv(0)/2’v(z)]2}  T(z)  .  (2.13) 


2.2  Microwave  Radiative  Transfer  in  Scattering  Cloudy  Atmospheres 

The  basic  equation  of  transfer  for  a  plane-parallel  cloud 

layer  consisting  of  absorbing  gases  in  local  thermodynamic  equilibrium 

may  be  written  in  the  form 

d  I  .(t.y)  53  1 

u - - - =  Iv(T,u)  -  -£■  {  Pv(y,vi,)Iv(T,u,)dy' 


-  (1-wv)Bv[T(t)]  .  (2.14) 

In  Eq.  (2.14),  1^  represents  the  monochromatic  radiance  of  frequency 
v,  y  the  cosine  of  the  emergent  angle  with  respect  to  the  zenith,  t 
the  optical  depth  for  cloud  particles  and  gases  within  the  cloud,  P 

v 

the  normalized  axially  symmetrical  scattering  phase  function,  T  the 
cloud  temperature,  and  the  single  scattering  albedo 
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(2.15) 


u)  =  — 

v  8 


s»v 


+  8  +  k  p 

s,v  a,v 


where  8.  „  and  8a  ,  denote  the  volume  scattering  and  absorption  cross 
sections,  respectively,  for  cloud  particles  at  frequency  v.  The  nor¬ 
malized  phase  function  can  be  expressed  as  a  Legendre  polynomial  of 
finite  terms.  By  approximating  the  integration  in  Eq.  (2.14)  utilizing 
Gauss'  quadrature  formula,  a  set  of  first  order  inhomogeneous  differ¬ 
ential  equations  can  be  derived.  Upon  searching  the  homogeneous  and 
particular  solutions  of  the  differential  equations  as  outlined  by 
Chandrasekhar  (1950),  the  complete  solution  of  the  scattered  intensity 
from  an  isothermal  cloud  having  a  temperature  Tc  at  a  given  discrete 
stream  can  be  expressed  as 


MT»V  *  ILm^i)exp(-kmx)  +  MV  »  (2J6) 

m 

where  m  is  the  number  of  discrete  streams  employed,  <|>m  and  km  repre¬ 
sent  eigenfunctions  and  eigenvalues  for  the  differential  equations, 
and  are  associated  with  the  scattering  phase  function  and  single¬ 
scattering  albedo,  and  Lm  are  a  set  of  constants  of  proportionality 
which  can  be  determined  from  the  radiation  boundary  conditions  above 
and  below  the  cloud  layer. 

In  the  previous  subsection,  we  pointed  out  that  in  the 
microwave  region,  the  Planck  function  is  normally  expressed  in  terms 
of  the  temperature  and  the  measured  radiance  is  given  by  an  equivalent 
brightness  temperature.  Thus,  using  Eqs.  (2.7)  and  (2.8),  Eq.  (2.14) 
may  be  written  in  terms  of  the  brightness  temperature  as  follows: 
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(2.17) 


d  Tg(v,x;u) 
u - d7~ 


S)  +1 

Tg(v,x;u)  -  -y  J  Pv(u»u' )Tb(v,t;u' )  du1 


*  (1”5v)T(t)  . 

The  solution  of  Eq.  (2.17)  can  be  derived  in  a  manner  similar  to  that 
of  Eq.  (2.14)  and  is  given  by 


Tg(v,t;ui) 


I  W'iW-V1 


+  T 


m 


(2.18) 


Here  the  constants  of  proportionality  differ  from  those  given  in 
Eq.  (2.16)  and  are  to  be  evaluated  from  the  radiation  boundary  condi¬ 
tions  described  below.  Note  that  Tc«  the  cloud  temperature,  is  in¬ 
dependent  of  frequency  v. 

The  inhomogeneous  and  nonisothermal  structure  of  a  cloud 
layer  may  be  taken  into  consideration  by  dividing  the  cloud  layer 
into  a  number  of  sublayers,  each  of  which  has  a  mean  isothermal 
temperature  and  is  homogeneous  with  respect  to  the  cloud  composition. 
At  the  cloud  top,  the  downward  brightness  temperature  is  equal  to  the 
brightness  contributions  from  every  point  in  the  atmosphere  above  the 
cloud  top.  This  can  be  expressed  by 


TB^v,zt’"pi^ 


oo 

/  T(z)  d2’v(z,zt;-ui)  , 
z=zt 


(2.19) 


where  zt  is  the  height  at  the  cloud  top  and  the  negative  sign  on  ^ 
simply  indicates  downward  transfer.  Note  here  that  we  change  the 
x-coordinate  to  the  height  coordinate  for  the  convenience  of  dis¬ 
cussion.  The  transmittance  in  this  equation  is 


8 


(2.20) 


i  z 

2’v(zt,z;p.)  =  2,v(z,zt;-ui)  =  exp[- —  /  k^z'Jpfz'Jdz'], 

1  zt 

and  In  terms  of  transmittance  with  respect  to  the  top  of  the  atmos¬ 
phere,  we  find 

ry(  z(.,z;yi)  =  2’v(zt,»;y.)/3’v(z,oo;Mi.) 

=  2’v(zt;yi)/3’v(z;yi),  zt>z  .  (2.21) 

Within  the  cloud  layer,  where  scattering  occurs,  continuity  of  bright¬ 
ness  values  from  all  directions  is  required.  Thus, 

TB^v,Z8,;Mi^  *  TB^v*ZX,+liMi^’  <-“l»2,...,N-l  ,  (2.22) 

where  N  is  the  total  number  of  sublayers  within  the  cloud. 

At  the  lower  boundary  of  the  cloud,  three  brightness  con¬ 
tributions  are  immediately  apparent.  These  include  (a)  the  surface 
contribution,  (b)  the  direct  atmospheric  contribution  from  below  the 
cloud,  and  (c)  the  reflected  atmospheric  contribution  from  below  the 
cloud.  A  fourth,  and  perhaps  less  obvious,  brightness  contribution 
at  the  lower  boundary  must  be  considered.  Since  the  reflection  of 
downward  brightness  by  the  earth's  surface  is  significant  for  micro- 
wave  radiation,  the  emergent  brightness  at  the  cloud  bottom  will 
contribute  to  the  lower  boundary  condition.  That  is  to  say,  the 
solution  to  the  radiative  transfer  through  the  cloud  affects  the 
boundary  conditions  used  to  obtain  the  solution.  This  suggests  that 
an  iterative  approach  to  the  correct  solution  is  required.  Practically, 
this  is  accomplished  by  assuming  initially  that  the  top-down  throughput 
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of  the  cloud  equals  one.  The  emergent  intensity  at  the  cloud  bottom 
obtained  by  this  assumption  can  then  be  attenuated  by  two  trips  through 
the  atmosphere  between  the  earth's  surface  and  the  cloud  base  and  by 
the  reflection  at  the  earth's  surface.  This  value  is  added  to  the 
lower  boundary  condition  and  the  solution  to  the  radiative  transfer 
equation  then  produces  a  new  value  for  the  emergent  brightness  temperature 
at  the  cloud  bottom.  From  this  solution  a  new  top/down  throughput  can  be 
calculated  and  the  process  is  repeated  until  the  new  throughput  varies 
from  the  old  by  less  than  one  tenth.  In  all  cases  examined  by  this 
study  the  iteration  halted  after  only  two  steps.  The  lower  boundary 
conditions  can  then  be  expressed  as 

zb 

TB(v,2b;+Mi)  =  SVv(0,zb;+Mi)  +  T(z)drv(z,zb;+ii.) 

zb 

+  ( l-ev)!Tv(0,zb;+ui)C  J  T(z)drv(0,z;-u.) 

+  Tgtv.Z^-y^^^Jr^Zh.Oj-u.)],  (2.23) 

where  T( p^ )  is  defined  to  be  the  top/down  throughput  of  the  cloud  for 
the  stream  defined  by  p..  and  zb  is  the  height  of  the  cloud  base.  It 
Is  given  by 

2\>(p.j)  *  ^B^v*zb’"p1^B^v,zt,_lii^  *  (2.24) 

In  Eq.  (2.23),  all  the  transmittances  can  be  expressed  in  reference 
to  the  top  of  the  atmosphere.  The  numerical  technique  for  solving  a 
set  of  simultaneous  linear  equations  denoted  in  Eqs.  (2.19),  (2.22) 
and  (2.23)  has  been  well  documented  in  our  previous  Infrared  studies 


(Llou  et  al.,  1977;  Feddes  and  Llou,  1977). 

The  microwave  transfer  program  for  cloud  layers  will  give 
the  upward  brightness  temperature  at  the  cloud  top,  i.e.,  Tg(v,2t;ui). 
Thus,  the  upward  brightness  temperature  at  the  satellite  point  of 
view  in  overcast  cloudy  conditions  may  then  be  written  as 

00 

Tgfv,®;^)  *  Tg(v,zt;u1)rv(2t,»;M^)  +  /  T (z)  drv(z,®;iJ.)  , 

1  (2.25) 

where  the  first  term  in  the  right-hand  side  of  Eq.  (2.25)  represents 
the  contribution  from  the  cloud  top  brightness  temperature  which  is 
being  attenuated  to  the  top  of  the  atmosphere,  and  the  remaining 
term  is  the  atmospheric  contribution  above  the  cloud  top. 


SECTION  3 


APPLICATIONS  OF  MICROWAVE  RADIATIVE  TRANSFER  PROGRAM 
TO  DMSP  BLOCK  5D  SSM/T  CH/WNELS 

3.1  Characteristics  of  DMSP  Block  5D  SSM/T  Channels 

The  Passive  Microwave  Temperature  Sounder  (SSM/T)  was 
launched  in  June,  1979,  in  a  sun -synchronous  polar  orbit,  by  the 
United  States  Air  Force  as  part  of  the  Defense  Meteorological 
Satellite  Program  (DMSP)  Block  5D  package.  Other  equipment  of 
meteorological  interest  on  board  the  space  craft  include  visual  and 
infrared  imagery  channels  and  a  scanning  infrared  spectroradiometer. 

The  SSM/T  sensor  is  a  cross  track  scanning  radiometer,  which  acquires 
data  at  32  second  intervals  and  at  seven  angular  positions  separated 
by  12  degrees.  The  "footprints"  of  a  scan  for  the  SSM/T  along  the 
satellite  subtrack  are  shown  in  Figure  1.  The  horizontal  resolution 
is  a  near  circle  of  174  km  diameter  in  the  nadir  direction,  while  it 
is  an  ellipse  with  a  major  axis  of  304  km  and  minor  axis  of  approxi¬ 
mately  213  km  at  the  maximum  scan  anqle  of  36°  from  nadir.  Seven  opera¬ 
tional  frequencies  were  chosen  in  the  vicinity  of  a  strong  oxygen 
absorption  band  In  the  50-60  GHz  region.  Figure  2  depicts  the 
weighting  functions  (3 T/dz)  of  the  SSM/T  channels  in  the  nadir  di¬ 
rection  for  a  surface  emisslvlty  of  0.97.  Channel  1,  50.50  GHz,  is 
a  "window"  channel  and  senses  at  or  near  the  earth's  surface.  Channels 
2,  3,  and  4  (53.20,  54.35  and  54.90  GHz)  have  weighting  functions  that 
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Figure  1.  The  scan  pattern  of  SSM/T. 


peak  at  increasing  heights  in  the  troposphere.  The  weighting  functions 
of  Channels  5,  6  and  7  peak  in  the  stratosphere  with  Channel  5  peaking 
highest  in  the  atmosphere  near  28  km.  The  weighting  functions  result 
from  consideration  of  atmospheric  absorption  due  to  water  vapor  and 
molecular  oxygen.  The  peaks  of  these  weighting  functions  indicate  the 
approximate  location  in  the  atmosphere  from  which  most  of  the  energy 
that  reaches  the  top  of  the  atmosphere  originates.  The  weighting  func¬ 
tions  presented  in  Figure  2  also  include  the  transmittance  corrections 
for  antenna  gain  characteristics.  The  DMSP  Block  5D  satellite  also 
carries  channels  in  the  infrared.  The  Special  Sensor  H  (SSH)  is  the 
infrared  temperature-humidity-ozone  sounder  consisting  of  six  channels 
in  the  15  pm  C02  band,  eight  channels  in  the  18-30  pm  H20  rotational 
band,  and  one  channel  in  the  9.6  pm  0^  band. 

3.2  Atmospheric  Absorption  and  Scattering  Parameters  for 

SSM/T  Frequencies 

Microwave  frequencies  between  50.0  and  60.0  GHz  are  affected 
by  several  different  atmospheric  constituents.  Gases  in  the  atmos¬ 
phere  absorb  and  emit  microwave  radiation.  Cloud  droplets  and 
rain  drops  not  only  absorb  and  emit  but  also  scatter  the  incoming 
radiation. 

3.2.1  Absorption  due  to  gaseous  constituents.  The  primary  gaseous 
absorbers  in  the  microwave  region  are  H20  and  02.  Atmospheric  nitrogen 
gradually  becomes  a  more  significant  absorber  at  frequencies  greater 
than  120  GHz,  but  has  little  influence  in  the  vicinity  of  60  GHz. 

Other  gaseous  constituents  occur  in  very  low  concentrations  and,  in 
general,  their  effects  are  very  small. 
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The  02  absorption  spectrum  in  the  microwave  region  is  due 
to  the  molecule's  magnetic  dipole  moment.  Figure  3  presents  the 
absorption  coefficient  of  dry  air  versus  frequency.  The  absorption 
coefficient  (decibels/km)  shows  a  dramatic  peak  near  60  GHz  and  a 
secondary  maximum  near  120  GHz.  Both  peaks  are  attributable  to  mole¬ 
cular  oxygen.  The  characteristics  of  this  absorption  spectrum  were 
first  investigated  by  Van  Vleck  (1947).  Other  investigators  have 
built  on  his  studies,  culminating  in  a  paper  by  Meeks  and  Lilley 
(1963).  For  this  study,  the  absorption  coefficients  for  02  were 
calculated  using  Meeks  and  Li  1 ley ' s  parameters. 

Water  vapor  absorption  in  the  microwave  is  chiefly  due  to  an 
electric  dipole  moment.  Initial  studies  were  published  by  Van  Vleck 
(1947),  Becker  and  Autler  (1946),  and  King,  et  al.  (1947).  Rosenblum 
(1961)  gave  an  excellent  summary  of  work  done  up  to  1961.  Figure  4 
presents  water  vapor  absorption  (decibels/km)  versus  frequency  in  the 
microwave  region.  It  is  seen  that  the  values  of  absorption  increase  gra¬ 
dually  through  the  range  of  frequencies  under  investigation.  It  is 
apparent  that  the  absorption  due  to  water  vapor  is  less  than  that  of 
molecular  oxygen  for  all  the  channels  chosen  for  the  SSM/T  sounder. 

Note,  however,  that  even  though  the  absorption  coefficients  for  02 
are  generally  higher  than  those  for  HgO,  the  two  values  are  roughly 
equivalent  for  the  window  channel  at  50.50  GHz.  H20  absorption 
coefficients  were  calculated  using  the  method  of  Barrett  and  Chung 
(1962).  Computer  programs  for  determination  of  both  of  these  coeffi¬ 
cients  were  kindly  provided  by  Falcone  (personal  communication). 
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Absorption  coefficient  (decibels/km)  of  dry  air  as  a  function  of 
frequency  (after  Meeks  and  Li 1  ley,  1963).  The  dots  are  measured 
points  by  Hans  Lei  be  during  the  period  from  1974  to  1976  (courtesy 
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Figure  4,  Water  vapor  absorption  (decibel s/km)  as  a  function  of  frequency  in  the 
microwave  region. 


3.2.2  Scattering  and  absorption  due  to  cloud  droplets.  The 
application  of  Mie  scattering  theory  enables  the  calculation  of  optical 
parameters  for  assumed  spherical  particles  given  the  complex  index 
of  refraction,  wavelength  of  incident  radiation,  and  particle  size 
distribution.  Although  larger  rain  drops  are  known  to  deviate  some¬ 
what  from  the  spherical  shape,  Mie  theory  still  remains  probably  the 
best  currently  available  theoretical  tool  for  analyzing  the  effects 
of  cloud  droplets  and  rain  in  the  atmosphere. 

The  Mie  program  used  for  the  calculations  in  this  study  is 
a  modification  for  Uni  vac  computers  of  the  program  by  Liou  and  Hansen 
(1971).  The  output  of  the  Mie  computations  include  the  coefficients 
of  the  Legendre  polynomial  approximation  to  the  phase  function,  the 
extinction  coefficient  (km-1)  and  the  single  scattering  albedo.  Note 
that  the  extinction  coefficient  and  single  scattering  albedo  determined 
here  involve  scattering  and  absorption  by  the  particles  as  the  only 
mechanisms  for  extinction.  These  values  are  modified  to  include 
absorption  by  the  gaseous  atmosphere  within  the  cloud  prior  to  their 
use  in  the  solution  of  the  radiative  transfer  equation. 

The  complex  indices  of  refraction  for  pure  liquid  water  at 
temperatures  of  +10°  and  -10°  C  have  been  determined  by  Savage  (1976) 
for  several  different  frequencies.  These  values  (Table  1)  are  based 
on  a  least  squares  fit  to  the  data  of  Saxton  and  Lane  (1962)  given  by 
Hoi  linger  (1973).  The  complex  indices  of  refraction  for  the  SSM/T 
channels  as  shown  in  Table  2  were  obtained  from  Table  1  by  linear 
interpolation  in  frequency.  For  the  purpose  of  this  study  the  values 
for  0°  C  are  used  throughout. 
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Table  1.  Refractive  index  m=mr-imr  of  pure  liquid  water. 


Frequency 

(GHz) 

T=10°C 

T= 

0°C 

T=-10°C 

mr 

mi 

mr 

mi 

mr 

mi 

19.35 

5.87 

2.96 

5.20 

2.94 

4.64 

2.81 

37.00 

4.32 

2.60 

3.83 

2.37 

3.46 

2.12 

50.30 

3.76 

2.29 

3.37 

2.03 

3.08 

1.78 

89.50 

3.00 

1.67 

2.76 

1.42 

2.60 

1.20 

100.00 

2.89 

1.56 

2.68 

1.31 

2.54 

1.10 

118.00 

2.75 

1.39 

2.58 

1.15 

2.46 

0.96 

130.00 

2.68 

1.30 

2.53 

1.07 

2.43 

0.89 

183.00 

2.49 

1.00 

2.39 

0.81 

2.33 

0.66 

231.00 

2.40 

0.82 

2.33 

0.66 

2.29 

0.53 

Table  2.  Refractive  index  of  pure  liquid  water  (SSM/T). 


Channel  Frequencies  (GHz) 

50.50  53.20  5  4.35  5  4.90  5  8.40  58.82  59.40 


mr  3.37  3.32  3.31  3.30  3.24  3.24  3.23 

mi  2.03  1.98  1.97  1.96  1.90  1.90  1.89 
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Utilizing  the  wavelength  and  the  associated  index  of  refraction, 
Mie  scattering  calculations  for  single  particles  may  be  carried  out.  Once 
the  absorption  and  scattering  efficiencies  for  single  particles  have  been 
determined  in  terms  of  effective  cross  sectional  area,  the  efficiencies  for 
a  size  distribution  of  particles  can  be  determined  by  integrating  over  the 
distribution.  This  area  per  volume  relationship  then  has  dimensions  of 
inverse  length  and  is  described  as  the  absorption  or  scattering  coefficient. 
The  extinction  coefficient  is  then  the  sum  of  the  absorption  and  scattering 
coefficients.  The  drop  size  distributions  used  for  this  study  are  neces¬ 
sarily  numerical  models.  Two  models  have  been  chosen  to  represent  nonpre¬ 
cipitating  and  precipitating  clouds,  respectively.  Deirmendjian's  L-Model 
cloud  (Deirmendjian,  1969)  drop  size  distribution  is  used  for  the  nonpreci¬ 
pitating  case.  It  is  given  by 

N(r)  =  ara  expC-br5)  ,  (3.1) 

where  r  is  the  drop  radius  in  mm  and  note  that  N(r)  is  in  units  of  cm-4. 

By  varying  the  constants  this  model  can  be  made  to  represent  a  variety  of 
cloud  types.  For  this  study,  a=4.975X107,  a=2,  b=l 5 . 1 186 ,  and  6=.5.  The 
selection  of  the  nonprecipitating  cloud  model  here  is  hypothetical  and  we 
fully  realize  that  it  is  not  a  representative  cloud  model  for  typical 
clouds  occurring  in  the  atmosphere  which  have  normally  insignificant  effect 
in  the  microwave  frequencies.  For  the  precipitating  case,  a  theoretical 
drop  size  distribution  based  on  rainfall  rate  is  used  (Marshall  and 
Palmer,  1948).  The  distribution  is  expressed  by 

N(r)  =  0.16  exp(-2r6)  ,  (3.2) 
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where  6=41  R  '  and  R  denotes  rainfall  rate  (mm/hr.).  This  exponential 
behavior  of  drop  size  distribution  has  been  experimentally  verified  by 
Gunn  and  Marshall  (1958)  and  Sekhon  and  Srivastaval  (1970). 

The  phase  function  in  terms  of  Legendre  polynomial  coefficients 
extinction  coefficients,  and  single  scattering  albedos  calculated  for 
each  SSM/T  channel  from  Mie  theory  are  tabulated  in  Tables  3-13.  Since 
the  extinction  coefficient  input  to  the  radiative  transfer  program  must 
be  in  units  of  (km"^),the  final  values  output  by  the  Mie  program  have 
been  multiplied  by  the  proper  constant  to  obtain  these  units. 

3. 3  Atmospheric  Profile  and  Cloud  Models  Used 

Two  climatological  profiles  were  used  for  the  sensitivity 
study.  The  first  profile  was  derived  from  the  northern  hemispheric 
mid-latitude  Spring/Fall  climatology  contained  in  the  U.S.  Standard 
Atmospheric  Supplements,  1966.  The  second  profile  was  derived  from 
the  climatology  for  30°  north  latitude,  July,  contained  in  the  same 
source.  For  both  profiles  temperature  and  height  values  were  inter¬ 
polated  to  40  standard  pressure  levels  used  for  this  study.  Mixing 
ratio  values  for  levels  below  250  mb  were  obtained  from  a  Skew-T  dia¬ 
gram  and  hence  are  saturation  values.  Above  250  mb  a  constant  mixing 
ratio  of  3  parts  per  million  was  used.  The  profiles  for  the  mid¬ 
latitude  Spring/Fall  and  the  30°  N  latitude,  July,  are  depicted  in 
Figure  5. 

To  account  for  the  total  water  within  the  cloud,  the 
vertical  liquid  water  content  is  customarily  used.  For  spherical 
particles  the  droplet  volume  per  unit  volume  of  atmosphere  is  given 
by 
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Table  3.  Optical  parameters  (Delrmendjlan  L*>Model) 
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Table  4.  Optical  parameters  (Marshall -Palmer  Imm/hr) 
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Table  5.  Optical  parameters  (Marshal 1 ^Palmer  2iw«/hr) 
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Table  6.  Optical  parameters  (Marshall -Palmer  3mm/hr) 
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Table  7.  Optical  parameters  (Marshall -Palmer  4mm/hr) 
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Table  8.  Optical  parameters  (Marshall -Palmer  5mm/hr) 


#3  '  *3  #3 


rs.  CM  O 

VO  00  vO 

CO  CM  O'  GO  r«s 

vo  rs.r-s.cD 

to  vo  ov  vo  co 

O  ro  vo  co 

00  ao  o  o  co 

r-  <n  r—  <NJ  CD 


•  1  1  '  •  a  11  it 

UI  Ui  UI  UJ  uj  uj  UJ  J  IjJ  laJ 

vo  co  vo  vo  «r  00  in  >0  wn 

*3"  VO  IAN  r-  rs.  <Ar- 

cd  » —  vo  <">  co  co  v  c  o  io»-ov  i/»<n<n 

tocvjr>*  cvj  vo  < —  or-.  cvj  vo  <o<M<r 

LOVOvO  GO  * —  O  NOn  vmr»  NU)(7« 

CO  XT  CVJ  HO>  V  CVJ  O*  CVJ  r-*-n  OU3  5T 

O- —  cr*  O*— o>  Ocr  O  oino  om  o 

• —  cvj  cd  f—  cvj  co  *—  cv»  « —  cvj  «t  r—  cvj  ^ 


O  K 
v  » —  a# 
13  #3  «a 


o 

1  • 

UJ  uj  Uj 
U3  i£)  rs 
to  <vi  rs 
<vj  tn  rv 
co  vo  cvj 
O'  rs.  cvj 

CVJ  O  O' 
o  CVJ  VO 

Mfv  vr 


O  VO  CVJ 

O  r-  •— 

»-  O  <T 

^«oo> 

«t  O  VO 
r-  von 
lO  CVJ  cr> 


r—  r~  O 

•  I  + 
UJ  UJ  UJ 
COOVfs 

vo  cvj  ao 

mu)rs 
r^cjvn 
O  CVJ  co 
cocorv 
CVJ  CVJ  o 
r-  rv  VJ3 


UJ  UI  UJ 
GO  CVJ  CVJ 
I —  CO  O 

O—  rn 
VO  o  vo 
vrio  io 
vo  O'  w 

•—  CVJ 

r-  h*vn 


CVJr-O 

•  •  4- 

UI  UI  UJ 

00  aa 
O'  O'  O' 
voovrv 
co  co  «r 

CM  CO  VO 

fvnvA 

orvin 


CVJ  r—  O 

I  I  4- 

UI  LtJ  UJ 
<n*-rv 
»-vra) 
rs.  r-  r— 

co  ao  co 
*—  <0  vo 
«r  ov  o 
fvnvo 
cvjrviA 


CVJ  O 

•  t 

UJ  UI  UJ 
O  0%  cn 
00  O'  vo 
Ooo^r 
nvoN 
VO  o  o 
CM*-  fs 

*-ivm 


in 

ro  a»  * — 

*3  *3  r3 


OO  vr  r- 
Cs  *—  vo 
^r 

r>*—  m 
vo*-  n 
kO  CM  VO 
mo** — 
O'  *—  •— 


o--  — 

+  1  • 

UJ  UJ  UJ 
r—  CO  <cr 
n  ov  m 
voois. 

O'  O  CM 

m  cd 
O  -  co 
*—  CM  •— 


O  *x  ^r 

«r  «A(Vi 

O'  co 
VO  GO  VO 
O'  CM  CM 

co  co  m 

O  CM  CO 
*-  <Mr- 


O  — 
4*  1  | 

ui  ui  uj 
O  is.  ^ 
CM  O  VO 

m  o  cm 

• —  CD  CO 

*—  cm  co 

r-CM^. 


O^r- 
«♦»  •  « 
UJ  UI  UI 
VO  O'  CO 
w—  *T 

aj«r  O' 

CM  O  CM 
rs. « —  o 
vrcnrs 
cvi  m  m 
*■—  cm  * — 


*tr  co  co 
ao  vo  co 
ao  •—  cm 

cm  *—  CVJ 
co  rs*  vo 
VO  —  CO 
cm  vo  m 
*~  CM  c— 


O*— *“ 
4*  •  I 

UJ  UI  UI 
O'  « —  CO 
to  O  O' 
r-  lA  sr 


CVJ  vo  to! 
—  cvi^l 


CM  CO  r- 

*3  13  t3 


o*— 

4-  •  • 

UI  UI  UI 

roo- 

OfvCVJ 
O  CO  CM 
»-  (vjrv 
CVJ  LO  . — 
CD  CVJ  CO 
CO  #—  CM 
vf  rvu) 


O  •—  *— 
+■  •  • 
UI  UI  u* 

m  00  cc 

co  co  co 
CO  o  o 
co  cr» 
O'  vo  CD 

mrs.s^ 
CO  —  CVJ 
c  rs*  10 


CD  cr*  *— 

»-o«) 

O'  • —  • — 
CO  O'  O' 

O  vo  VO 

ao  O'  co 

CO  *—  CVJ 
vr  fs  in 


O  — 
4-  1  • 

ui  ui  uj 
vo  cd  m 
vn  m  cm 

cvj  rs.  <-> 

coon 

O'  VO 
CO  O  O' 
O'  CVI  CM 
fs»m 


O  r-  *- 
4*  I  I 
UJ  UI  UI 

cr*  cr* 

OvOn 

KT  CVI  O 

vo  cvj  rs 
ao  ^  cvi 
O'  O'  vO 
O'  CM  CO 
vj-rvvn 


O  •—  r— 
+  •  • 
UI  UI  UJ 
mrsio 
CM  CO  CO 
O  CO  *— 
co  vr  n 
VO  CM  O' 
O  O  VO 
O  CD  CO 
inrvin 


o.—  *— 

4>  I  I 
UJ  UI  UI 

Ovr  vo 
co  vrn 
ao  cd  ao 
conrv 

*ooo<- 

cm  — -  ao 

O  CO  CD 

inrv  in 


#3  #3  c3 


O  — 
4*  «  I 

UJ  UJ  UI 
fs  «T  CD 

cd  m  vo 

CD  O'  O 

•—  O  vo 
O  rs.  cm 
r**  ao  co 

O'  rs.  ud 
CM  *—  *— 


O*-^* 
4  »  1 

UJ  UI  UJ 
O'  C*  VC 
nc  vt 
fs*  fs* 

5T  O'  Lf> 

vo  *r  »— 

CD  O'  rs. 
n  O'  rs 
CD  *—  » — 


vr  n  m 
r*  OM 
cd  ^r  O' 
CM  VO  vo 
O  ao  O' 
VO  vo  CD 
*T  OCD 
m  cm  *— 


O^-  r- 

4  I  I 

UJ  UJ  Uj 

* —  CD  vo 
CD  <M  VO 
CVJ  o  — 
CVJ  O'  UD 
VO*-  OS 
cd  -  rs 
in  *-  qq 
CD  CM  r— 


O 

4-  I  I 
UI  UI  UJ 
OvOCQ 
.O  — 

CVJ  vo  vf 
TT  UD  UD 
to  VO  UD 
r-  CO  CD 
OlDr- 
CM  CM 


O  *“  r— 
-4  •  I 

UI  UJ  UJ 
• —  |S»  lO 
OCViC 
O'  O  in 

vr  Nr- 
*—  CO  r» 

mom 

O  V  r- 
V  CM  CVJ 


CD  CVJ  UD 

O  VO  CO 

CM  CD  CM 
'T  O'  O' 
cvj  vr  o 
cd  m  o 

*~  tr  CVI 
vr  cvj  cvj 


JZ 

O  CM 

-  O  lO  ■ — 

cr  <3  «3  » 3 

a> 


|  UI  UI  UI 

:  O  —  co 

I  O  O'  CD  l 

>  o  vo  rs.  1 

>  o  <n 

>  O  VT  Os 
iOcvj  vr 

I  CD  CVJ  UD 

f—  vo  to 


.♦.•l|  4-  •  •  4*  •  I  I  4-1 

I  UJ  UJ  UJ  I  rj  UI  UJ  UJ  NUJUJUJ  NUIUJI 

;  2  2?  1 3=  o  ud  to  3:  o  vo  co  r;  o  to  < 

2SSy  Or^^OON  <33  0  O  *—  • 
OO'O'OOvOvTCVJOOrlOOM' 
v9r_N®00  O'wr  O  CM  V 
CD  ao  O'  •  CD  UD  **—  .  o  CD  0*1  *  o  CD  I 

on<Mc30vrr>o30<r  *-  ps  o  «cr  < 
,°.ci'au>c><>i,s‘M>OCvjrs  in  O  cvj  j 
f—  VO  VO  r—  VO  vo  r—  vo  loI  r~  VO  ■ 


Table  *9.  Optical  parameters  (Marshall ^Palmer  lOmm/hr) 
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Table  10.  Optical  parameters  (Marshall -Palmer  15mm/hr) 
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Table  11.  Optical  parameters  (Marshall ^Palmer  20mm/hr) 
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Table  12.  Optical  parameters  (.Marshall -Palmer  25mm/hr) 
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Table  13.  Optical  parameters ' (Marshall-Palm6r  30mm/hr) 
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Figure  5.  Climatological  temperature  and  water  vapor  profiles  for  the 
northern  hemispheric  mid-latitude  Spring/Fall  and  30°  north 
latitude,  July. 
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m 


V  =  4  IT  7  r3  N(r)dr  . 
J  0 


(5.1) 


Therefore  the  total  liquid  water  content  per  unit  volume  of  atmosphere 
is 


W  *  4  tp  7  i’''*  N(r)dr  (gm  cm"8)  , 


(5.2) 


where  r  is  in  units  of  cm,  N(r)  is  in  terms  of  cm ”4  and  the  density  of 

_3 

liquid  water  p  is  1  gm  cm  .  For  the  Deirmendjian  L-Model  we  have 


W  =  1.16678(10~2)  (gm  cm"2  km-1)  , 


(5.3) 


while  for  the  Marshall-Palmer  size  distribution  we  find 


W  «  1 . 77883( 10”2)  R'84  (gm  cm"2  km"1)  .  (5.4) 

Table  14  contains  the  total  mass  of  liquid  water  within  a  column  of 
2 

1  cm  cross  section  for  the  cloud  models  used. 


Table  14. 


Total  liquid  water  content  (10 


-2 


O 

gm  cm  ) . 


Thickness  (km) 

1 

2 

3 

4 

5 

Dei  rmendjian 

1.16678 

2.33356 

3.50034 

4.66712 

5.83390 

Marshall-Palmer 

R( mm/hr) 

1 

0.  88941 

1.77882 

2.66823 

3.55764 

4. 44705 

2 

1.59209 

3.18418 

4.77627 

6.36835 

7.96044 

3 

2.23812 

4.47625 

6.71437 

8.95249 

11.1906 

4 

2.84992 

5.69984 

8. 54975 

11.3997 

14.2496 

5 

3.43745 

6.87491 

10.3124 

13.7498 

17.1873 

10 

6.15321 

12.3064 

18.4596 

24.6129 

30.7661 

15 

8.65005 

17.3001 

25.9502 

34.0020 

43.2503 

20 

11.0146 

22.0291 

33.0437 

44.0583 

55.0728 

25 

13.2853 

26.5706 

39.8559 

53.1412 

66.4266 

30 

15.4840 

30.9681 

46.5421 

61.9361 

77.4201 
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SECTION  4 


SENSITIVITY  ANALYSES  ON  THE  EFFECTS  OF  CLOUD  WD 
PRECIPITATION  ON  THE  DMSP  SSM/T  CHANNELS 

In  this  section,  we  present  the  results  computed  from  the 
microwave  radiative  transfer  program  for  cloudy  atmospheres  described 
in  Section  2.  The  computations  utilize  the  single-scattering  para¬ 
meters  presented  in  Section  3  for  input  to  the  transfer  program.  In 
view  of  the  weighting  function  diagram  depicted  in  Figure  2,  it  is 
anticipated  that  only  50.50,  53.20  and  54.35  GHz  (channels  1-3)  whose 
peaks  of  the  weighting  function  are  below  the  tropopause  will  be 
affected  by  clouds  and  precipitation.  However,  in  the  following 
analyses,  54.90  GHz  (channel  4)  is  also  included  in  the  computations. 
Effects  of  the  rainfall  rate  and  cloud  layer  thickness  on  the  up- 
welling  (p-1)  brightness  temperature  over  land  and  ocean  surfaces  are 
first  discussed.  Computational  results  due  to  variations  of  the  cloud 
location  and  atmospheric  profile  are  then  reported.  Physical  inter¬ 
pretation  of  the  results  are  given  in  the  last  Subsection. 

4. 1  Dependence  on  Layer  Thickness  and  Rainfall  Rate 

To  determine  the  effect  of  varying  cloud  thicknesses  on  the 
brightness  temperatures  reaching  the  top  of  the  atmosphere,  cloud  models 
were  inserted  into  the  atmosphere  with  a  constant  cloud  base  of  1  km. 


36 


Thickness  cases  of  1,  2,  3,  4  and  5  km  were  examined.  Variations  of 
the  brightness  temperature  over  land  with  respect  to  cloud  thickness 
and  rainfall  rate  for  channels  1  through  4  are  displayed  in  Figures  6 
through  10.  The  surface  emissivity  over  land  is  taken  to  be  0.97  for 
all  channels  in  this  sensitivity  study.  The  energy  sensed  by  channels 
5,  6  and  7  originates  sufficiently  high  in  the  atmosphere  that  the 
brightness  temperatures  are  insignificantly  affected  by  the  cloud  models 
as  pointed  out  previously. 

The  differences  due  to  thickness  for  channels  3  and  4  may  be 
misleading.  It  is  more  likely  that  a  positional  dependence  for  the 
cloud  is  being  exhibited.  Since  the  cloud  base  was  fixed  at  1  km,  only 
the  thicker  clouds  reached  into  the  energy  peak  source  regions  for 
channels  3  and  4  (weighting  functions  peak  at  about  7  and  10.5  km, 
respectively).  For  example,  a  5  km  thick  cloud  with  a  1  km  base  at  a 

2  mn/hr  rainfall  rate  results  in  brightness  temperatures  for  channels 

3  and  4  of  237.3°  K  and  228.9°  K.  The  Deirmendjian  L-Model  cloud  used 
in  the  position  study  has  liquid  water  content  approximately  equal  to 
that  of  the  2  mm/hr  rain  model  and  yet  a  2  km  thick  L-Model  cloud  with 
base  at  4  km  results  in  nearly  the  same  brightness  temperatures  for 
channels  3  and  4  as  the  5  km  thick  cloud  discussed  above  (237.2°  K  and 
228.7°  K,  respectively). 

Channels  1  and  2  (peaking  near  the  surface  and  2  km)  show 
decreasing  brightness  temperatures  for  increasing  cloud  thicknesses. 

The  decrease  is  greater  for  clouds  modeling  higher  rainfall  rates. 

For  channel  1  at  a  rainfall  rate  of  1  mm/hr  the  decrease  ranges  from 
2°  K.  for  a  1  km  thick  cloud  to  16°  K  for  a  5  km  thick  cloud.  For  a 
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275 


Figure  6.  Brightness  temperature  versus  rainfall  Figure  7.  Brightness  temperature  versus  rainfall 
rate  (or  liquid  water  content)  for  1  km  rate  (or  liquid  water  content)  for  2  km 

thick  cloud  over  land  (mid-latitude  thick  cloud  over  land  (mid-latitude 

Spring/Fall  profile).  Spring/Fall  profile). 
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Figure  8.  Brightness  temperature  versus  rainfall  Figure  9.  Brightness  temperature  versus  rainfall 
rate  (or  liquid  water  content)  for  3  km  rate  (or  liquid  water  content)  for  4  km 

thick  cloud  over  land  (mid-latitude  thick  cloud  over  land  (mid-latitude 

Spring/Fall  profile).  Spring/Fall  profile). 
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water  content)  for  5  km  thick  cloud  over  land  (mid¬ 
latitude  Spring/Fall  profile). 


rainfall  rate  of  30  mm/hr  the  channel  1  brightness  temperature  decreases 
14°  K  for  a  1  km  thick  cloud  and  43°  K  for  a  5  km  thick  cloud.  Channel 

2  displays  the  same  trends  as  channel  1  but  to  a  lesser  degree.  For  a 
rainfall  rate  of  1  mm/hr  the  decrease  in  channel  2  brightness  tempera¬ 
ture  as  a  function  of  cloud  thickness  ranges  from  less  than  0.5°  K  to 
5°  K.  For  a  rainfall  rate  of  30  mm/hr  the  range  is  4°  K  to  22°  K. 

The  results  for  channels  3  and  4  show  relatively  constant 
brightness  temperatures  unaffected  by  rainfall  rates  less  than  about 

3  irm/hr.  For  rainfall  rates  greater  than  3  mm/hr  the  brightness 
temperature  decreases  with  increasing  cloud  thickness.  The  maximum 
decreases  are  7°  K  for  channel  3  and  2°  K  for  channel  4;  both  at  rain¬ 
fall  rates  of  30  mm/hr  and  cloud  thicknesses  of  5  km. 

Figures  11  through  14  depict  the  variation  of  brightness 
temperatures  over  ocean  having  an  emissivity  of  0.51  with  respect  to 
cloud  thickness  and  liquid  water  content  for  channels  1-4.  The  results 
over  the  ocean  are  quite  different  from  the  results  over  land.  Note 
that  the  only  physical  difference  between  these  two  cases  is  that  the 
surface  emissivity  for  all  channels  is  taken  to  be  0.97  for  a  land 
surface  and  0.51  for  an  ocean  surface.  The  difference  in  the  surface 
emissivity  reveals  that  the  surface  contribution  over  ocean  is  approxi¬ 
mately  one  half  of  the  surface  contribution  over  land,  and  that  the 
atmospheric  contribution  reflected  from  the  earth's  surface  is  approxi¬ 
mately  16  times  greater  over  ocean. 

Figures  11  and  12  show  that  channels  1  and  2  follow  similar 
trends  over  the  ocean  but  to  different  degrees.  For  all  cloud  thick¬ 
nesses,  an  Initial  reduction  In  the  brightness  temperature  is  observed 
for  clouds  with  low  liquid  water  content.  For  channel  1  the  maximum 
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Figure  12.  53.20  GHz  brightness  temperature  as  a  function  of  cloud 
thickness  and  liquid  water  content  (mid-latitude  Spring/ 
Fall  profile  over  ocean).  Values  of  the  cloud  liquid  wa 
content  in  this  figure  should  be  multiplied  by  10"^. 
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Figure  13.  54.35  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (mid-latitude  Spring/ 
Fall  profile  over  ocean). 
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Figure  14.  54.90  GHz  brightness  temperature  as  a  function  of  cloud 

thickness  and  liquid  water  content  (mid-latitude  Spring/ 
Fall  profile  over  ocean). 
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reduction  is  2°  K,  while  for  channel  2  the  maximum  reduction  is  nearly 

-2  -2 

13°  K.  At  some  liquid  water  content  less  than  5X10  gm  cm  the 
brightness  temperature  begins  to  increase.  The  turning  point  occurs 
at  lower  liquid  water  content  for  thinner  clouds.  The  brightness  tem¬ 
perature  then  increases  until  the  liquid  water  content  reaches  about 
-2  -2 

15X10  gm  cm  .  For  channel  1  this  increase  results  in  a  maximum 

brightness  temperature  more  than  35°  K  higher  than  the  clear  column 

value.  For  channel  2  the  maximum  brightness  temperature  is  about  3°  K 

greater  than  the  clear  column  value.  For  liquid  water  contents  greater 
-2  -2 

than  15X10  gm  cm  ,  the  brightness  temperature  decreases  for  both  channels 
and  all  cloud  thicknesses.  Note  that  the  maximum  brightness  tempera¬ 
ture  observed  over  ocean  is  very  nearly  the  same  as  the  minimum  bright¬ 
ness  temperature  observed  over  land  (within  about  4  -  7°K). 

Figures  13  and  14  indicate  that  channels  3  and  4  are  relatively 
unaffected  by  a  1  km  thick  cloud  with  a  cloud  base  at  1  km.  Increasing 
the  cloud  thickness  above  1  km  causes  steadily  decreasing  brightness 
temperatures.  For  a  constant  cloud  thickness,  increasing  the  liquid 
water  content  causes  a  decrease  in  the  brightness  temperature.  Note 
that  channel  4  is  affected  by  less  than  1°  K  for  all  of  the  cases 
studied. 

4.2  Dependence  on  Layer  Location  and  Atmospheric  Profile 

In  this  subsection,  we  first  investigate  the  effect  of  the 
position  of  the  cloud  layer  in  the  atmosphere  on  the  upwelling  bright¬ 
ness  temperature.  Figures  15  through  20  display  the  importance  of  cloud 
position  relative  to  the  peak  of  the  weighting  function  for  each  channel 
over  the  land  surface. 
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dependence  on  cloud  base  height  for  dependence  on  cloud  base  height  for 

a  2  km  thick  Deirmendjian  L-Model  a  2  km  thick  Deirmendjian  L-Model 

cloud  (mid-latitude  Spring/Fall  cloud  (mid-latitude  Spring/Fall 

profile  over  land).  profile  over  land). 


Figure  77.  54.35  GHz  brightness  temperature  Figure  18.  54.90  GHz  brightness  temperature 

dependence  on  cloud  base  height  for  dependence  on  cloud  base  height  for 

a  2  km  thick  Deirmendjian  L-Model  a  2  km  thick  Deirmendjian  L-Model 

cloud  (mid-latitude  Spring/Fall  cloud  (mid-latitude  Spring/Fall 

profile  over  land).  profile  over  land). 


Channels  1  and  2  show  the  effects  on  a  channel  when  the  cloud 
is  at  or  above  the  peak  of  the  weighting  function,  i.e.,  the  prime  energy 
source  region  for  the  channel.  The  brightness  temperature  decreases 
rapidly  as  the  cloud  is  moved  higher  putting  more  and  more  of  the  energy 
source  for  the  channel  beneath  it.  Eventually,  the  brightness  temperature 
becomes  near  constant  as  all  of  the  significant  source  region  for  channel 
energy  is  below  the  cloud  and  therefore  raising  the  cloud  higher  has  only 
slight  effects. 

Channels  3  and  4  show  the  results  of  moving  a  cloud  from  below 
the  prime  energy  source  region  up  through  the  source  region.  Far  below 
the  energy  source  the  brightness  temperature  remains  nearly  constant. 

Then,  as  the  cloud  moves  into  the  source  region,  the  brightness  tempera¬ 
ture  decreases  significantly,  becoming  near  constant  above  the  source 
region  as  for  channels  1  and  2. 

Channels  6  and  7  indicate  that  for  channels  peaking  high 
enough  in  the  atmosphere  to  be  free  of  surface  effects,  brightness 
temperatures  increase  slightly  as  the  cloud  approaches  the  energy 
source  region  from  below. 

Sensitivity  analyses  on  the  positional  dependence  were  also 
carried  out  over  the  ocean  surface.  It  was  found  that  the  same  positional 
dependence  exists  over  the  ocean  as  over  the  land  with  no  change  in  trends 
for  any  of  the  channels. 

Figures  21  through  28  display  the  results  when  the  sensitivity 
analysis  is  duplicated  for  the  30°  N  latitude,  July,  profile.  Although 
the  different  profile  naturally  results  in  different  brightness  tempera¬ 
tures,  the  trends  noted  for  the  northern  hemispheric  mid-latitude 

Spring/Fall  profile  for  liquid  water  content  greater  than  about  3X1 0~2  gm 
-2 

cm  are  observed  to  persist. 
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over  land. 
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L3  Interpretation  of  Results 

In  order  to  aid  in  the  interpretation  of  results  Figure  29  is 
included  which  graphically  depicts  the  terms  contributing  to  the  cloud 
boundary  conditions  and  the  mechanisms  for  energy  loss  and  gain  within 
the  cloud.  There  are  three  mechanisms  for  energy  gain  within  the  cloud 
caused  by  gaseous  emission,  droplet  emission, and  multiple  scattering. 

On  the  other  hand,  there  are  three  mechanisms  for  energy  loss  due  to 
gaseous  absorption,  droplet  absorption, and  single  scattering.  The 
balance  of  these  mechanisms  determines  whether  the  emergent  energy  at 
the  cloud  top  and  bottom  is  greater  or  less  than  the  respective  boundary 
condition  at  the  cloud  bottom  and  top.  Of  course,  the  boundary  condi¬ 
tions  themselves  are  an  important  contribution  to  the  emergent  energy 
at  the  opposite  side  of  the  cloud.  The  upper  boundary  condition  is 
totally  unaffected  by  the  earth's  surface.  For  the  lower  boundary 
condition,  quite  the  opposite  is  true.  When  the  emissivity  (and  hence 
the  reflectivity)  of  the  earth's  surface  varies,  three  of  the  four 
contributing  terms  for  the  lower  boundary  condition  are  influenced. 

These  include  (1)  the  surface  emission  term,  (2)  the  atmospheric  con¬ 
tribution  reflected  from  the  surface  of  the  earth,  and  (3)  the  emer¬ 
gent  energy  from  the  cloud  bottom  reflected  from  the  earth's  surface. 
Finally,  it  is  important  to  note  that  the  emission  terms  are  generally 
stronqer  near  the  cloud  bottom  due  to  the  temperature  gradient  within 
the  cloud;  and  throughout  the  cloud,  emission  by  water  droplets  is 
greater  than  emission  by  atmospheric  molecules. 

Introducing  a  cloud  into  the  atmosphere  will  then  result  in 
increased  emission  within  the  cloud  layer  and  decreased  transmission 
through  the  cloud  layer.  The  three  contributions  to  the  cloud  top 
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Figure  29.  Radiative  transfer  through  a  cloud  layer. 
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brightness  temperature  are  emission  by  the  earth's  surface,  atmospheric 
emission  and  emission  from  within  the  cloud  layer. 

To  remove  the  surface  term,  atmospheric  contribution,  or 
cloud  emission,  we  simply  set  the  surface  temperature,  boundary 
conditions, or  cloud  temperature  equal  to  zero,  respectively.  Each 
contribution  can  be  analyzed  separately  by  removing  the  other  two. 

In  Figures  30  through  33  the  effects  of  cloud  thickness 
and  liquid  water  content  on  the  three  contributions  to  the  cloud  top 
brightness  temperature  over  the  ocean  are  displayed.  The  solid  line 
is  the  graphical  sum  of  these  components  and  the  X's  are  computed  cloud 
top  temperatures  directly  from  the  microwave  transfer  program  for 
randomly  selected  cases.  Clearly,  for  channel  1  the  addition  of  energy 
by  droplet  emission  is  greater  than  the  energy  loss  due  to  droplet 
absorption  and  single  scattering  for  liquid  water  contents  greater  than 
0.01  gm  cm  .  This  explains  the  region  of  increasing  brightness  tempera¬ 
ture  for  channel  1  observed  in  Figure  11.  For  liquid  water  content 
greater  than  0.20  gm  cm  the  droplet  emission  reaching  the  cloud  top 
is  nearly  constant  and  represents  almost  all  of  the  energy  reaching 
the  cloud  top.  Although  the  surface  emission  and  atmospheric  contri¬ 
butions  to  cloud  top  brightness  temperature  are  small  in  this  region, 
they  are  decreasing  with  increasing  liquid  water  content  so  the  cloud  top 
brightness  temperature  decreases  slightly  also.  This  explains  the  slight  de 
crease  in  the  brightness  temperature  for  channel  1  in  this  region.  The 

fact  that  cloud  emission  accounts  for  almost  all  of  the  cloud  top  bright- 

_o 

ness  temperature  for  water  content  greater  than  0.20  gm  cm  explains 
why  the  brightness  temperatures  over  ocean  and  over  land  are  nearly  the 
same  in  this  region.  The  cloud  emission  must  be  less  than  the  reduction 
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cloud  top  brightness  for  a  1  km  thick  cloud  top  brightness  for  a  3  km  thick 
cloud  over  ocean  (mid-latitude  Spring/  cloud  over  ocean  (mid-latitude  Spring/ 
Fall  profile).  Fall  profile). 
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in  the  surface  and  atmospheric  terms  for  liquid  water  content  less  than 

_2 

0.01  gm  cm  in  order  to  explain  the  slight  decrease  in  brightness  tempera¬ 
tures  in  this  region. 

The  cloud  emission  curves  remain  the  same  over  land.  As  previously 
noted,  the  surface  term  is  nearly  twice  as  large  over  land  while  the  reflected 
terms  are  only  one  sixteenth  as  large.  The  direct  atmospheric  contribution 
remains  the  same.  The  net  result  is  that  the  upwelling  brightness  tempera¬ 
ture  at  the  cloud  bottom  over  land  is  significantly  larger  than  that  over 
ocean.  For  example,  our  numerical  experiments  show  that  the  lower  boundary 
condition  over  land  for  a  1  km  thick  cloud  ranges  from  280°  K  to  285°  K  for 
rainfall  rates  of  1  to  30  mm/hr.  Over  ocean  the  range  is  200°  K  to  240°  K 
for  the  same  cloud  thickness  and  rainfall  rates.  Since  the  loss  of  energy 
within  the  cloud  layer  is  proportional  to  the  energy  incident  at  the  cloud 
base,  a  greater  reduction  takes  place  over  land.  This  increased  reduction 
is  associated  with  the  continuously  decrease  of  brightness  temperatures  with 
increasing  liquid  water  content  over  land  for  channel  1  as  shown  in  Figures 
6  through  10.  In  these  cases,  effects  of  the  cloud  transmittance  is  larger 
than  the  cloud  emission. 

For  channel  2  the  surface  contribution  to  the  cloud  top  brightness 
temperature  over  ocean  is  less  than  that  for  channel  1  due  to  the  differences 
in  transmittance  as  illustrated  in  Figures  30  and  32  and  Figures  31  and  33 
for  1  km  and  3  km  cases,  respectively.  The  atmospheric  terms  are  increased, 
however,  due  to  the  shape  of  the  weighting  function  for  this  channel  and  the 
high  reflectivity  of  the  ocean  surface.  Although  the  cloud  emission  contri¬ 
bution  for  the  3  km  case  is  slightly  greater  for  channel  1  over  channel  2, 
the  net  result  is  a  higher  cloud  top  brightness  temperature  for  channel  2. 

The  lower  boundary  condition  is  also  higher  for  the  same  reasons,  resulting 
in  a  greater  reduction  in  the  surface  and  atmospheric  terms  as  they  traverse 
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the  cloud.  Thus,  for  liquid  water  contents  less  than  about  0.05  gm  cm  (see 
Figure  12)  at  which  the  cloud  emission  is  small,  we  observe  the  deepening  of 
the  reduction  in  the  brightness  temperature.  The  greater  reduction  in  the 
surface  and  atmospheric  terms  also  explains  why  the  increase  in  brightness 
temperature  for  liquid  water  content  of  0.05  through  0.15  gm  cm  is  not  as 
large  as  the  corresponding  increases  for  channel  1.  Note  that  the  cloud 
emission  is  only  slightly  larger  for  channel  1  than  channel  2  in  the  3  km  case 
and  for  the  1  km  case  it  is  approximately  the  same  for  both  channels. 

Based  on  the  numerical  experiments  carried  out  for  the  study  of  the 
effects  of  the  liquid  water  content  on  the  surface,  atmospheric,  and  cloud 
emission  contributions  to  the  cloud  top  brightness  temperature  over  land,  we 
find  that  nearly  the  same  discussion  applies  for  channel  2  over  land  as  for 
channel  1.  Thus,  no  plots  similar  to  Figures  30  -  33  will  be  presented  but 
rather  some  general  discussions  are  made  here.  The  surface  terms  increases  due 
to  the  change  in  emissivity  of  the  earth's  surface.  The  increase  is  smaller  than 
for  channel  1,  however,  due  to  transmittance  differences.  The  atmospheric  terms 
decrease  due  to  the  change  in  reflectivity  of  the  earth's  surface.  The  decrease 
is  larger  than  for  channel  1  due  to  the  shape  of  the  weighting  functions.  The 
net  result  is  that  the  lower  boundary  condition  increases  for  channel  2  over 
land,  but  the  increase  is  much  reduced  for  channel  1.  For  example,  for  a  1  km 
thick  cloud  and  rainfall  rates  of  1  and  30  mm/hr  the  lower  boundary  condition 
for  channel  2  is  nearly  constant  at  285°  K  over  land  but  ranges  from  250°  K  to 
265°  K  over  ocean.  Once  again,  this  increase  in  the  lower  boundary  condition 
results  in  increased  reduction  through  the  cloud.  Thus,  we  find  that  the  bright¬ 
ness  temperatures  are  continuously  decreasing  with  increasing  rainfall  rates. 

Channels  3  and  4  are  relatively  unaffected  by  the  emissivity  of  the 
earth's  surface.  This  is  due  to  the  fact  that  the  transmittance  from  the  earth's 
surface  to  the  top  of  the  atmosphere  is  only  0.023  for  channel  3  and  0.003  for 
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channel  4  based  on  the  results  from  the  transmittance  program.  This  fact  ex¬ 
plains  why  the  brightness  temperatures  for  channels  3  and  4  do  not  vary  much 
between  land  and  ocean  surfaces.  In  this  study  the  highest  cloud  top  modeled 
was  set  at  6  km.  This  is  just  below  the  peak  of  the  weighting  function  for 
channel  3  (7  km),  but  still  well  below  the  peak  for  channel  4  (10.5  km).  As  the 
cloud  top  moves  into  the  energy  source  regions  for  these  channels,  emission  by 
the  cloud  is  lower  since  the  temperatures  are  lower  at  higher  altitudes.  Fur¬ 
thermore,  the  transmittance  from  6  km  to  the  top  of  the  atmosphere  is  less  than 
0.3  for  both  channels.  Therefore,  within  the  energy  source  region  occupied  by 
the  upper  portion  of  the  cloud,  the  extinction  mechanisms  dominate  over  the 
mechanisms  for  the  addition  of  energy,  leading  to  slight  decreases  in  the  bright¬ 
ness  temperature. 

The  interpretation  of  the  cloud  position  dependence  illustrated  by 
Figures  15  through  20  is  quite  straightforward.  Note  that  the  lowest  cloud  base 
examined  is  at  1  km  where  the  brightness  temperatures  over  land  for  channels  1 
and  2  are  265.5°  K  and  255.2°  K,  respectively.  These  values  are  higher  than 
clear  column  values  due  to  the  fact  that  cloud  emission  exceeds  extinction  for 

a  2  km  thick  cloud  with  a  base  at  1  km  for  liquid  water  contents  greater  than 

-2 

0.01  gm  cm  (see  Figures  30  through  33).  The  liquid  water  content  for  the  2  km 
thick  L-Model  cloud  is  0.23  gm  cm  .  Thus,  as  the  cloud  is  moved  higher  in  the 
atmosphere,  two  changes  take  place.  The  emission  by  the  cloud  reduces  due  to 
the  decrease  in  the  cloud  temperature.  Meanwhile,  the  lower  boundary  condition 
increases  because  a  greater  portion  of  the  energy  source  region  for  channels  1 
and  2  is  below  the  cloud  base.  This  results  in  the  increased  energy  reduction 
within  the  cloud  layer.  The  net  result  is  that  brightness  temperatures  for 
channels  1  and  2  decrease  as  the  cloud  is  moved  higher  up  in  the  atmosphere. 
Eventually,  the  cloud  is  raised  high  enough  that  nearly  all  of  the  energy  source 
regions  for  channels  1  and  2  are  already  below  the  cloud  and  cloud 
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temperatures  are  sufficiently  low  that  emission  by  the  cloud  is  very 
small.  Then  moving  the  cloud  still  higher  reveals  very  little  effect 
so  that  brightness  temperatures  are  nearly  constant.  The  same  discus¬ 
sion  applies  over ocean  except  that  over ocean  even  the  1  km  cloud  base 
shows  a  reduction  in  brightness  temperatures  as  compared  with  the  clear 
column  values. 

Channels  3,  4,  6  and  7  have  weighting  functions  which  peak 
higher  in  the  atmosphere,  reducing  the  effects  of  surface  emissivity. 
Therefore,  the  effects  of  a  2  km  thick  cloud  are  basically  the  same 
over  ocean  and  land  surfaces.  For  clouds  well  below  the  peak  of  the 
weighting  functions  they  show  no  effect  on  the  brightness  temperature. 
This  is  because  there  is  no  significant  energy  source  below  the  cloud 
to  be  reduced  by  the  cloud  extinction.  Also,  the  transmittance  from 
the  cloud  top  to  the  top  of  the  atmosphere  is  sufficiently  small  that 
emission  by  the  cloud  may  be  negligible.  As  the  cloud  moves  into  the 
energy  source  region,  an  intricate  trade-off  of  cloud  emission  and  ex¬ 
tinction  takes  place.  For  channels  3  and  4  this  trade  off  holds  the 
brightness  temperatures  very  nearly  constant  until  the  cloud  is  well 
into  the  energy  source  region.  For  channels  6  and  7  a  very  slight 
increase  in  the  brightness  temperatures  is  noted.  This  increase  is 
less  than  0.2°  K,  however,  which  is  below  the  noise  level  of  the  SSM/T. 
Finally,  as  the  cloud  is  moved  above  the  peaks  of  the  weighting  func¬ 
tions,  brightness  temperatures  decrease,  eventually  becoming  near 
constant  as  for  channels  1  and  2. 
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SECTION  5 


TEMPERATURE  PROFILE  RETRIEVAL  EXERCISES 

The  main  objective  of  the  DMSP  SSM/T  microwave  sounders  has 
been  to  derive  accurate  temperature orof lies  in  all  weather  conditions 
for  operational  use.  The  prime  advantage  of  microwave  temperature 
sounders  over  infrared  sounders  is  that  the  longer  microwaves  are 
much  less  affected  by  clouds  and  precipitation.  It  seems  therefore 
important  to  investigate  the  effects  of  clouds  and  precipitation  on 
the  temperature  retrieval  program.  In  this  section,  we  first  present 
the  hypothetical  temperature  inversion  in  precipitating  atmospheres 
using  the  brightness  temperatures  calculated  from  the  microwave  ra¬ 
diative  transfer  program.  The  retrieval  program  adopted  in  this  study 
is  the  statistical  method  developed  for  the  Air  Force  Global  Weather 
Central  intended  for  operational  use.  We  also  report  a  number  of 
case  studies  using  the  real  OMSP  SSM/T  data  that  we  recently  secured 
for  clear,  cloudy  and  precipitating  cases. 

5.1 

The  program  that  we  use  for  temperature  retrieval  exercises 
is  based  on  the  statistical  method  described  by  Rigone  and  Stogryn 
(1977).  The  computer  package  for  the  statistical  method  was  kindly 
provided  to  us  by  the  Air  Force  Global  Weather  Central. 

In  the  statistical  method,  the  surface  emissivity  effect  is 


Temperature  Profile  Retrieval  Using  Simulated  Brightness 
temperatures 
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first  removed  so  that  the  retrieval  method  could  be  applied  to  all 
surface  conditions.  For  the  purpose  of  outlining  the  method,  we 
def i ne 

*  M*  (z,») 

Tu(v)  =  /  T(z)  - dz  ,  (5.1) 

and 

0  M*( O.z) 

Td(v)  =  /  T(z)  — ^ - dz  ,  (5.2) 

oo 

so  that  Eq.  (2.9)  can  be  rewritten  as  follows: 

tB(v)  =  %Vv(0)  n-Td(v)/Ts]  +  Ta(v)  ,  (5.3) 

where 

Ta(v)  =  Tu(v)  +  Td(v)  Tv(0)  .  (5.4) 

In  Eq.  (5.3),  the  second  term  in  the  right-hand  side  denotes  the  con¬ 
tribution  to  the  upwelling  brightness  temperature  caused  by  the  atmos¬ 
phere  only,  and  the  surface  effects  are  contained  in  the  first  term. 
Since  channel  1  centered  at  50.5  GHz  has  a  weighting  function  peak  at 
the  surface,  it  is  utilizied  in  the  context  of  removing  the  surface 
contribution  for  other  channels.  Based  on  Eq.  (5.3),  we  may  define 
the  contribution  to  the  brightness  temperature  caused  by  the  atmosphere 
only  for  channels  2-7  in  the  form 

W  =  W  '  EW  '  Ta(vl)]  a(vj>’  j=2 ,3, ' ' '  ,7  ,  (5.5) 


where 
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a^) 


%1Ts%(Ps)C,-V''i)/TsJ 

’  a"d  a(v>’  *  ’  ' 


In  the  statistical  approach,  it  is  generally  assumed  that 
the  derivation  of  the  predicted  parameter  Ti  from  the  climatological 
mean  may  be  expressed  as  a  linear  combination  of  the  deviation  of 
the  measured  data.  Upon  finding  a  linear  operator  D  which  will  yield 
a  minimum  mean  square  deviation  of  the  predicted  tenperature  profile 

A 

from  the  true  temperature  profile  in  a  statistical  sense,  the 
predicted  temperature  profile  may  be  obtained.  The  linear  operator, 
called  the  predictor  matrix,  may  be  expressed  in  terms  of  a  covariance 
matrix,  which  can  be  constructed  experimentally  by  collecting  coinci¬ 
dences  of  radiances  derived  from  remote  sounders  with  temperature 
values  obtained  from  direct  soundings.  The  measured  data,  in  the 

A 

present  case  is  T  given  by  Eq.  (5.5).  Thus,  we  write 

d 

<VV  ■  IVVV 

J 

J 


0ij\j^Bl  ?°ijaj  +  lDij(^alaj'Taj)  * 


(5.6) 


where  Ta  denotes  the  mean  measured  value  and  we  note  that  TB1  is  not 
defined  in  Eq.  (5.5),  and  so  the  first  terms  contain  j=2,***,7.  In 
matrix  notations,  we  find 

T  =  D'  Tb  +  R  ,  (5.7) 
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where 


R  =  T  +  Tai  D  a  -  D  ^  , 

and  the  linear  operator  D'  is  a  matrix  whose  first  column  is  -D  a  and 
whose  remaining  columns  are  the  columns  of  D.  It  is  clear  that  the 
retrieval  technique  contains  elements  depending  mainly  on  the  atmos¬ 
phere  but  not  on  the  surface,  and  so  it  should  be  valid  over  land, 
water,  or  mixed  surface  conditions.  As  pointed  out  previously,  the 
D  and  R  may  be  determined  from  a  large  number  of  upper  air  soundings 
for  a  wide  range  of  meteorological  conditions  which  have  been  achieved 
over  the  years  and  the  brightness  temperatures  calculated  for  a  given 
atmosphere. 

Shown  in  Figure  34  is  an  exercise  of  temperature  retrieval 
using  the  statistical  covariance  method.  The  mid-latitude  Spring/Fall 
profile  of  a  standard  atmosphere  (solid  curve)  is  used  and  the  observed 
brightness  temperatures  used  for  the  seven  SSM/T  channels  are  values 
theoretically  calculated.  The  exercise  has  been  carried  out  for  cases 
over  ocean  and  land.  It  is  apparent  that  the  procedure  outlined  above 
has  very  successfully  removed  surface  effects  from  the  temperature  re¬ 
trieval.  Also  shown  are  the  temperature  retrievals  when  a  2  km  thick 
precipitation  layer  with  a  base  height  set  at  1  km,  having  various 
rainfall  rates,  have  been  added  to  the  atmosphere.  It  is  seen  that 
the  surface  temperature  suffers  increased  degradation  as  the  rainfall 
rate  increases.  Based  on  these  analyses,  it  seems  that  large  errors 
in  the  recovered  temperature  profile  may  be  anticipated,  even  with 
microwave  sounders,  when  the  atmosphere  within  the  satellite  field-of- 
view  contains  precipitation  and  heavy  clouds. 
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TEMPERATURE  (°K) 

Figure  34.  Hypothetical  temperature  retrieval  exercise  over  land  and  ocean  using  the  statistical 
covariance  method. 


5.2  Temperature  Profile  Retrieval  Using  DMSP  SSM/T  Data 

The  statistical  method  for  temperature  retrieval  described 
in  Section  5.1  was  applied  to  a  number  of  cases  where  SSM/T  and  radio¬ 
sonde  data  were  both  available.  Two  days,  30  October  1979  and  23 
November  1979,  were  chosen  during  which  significant  cloud  and  precipi¬ 
tation  events  were  present  over  the  continental  United  States.  The 
cases  selected  are  presented  in  Table  15  where  the  satellite  pass 
times  are  the  actual  observation  times  of  the  SSM/T  instrument.  We 
find  that  the  scan  times  are  generally  between  0000  Zulu  and  0600  Zulu 
on  these  two  days.  Thus,  the  0000  Zulu  radiosonde  observations  were 
deemed  most  representative  and  were  used  in  the  comparisons. 

Figure  35  shows  the  retrieved  temperature  profiles  for  the 
four  clear  cases;  two  on  30  October  1979  and  two  on  23  November  1979. 
Except  for  the  Greensboro  case,  the  retrieved  temperature  profiles 
when  they  are  compared  with  nearby  radiosonde  data  appear  to  be  reason 
ably  good  in  view  of  the  statistical  method  used.  The  failure  in  the 
retrieval  for  the  Greensboro  case  seems  largely  due  to  the  fluctuated 
temperature  profile  that  occurred  in  the  atmosphere.  Generally,  we 
found  that  the  statistical  method  is  working  properly  when  the  actual 
profile  is  smooth  and  when  no  inversion  is  present. 

The  retrieved  temperature  profiles  under  cloudy  conditions 
are  illustrated  in  Figure  36.  Apparently,  the  temperature  profile 
retrieval  program  using  the  microwave  frequencies  in  the  60  GHz  oxygen 
band  is  affected  insignificantly  by  non-precipitating  clouds.  Com¬ 
pared  with  the  temperatures  obtained  from  radiosondes,  the  retrieved 
patterns  involving  clouds  are  similar  to  those  under  clear  conditions. 
Note  that  in  each  diagram,  the  percentage  of  cloud  cover  is  depicted. 
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Table  15.  Selected  cases 


Station 

Name 

Satellite 
Pass  Time 

Latitude 

(°N) 

Longi  tude 

(°w) 

Case  Type 

Centerville 

Alabama  (AL) 

0346  Z 

30  Oct  79 

32.54 

87.15 

Clear 

Little  Rock 

Arkansas  (AR) 

0346  Z 

30  Oct  79 

34.44 

92.14 

Clear 

Greensboro 

North  Carolina  (NC) 

0251  Z 

23  Nov  79 

36.03 

79.57 

Clear 

Glasgow 

Montana  (MT) 

0435  Z 

23  Nov  79 

48.13 

106.37 

Clear 

Bismark 

North  Dakota  (ND) 

0350  Z 

30  Oct  79 

46.46 

100.45 

Cloudy 

Medford 

Oregon  (OR) 

0531  Z 

30  Oct  79 

42.22 

122.52 

Cloudy 

Green  Bay 

Wisconsin  (WI) 

0253  Z 

23  Nov  79 

44.29 

88.08 

Cloudy 

Monterrey 

Mexico  (f€X) 

0429  Z 

23  Nov  79 

25.52 

100.12 

Cloudy 

Dodge  City 

Kansas  (KS) 

0348  Z 

30  Oct  79 

37.46 

99.58 

Precipitating 

Omaha 

Mebraska  (NE) 

0349  Z 

30  Oct  79 

41.22 

96.01 

Precipitating 

Pittsburgh 
Pennsylvania  (PA) 

0252  Z 

23  Nov  79 

40.32 

80.14 

Precipitating 

Spokane 

Washington  (WA) 

0435  Z 

23  Nov  79 

47.38 

117.32 

Precipitating 

TEMPERATURE  (*K) 


!  36.  Comparisons  of  the  retrieved  temperature  profiles 
(dots)  with  the  radiosonde  data  (solid  lines)  for 
the  four  cloudy  cases. 
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Generally,  the  retrieved  and  observed  (radiosonde)  temperatures  in  clear 
and  cloudy  conditions  are  within  about  5°  K. 

In  the  final  figure  (Figure  37)  we  show  the  retrieved  tempera¬ 
ture  profiles  under  precipitating  conditions.  Again,  four  cases  are 
presented  in  this  study.  Two  cases  are  selected  from  30  October  1979; 
both  have  a  5  mm/hr  rainfall  rate  with  30%  cloud  cover  in  the  field  of 
view  of  the  SSM/T.  In  the  other  two  cases,  selected  from  23  November 
1979,  both  indicate  a  1  mm/hr  rainfall  rate  but  with  cloud  covers  vary¬ 
ing  from  50%  to  80%.  The  most  distinct  feature  in  the  retrieved  tem¬ 
perature  profiles  using  the  statistical  covariance  method  for  precipi¬ 
tating  cases  is  the  significant  and  consistent  deviation  from  the  radio¬ 
sonde  data  in  the  lower  boundary  layer  where  precipitation  takes  place. 
In  the  moderate  5  m/hr  rainfall  rate  cases,  the  differences  between 
the  retrieved  and  radiosonde  temperature  profiles  near  the  surface  are 
as  large  as  10  -  15°  K.  It  should  be  noted  that  precipitation  in  these 
two  cases  cover  only  about  30%  within  the  field  of  view  of  SSM/T.  As 
for  the  cases  involving  1  mm/hr  rainfall  rate,  about  5  -  10°  K  differ¬ 
ences  near  the  ground  are  observed.  The  findings  for  these  precipita¬ 
tion  cases  using  the  real  SSM/T  data  are  in  general  agreement  with 
those  described  previously  in  the  hypothetical  temperature  retrieval 
exercises.  As  shown  in  Figure  34  the  surface  temperatures  suffer 
increased  degradation  from  10°  to  20°  K  as  the  rainfall  rate 
increases  from  1  mm/hr  to  20  mm/hr.  These  hypothetical  analyses,  of 
course,  assume  that  the  cloud  covers  the  entire  field  of  view. 
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SECTION  6 


CONCLUSION 

In  this  study,  we  have  developed  a  microwave  radiative 
transfer  program  for  cloudy  atmospheres  applicable  to  the  DMSP 
SSM/T  channels.  The  transfer  program  takes  into  account  the  simulta¬ 
neous  effects  of  multiple  scattering  and  absorption  by  hydrometeors 
and  absorption  due  to  molecular  oxygen  and  water  vapor  in  the  atmos¬ 
phere. 

Sensitivity  analyses  were  carried  out  to  investigate  the 
effects  of  the  rainfall  rate,  cloud  thickness,  and  cloud  location  on 
the  upwelling  brightness  temperature  over  land  and  ocean  surfaces  for 
two  atmospheric  profiles.  The  effects  of  precipitation  on  SSM/T 
channels  1  and  2  are  shown  to  be  significant  depending  on  the  cloud 
liquid  water  content  (or  rainfall  rate),  thickness,  and  surface 
emissivity.  Over  the  land  surface,  increasing  the  cloud  liquid  water 
content  and  thickness  reduces  the  upwelling  brightness  temperature  for 
channels  1  and  2.  For  channels  3  and  4,  unless  high  rainfall  rates  are 
involved,  the  reduction  in  the  brightness  temperature  is  normally  in¬ 
significant.  Over  the  ocean  surface,  however,  the  increase  of  the 
liquid  water  content  and  thickness  may  increase  or  decrease  the  up¬ 
welling  brightness  temperature  relative  to  the  clear  column  value. 

A  significant  variation  for  the  brightness  temperature  is  shown  for 
channel  1.  This  conclusion  is  in  general  agreement  with  that 


described  by  Wilheit  et  al.(1977)  using  a  simple  scattering  program 
for  a  frequency  of  37  GHz.  Moreover,  sensitivity  analyses  also  reveal 
the  importance  of  the  position  of  the  cloud  layer  in  the  atmosphere 
and  the  atmospheric  temperature  profile  on  the  upwelling  brightness 
temperature  values. 

In  addition,  investigation  of  the  effects  of  precipitation 
on  the  temperature  profile  retrieval  using  both  the  theoretically 
simulated  values  and  real  data  was  carried  out  utilizing  the  seven 
SSM/T  channels.  The  retrieval  method  adopted  in  this  study  is  the 
statistical  method  developed  at  the  Air  Force  Global  Weather  Central 
in  which  the  surface  effect  is  removed  in  the  recovery  program.  The 
hypothetical  retrieval  exercises  show  that  the  temperatures  close 
to  the  surface  suffer  increased  degradation  as  the  rainfall  rate 
increases.  This  finding  is  also  supported  by  the  analysis  employ¬ 
ing  the  real  SSM/T  data  for  a  number  of  case  studies  in  which  tempera¬ 
ture  profiles  from  radiosondes  are  available  for  comparison  Further¬ 
more,  the  latter  study  indicates  that  nonprecipitating  clouds  have  an 
Insignificant  affect  on  the  microwave  temperature  retrieval  with  accuracy 
generally  on  the  same  order  as  in  clear  conditions. 

Although  the  current  study  employs  only  four  precipitating 
cases  in  the  analysis  and  may  not  be  conclusive  in  view  of  the  limited 
sample  used,  it  appears  that  the  effect  of  precipitation  on  the  tem¬ 
perature  profile  retrieval  using  microwave  frequencies  is  substantial 
and  significant.  Of  course,  the  reliability  of  the  statistical  method 
for  the  temperature  profile  retrieval  in  clear  atmospheres  should  be 
examined  comprehensively  and  completely  utilizing  the  data  that  are 
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available  in  different  localities  and  seasons.  In  addition,  in  order 
to  derive  the  temperature  field  over  the  global  space,  further  studies 
concerning  the  influence  of  precipitating  clouds  on  the  temperature 
retrieval  in  the  microwave  region  seem  also  warranted. 
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APPENDIX 


LISTING  OF  MICROWAVE  TRANSFER  PROGRAM 


MAIN  PROGRAM 


PARAMETER  NUMCHN=6,NUMLVs40.NG=16.  NMOST=NUMLV+1c 

U IMENS I On  T (NMOST ), W ( NMOST ). T  AUM ( NMOST . NUMCHN ) , T  AUS  ( NMOST  . NUMCHN ) • 

•  EMlS(NUMcHN) #TANT (NUMCHN) 

COMMON  /ANGLE/  UM <NG) , A <NG ) .PI 
COMMON  /SFCEMS/  EMiS 
COMMON/FREQ/FMU (NUMCHN ) 

COMMON/HANOU/H  ( NMOST  )  *0  (NMOST  .NUMCHN) 

COMMON  /TAUVAL/  TAuM 
COMMON/TaNOM/T.N 

COMMQN/PkES/P (NMOST ) 

COMMON/NtWLEV/AOOHT (10) »NOMNEMf LVNUM 
REALwfl  SFC (2) 

data  sfc/*lano'.* water*/ 

data  KTau1*kTAU2/*TAU( ' . 'GHZ) */ 

CALL  PROFIL 

CALL  TRANMW <  T  AUM , T aUS , 0 . , 1 » NUMCHN . 1 , 1 . ) 

WRITE(6,900)  (KTaU1.FMU(  J  ) ,KTAU2,J=1 , NUMCHN) 

00  20  Isl.LVNUM 

WRITE.  (6 .901)  P(l)*T(l).W(I)*tTAUM(I,J)  rj=l, NUMCHN) 

20  CONTINUE 

REAU(5.903)  (EM1S ( 1 ). Isl, NUMCHN) 

60  CALL  ANTEMP(TAUM,TaUS,T,T<LVNUM>,EMIS»TANT,lVNUM,1, NUMCHN) 
WRITE(b,902)  SFC(H,EMIS.TaNT 
CALL  oUFF(P.T.w.H.U) 

900  F0KmaT<*1  P(MB)  T(K.)  W(Q/KGM> • ,2X.7(2X,A4»F5.2.A4) ./) 

901  FORMAT (2(2X,F6«l).2X.E6a3.5X.6(E10a5.5X) .E10.5) 

902  FORMAT UhO, AO. *...  EMISSIVlTT* .4X.6F15.3./. 1H0.12X. 

•  ’ANTENNA  T£MP,»,1X.7F15.3) 

903  FORMAT  (  ) 

STOP 

END 


SUBROUTINE  PROFIL 

c*****«»*»**»*. 

C  THIS  SHBROUTINE  is  to  READ  OBSERVED  PROFILES  (TEMPERATURE  k.  mixing 
C  RATIO  6/*G.  PRESSURE  Mq.  AND  HEIGHT  KM)  OF  ATMOSPHERE.  AND 
C  INTERPOLATE  THE  PROFILtS  TO  THE  DESIRED  PRESSURE  LEVELS  (40  STANDARD 
C  cEVELS). 

C  THE  LEVELS  Can  BE  AOUEu  UP  TO  10  LEVELS  OF  YOUR  CHOICE. 

€****•••*•*«*•«•*••**•••**••*»••*•*••*•••*••**•••*•******•**••*•*,»••» 


PARAMETER  NUMLVS40.  NM0ST=NUMLV+10 
DIMENSION  T ( NMOST ). W ( NMOST ). H ( NMOST ) . OBPR ( 40 ) 

•  »0BH(40) »ObTEMP(4o) .OBW(40) 

COMMON/pRES/P (NMOST ) 

COMMON/H aNDU/H  » u 
COMMON/TaNDW/T.W 

CoMMON/NlwLEV/ADDHT (10) .NUMNEW.LVNUM 

200  FORMAT  (  ) 

201  F0RM«T(1x'NL»0BS  To  high.  MAX  VALUE  EQUALS  40.*) 
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REAL  (5.200)  NLVObS 
IF  (4G-NL v VBS>202, 205.203 
202  PRINT  201 
RETURN  0 


203  IF  (NuVObS) 10b« 105,204 

204  N=NlVOBS 

REACH5.2Q0) (ObTEMPiI>»I=N»l»-l> 

Rc.AU(5,2oO)  (Ob«(I)  ,I=N,1»-1) 

READ(5,2uO)  (UbPR(I) »1=N.1»-1> 

READ(5»2oO) <OBH(l) ,1SN,1.-1) 

*KITE(6.301)  ObTEMR 
4RXTE (6, 301 )  Ob • 

»RlT£(6,301)  OBPR 
»RITE(6,301 )  OBh 

301  FORMAT (10F6.3) 

302  FORMAT  (lx.  *PKOFIL  USED  CLXMO  DATA  (MID  LAT  SPRING/FALL)  ABOVc. 

»  F7,2>'  Mb*) 

K-l 

00  100  1=1,NUMLV 
IF  (ObPR(K) .tS.0.0)  60  TO  100 

10  IF  (ObPRiM  .fcU.PU ) )  60  TO  90 
IF  (ObPR(N) ,lT ,P{ I) )  60  TO  101 
IF  (OBPRlA).tT.PiI))  60  TO  59 
XF  (K.EO.l)  60  TO  100 
DELI =OBTeMP ( k ) -OoTfcMP (K-l ) 


ObLbsObb (K ) -0B4 (K-l ) 

DLNP=ALOb(OePR«KJ/OBPR(K-l) ) 

T ( I J  =OELT /DlnP* ALOb ( P ( 1 ) /OBPR (K-l ) ) +0BTEMP ( K-l ) 

N ( X ) =UEL«/DU*P*AL0& (P ( X ) /ObPR (K-l ) ) +0BW (K-l ) 

CALL  H1NTHP(OBPR(Kj»OBPR(K-1>»P(I)»OBM(K)»OBH(K-1)*H(I)) 

60  TO  100 


69  K=K+1 

IF  (K.E0.2J  PRINT  302, P(X) 

GO  TO  10 

96  T(l)=OBTfcrtP(K) 

4(1)S0B4(K) 

rt(I>sOBH(K) 

99  K=K+l 

100  CONTINUE 

101  continue 

10b  READ (5,200)  NUMNEM 
NE»=NUMNe« 

LWNUM=NUML  V  +NUMNE ■ 

IF  (NUMNEW.EG.O)  60  TO  210 
HEAD ( 5 , 200 ) ( ADDhT ( 1 » , XSl , Ng*  > 
NTOPsNUMlV 
K-l 

00  206  IsNUMLV, 1,-1 
110  lF(M(I).Lt.AOOMT(K))  60  TO  206 


0=1+1 

00  175  L=NT0P,0,-1 


P(L+l)sP(LJ 
T (L+l JsT (L) 

H(L+1)SN(U 

H(L+l)sH(L> 

175  CONTINUE 

NTOPsNTOP+1 

h(J)=AODhT(K) 

CALL  NE»LVL(0) 

KSK+1 

1F(K.GT.NUMNl*0  60  TO  210 
60  TO  110 
206  CONTINUE 
210  CONTINUE 
RETURN 
END 


9 
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noon 


SUbAOUT  INt  F.INTKP  (SP10P.  SPFIOT  .  SPP . S2T..P . S2BOT ,  SZZ  > 


lblS  SUnhOUl IbE  IS  TO  INTERPOLATE  THE  OuEAVED  Htl&HTS  To  THE  HEIGhTS 
oF  uLSIltn  PRESSURE  LkYtHS. 


IMPLICIT  DOUBLE  PkElISIOU  (A-R.O-Z) 

DATA  TOi./ 1  .D-a/R At-Kb/fi .  37 122f  >3/N/ 11/ 

Fic.  Au  SP ToP .  SPhOT  .  SPP .  S2I <|P .  S?l  OT  ,  SZZ 

DIMENSION  ZmaIKa  (  SO  •  bU  )  tP  (5o)  .  ZMAT  ( 2SuO ) ,  COLF  t  SO ) 

r.iMtNSiui*  array (  id/ 1 .ou-i •  1 .2su-i » i.ao-1  . i .eo-i *  1  .ssn-i * 

»  1.7. 1-1. 1.5SD-1. 1 . 47u— 1 . 1  .  42I>-1 . 1.320-1*  1.210-1 .1.3150-1/ 
ifnifr=i 

1  2F.01=l)HLE<Sx.boT> 

2T0P=DHlE(SeT0P) 

PHOT  -OHLt  ( SPbOT ) 

P  I  OP— t)Hot  (SPTGP ) 

PP=DBLE(SPP) 

ZTSTi/HoT 

IF (iTST.bl.2T0P)  _TST=ZToP 

IF  1 1 .0Ub.LE.2TST)  60  TO  P 

CoNS= ( S2  roP-S/HOT ) / aL06 ( SPTOP/ SPROT ) 

S/Z-aLOG l SPP/ SPHOT I *C0NS+S2hC7 
Go  1  (i  94 

2  bK=<PTOP— PB0T>/<l.Cu/<2ToP+KALKM»-l.LiiJ/(ZR0T*RADKM)  ) 

RNUML  1 . UO-RADKM*  ( I  P-PhOT  )  /RK-R  ADKM/  ( ZrtOT  +R AUK  v ) 
RDEN=(PP-PhoT  )/kK  +  l.UO/(7HOT«-RAUKP) 

IF (IENTEH.Eu.1)  SZZ=SNGD ( RNU* /ROEN ) 

IF (2TST-1. 00)49. h.b 
S  LVLA-.I 

21  SI =f)FtLE(SZ2 1 

IF  (2lST.Gl.ti.nOJH  lVLA=LVLA+1 
IF (21ST.6T.4.75U11  LVLA=LVLA'H 
IF UIST.6T.4.25D1)  LVLA=LVUA*1 
IF  1/TST.oT. 3. 751)11  U.LA=(  VLAF1 
IF  1/1  ST. oT. 3.2501)  LVLA=LVLA41 
IF  (21  ST .GT. 2. 450 II  LVLA=LVUA-»1 
IF (21S1.GT.2.15U1)  lYLA=IVLAF1 
IF  UlsT.faT.l. ASCII)  lVLA=LVLhF1 
IF (2TST.OT.1.45U1I  lVLA=LVLA+1 
IF (Z15T.&T.1.25U1)  LVLA=LVLA«-i 
IF (21  ST .oT . b . SOU  1  LYLA=DuLA+l 
IF  (21ST.GT.0. 1(00 )  LvLASLvLAfI 
A=AaraY(lYLA> 

2 1 NCS ( 2TOP-2ROT ) /F  Lo AT ( N- 1 > 

[10  10  1  =  1. N 
2MATRX (1.11=1. 

io  Continue 

DO  20  1=1 > N 

2MA1RX  <1.21=1./ (ZuOT+FLOaT (1-1 )*ZINC) 

UC  IS  J=3*N 

/F.ATh  X  ( I  .  J )  =/MATR A ( 1 «  J-l )  *2MATRX  ( I  >2 ) 

IS  CoNl 1NUE 

22=1 ./2MATRX (1.2) 

6k Ao=7Z+Ft  ADKM 

F ( I )=0EXP (-A»22) / (&AAD*GRAO) 

COEF ( I ) =F ( I > 

20  C0NU.4UE 
k=0 

OO  23  1=1. N 
DO  21  o=l. N 
K=k-M 

2F.A1  ( K ) S2MATRX  (O .  i  ) 

21  CONTINUE 
23  coni inue 

CALo  SIMo(2MAT .COLF . N. KS I 

SOMbO  T=CoF  F  ( 1  )  *ZBlT+COEF  ( 2  >  *r>LOfa  (  2B0T  I 

SUM  I  uP=COEF  ( 1 )  *2T^P»lOEF  ( 2 1  *f  iLOG  ( 2T0P ) 

FU)T=1  ./ZbOT 
ToP=l .//TOP 
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No  no  US5.N 

toE*- 1  Jls-tOcF  ( J)/FU*AT(J-2» 

SUMmoTSSUMROT+CNEF  (o)*BOT 
aci=eoi/ZHOt 

SU*»TOP=SUK  ToH  +  CoEf  ( J)  *TOP 
TCPslOP/ZTOF 
6 A  CONilNUE 

AlNlsSUp  foP-SUMFiCI 
nn=lPHUf-PTOF)/AIi.T 
60  10  tnl.96» .ItNIEn 

61  CONSS  (PmOT"PP)  /WK+SllF'FiOT 
Ft.sCONb-SlJMHOT 
Flsteivb-aUMIOP 

62  ZF'E»O=(2H0T-*210P>/2. 

IF  U«S<ZTOP-2B0T)-TOL)95, 95,65 

65  SOBtAuSCGEFU  )*2KLAN+C0£F(2>*ljL0G<ZMtMN> 

EAN= J ./2MEAN 

f)0  66  J=A,N 

S,UMEAFI=SUF.£AI'.*COEF  (J)*EAw 
LKN-i  AN/ZMEaN 

66  CONI INUt 

FW£ AMZCoNS“60MEAN 
IF  <FR*Flwt.AWl7(i,  95.60 
70  ZTOPsZMtAu 
F  1st- MEAN 
faO  10  o2 
60  2ilOI=2Mt  «K 
FR=FnFAN 
60  1 0  62 

9b  SZZsSlJ&L  IZMtAU) 

96  Z/SLtnLElbZZl 

SoMloPsCotF  1 1  )  *72-*CoEF  1 2 1  *DtOG  ( 22 ) 

E22SI./22 
NO  97  JSJ.N 

Sl.r/JoM=SUF-TOP-*CoeF  ( yj )  *EZ7 
EZ2=tZ2/22 

97  CONI  1 NUt 
AlNl=SU«TGP-SnMnOT 
PF=-6K*aINT4PH0T 
SPP=S.*6clPP) 

60  10  99 

EhTaY  PInTRPISP101-,SPF»0T,SPP.SZT0P,SZ60T,SZ2> 
Ii  NlfWS? 

6i>  IO  1 
99  CONUNUE 
At-TuAN 
ENO 
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SUBROUTINE  S1MQ  (A.0.N.KS) 


obtains  solution  lf  a  set  of  simultaneous  linear  eolations  ... 
a  -  matrix  of  coefficients  stored  columnwise.  these  are 

DESTROYED  111  THE  COMMUTATIONS.  THE  SIZE  OF  MATRIX  *A* 

IS  N  X  N. 

B  -  VECTOR  OF  ORIGINAL  CONSTANTS  (LENGTH  N)  WHICH  IS  RFPLACFO 
HY  FINAL  SOLUTION  VALUES.  VECTOR  X. 

N  -  NO.  OF  EQUATIONS  AND  VARIABLES. 

KS  -  OUTPUT  DIGIT 

(1  FOk  NORMAL  SOLUTION 
1  FOR  m  SINGULAR  SET  OF  EQUATIONS 


DOUBLF  PRECISION  A.H.BIbA.TOL. SAVE 
UlMtNSION  A(] ).R(I) 

TOLSU.ODO 

KSSO 

Oj=-N 

DO  81)  JS1.N 
JY=U*1 
k>u- UU+N+ 1 

OIGasO.uIiO 
1 1  — 

DO  2 0  lsu.N 
Io=iT+l 

IF  (NABS (blbA I  — DAnS ( A ( IU) ) )  10.20.20 
10  H1Gh=A(1J) 

IMAXSl 
20  CONI  1 HUE 

IF  (DAHS(BlGA)-TOL)  30. 30.40 
30  KSS1 
RETURN 

AO  I1=U4N*( J-2) 

I1=IMAX-J 
DO  SO  K=J,N 
11S11+N 
12SI1+IT 
SAVtiA(H) 

AllDSAd?) 
a ( 12 ) —Save 
SO  A ( 1 1 ) SA ( 1 1 ) /RIGA 
SA Vt— h ( IMAX ) 
hllMAXtSH(J) 

MJ)=SAvE/BlGA 
IF  (U-N)  b0.9o.80 
80  IuS=N*(J-l) 

UO  80  IX=JV.N 

lxjsias+ix 

ITSJ-IX 
DO  70  JXSJY.N 
1XJx=N*(JX-1)+IX 
UJXSIXJX+IT 

70  a(1xjx):a(IaJX)-<a(1XU)*a(JJX)  ) 

80  h( IX>=Hl IX)-(8( J)*A( 1XJI ) 

90  NYSN-I 
ITSN»N 

DU  100  u=l.NY 
IAS1T-J 
I8=u-J 
1C=N 

DC  1O0  KSl.u 
B( IB)=R(IR)-Al IA>*B11C) 

1 A— I .“N 
Kill  ltSiC-1 
RETURN 
Enu 
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bURHOUlluE  Nt*LVL(K> 

PmHmiwETEr  LtVS4G.  1nT=6.  (.MOST SLEV+lb 
C0MmoN/P«ES/PwE(NF0S1  ) 

CCiMmuN/HANOU/ ALT  (I.KuST)  .IUNMOST) 

CCM*oN/T AND*/ TEMP  ( NF.OST ) .  H20 (NMOST ) 
bPTuPSPRE (K*  1  ) 

SPrtoT  sPRt(K-l) 
bZTo.  SALTOUl ) 
bZRo  l=AuT(K-l  ) 
b/Z=«LTiK) 

TtSl=SZ8uT 

IF (TtST.bT.SZTOP)  TEST=b?TOP 

ipntsi.bE .i)  go  io  io 

KH0b= (SfTGP-SPRuT 1 / (SZTOP-SZROT > 

Pk£  (K ) =SPH0T ♦KHCfa* ( bZZ-SZHOT  > 
uu  In  11 

lb  C*LL  PlhTRPlSPTuP.SPI  OT.PkED1  > .SZT0P.SZH0T.SZ7) 

11  SFPafRE(a) 

TEMP(K)  A(TEr,P  (K  +  1  )-ltMP(K-l)  >*AL06(SPP/SPBOT)/ALO<;(SPTOP/SPBOT) 

•  ♦TEMP(K-l) 

HiO  <*.  IKH?0  tK4 1 1  “U2o (K-  1 ) ) * ALOG 1 SPP/SPBOT  )  /ALOG (SPToP/SPBOT ) 

*  ♦H*o(K-l) 

RETURN 

ENO 


SUftROUT  ll.t  I K  aNMW  ( T  a  W .  TAWSQ .  7 .  NU 1 .  NU2 .  KOUNT •  UMP ) 


COMPUTE  fHANSMlSSIVlTltb  TO  THE  TOP  OF  THE  ATMOSPHERE  FrO"  EACH  LEVEL 
OF  ThL  PROFILE  FOR  EACH  CHANNEL  (ANTENNA  GAIN  CHARACTERISTICS  INCLUDED) 


PaRmk.E TEr  NoFCHN=f..NOMLV=40.Nb=16»  NPoSTsNUMLV+10 
(JOUoLE  PRfClSlON  I  An. TAX. TO*. T1  .T2 

COMMON/ INPUT /ALT (NKOST) .  TEMP  INI'OST  ) »PrE INMOST) .H201NM0ST) . 

*  THanS(S)  .NIhICK.NiiASE.TRN(NFiOST.NUFCHN)  .N*L£V(10> 

COMmbN  /STOkAL/  ALFm 

COMMON  /TANUW/  T.W 

COMmON/hAU0U/H< NMOST ) .U(nMOST.NUMCHN) 

COMMOl./NEWLtV/ADnt.TdO)  .NUMNE  h.LVNUM 

Cilrttl.SION  T (NMOST) •  » INMOST) .ALFA (NMOST .I.UMCHN) . TAW (NMOST.NUMChN ) . 

•  TA*SG<N.  .CST  .NUIUCHN)  » GbUM (NUMCHN) 

DIMENSION  DCS! lib) . IO» ( 1?7> .SA (127) . bZ ( 127) »G ( 127.NUMCMN) 
COMF.oD/««1N/GANT  1  .27  »  NUMCHN) 

COMMON/PKKS/  P(NMoSI) 

Lf Vl-LVNOk 

NPSnDMlV 

ZLASTsSVR. 

01  Hs  •  (il  7453a 
NANb— 127 
ISW=1 

CKHSJ .17768 
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IF (IS*)  10*111.1 

1  I.1NISNAN6/2 
KI.=n1IJT 
hlBUItl.l 
LA=1. 

ABO. 

00  2  I=l.NlNT 
A=A*1 . 

Chzcost  A*C'Th) 

IiLBCm-CH 
,)!/=.  S.UC 
DCS  ( 1 4>iN  I  =00 

ucsiMM-i i=oo 

2  casco 

DO  7  KK=NUl  .NU2 
no  o  I*tt4il27 

3  01  I  >kK  )=i,ANT  (1  .KK  >*i,AM  (  1  .KK) 

no  4  1=1.63 

4  0(I.KK>=6(12H-1.KK) 

6S=0  • 

01=61 l.KK) 

00  fe  1=2.127 

02=6 1 1 »KK ) 

6S=GS+ (61*62 ) *DCS ( 1-1 > 

6  61=62 

OSUk ( KK ) =6S 
7  Coni inue 

SMAASASIn ( 1 . /CRH ) /DIR 
1SW=0 

10  IF (2 .FO.ZLAST )  60  To  20 
OO  1?  K=l » 127 
UOSFlOATIK-F,*) 
sa<k>=ahs  12-Fur.) 

S2(k)=1. 

If (S« (K) .GT.SmAA)  60  TO  12 
A=ChH«SIMDTR»SA(K)  1 
CUSASINI A) 

S2tM  =  l./COStOU> 

12  CONTINUE 
2LAS1 =2 

?U  IF (2.LE.K0UNT )  60  To  30 
CaLL  RAsEU(LVNUK) 

3l)  NO  120  JSNU1.NU2 
T AUwOsO. 
no  no  i=i. LEvi 
unsud.ol/UMO 
loo  iF(nu)  1oi.101.10*. 
ini  TaUsi. 

60  lo  103 

102  IF  (nu.GT .  .flSt  1  )  GO  TO  53 
00  1 O  34 

53  OU=.onEl+TAONN 
TAUNO=TAoNUf.2 

54  TAU=tXP(-l/U) 

103  00  104  6=1.127 
TAX=T AU 

IF(Sa(K) .GT.SMAX)  Go  TO  1 04 
if  Max. lt. .uoai)  oo  To  me 

IF  (1AX.1.1 .  .R999)  60  TO  1()4 
IF  IS* <K) .EQ.l. )  60  TO  104 
TaX=T AX**S2(K) 

104  TOW  IK )=TAX 
65=0  • 

G1S61 1 . J) 

TlsTow(l) 
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Uvj  1  Uf>  K=2«127 
b2=to(K»J> 

T2=1om(K) 

IF (Sm (K) .GT.SmAX )  GO  TO  1 05 
bb=6S+(Gl»Tl+b2*T2)*DCS(K-l» 
10b  61=G2 
lOo  T1=T2 

TAWl 1 • Jl=GS/faSUM( J) 
TAU2=T4w(I«J) 

T  AWR- .  t>*  ( T  AU1+T  AU2 ) 

T  uMSO ( I > J)=TAwR»TAWh 
TAU1=TAU2 
110  CONTINUE 
120  CONTINUE 
HETUKN 
END 


SUBROUTINE  DASEU (LEVI ) 

PARAMSTER  NUMCHN=A»NUMLV=40.  NMOST =NUML V+ 1 0 
DIMENSION  KERNEL (NMOST, NUMCHN > .GAMT(NMOST) .RHOW(NMOST) 

REAL  KERNEL 

COMMON /FRrQ/FRGHZ< NUMCHN) 

COMMON /HAN DU/M'NMOST) > U < NMOST r NUMCHN > 

COMMON/PRTS/PfNNOST  ) 

COMMON /T  ANDW/T ' NMOST ) ,W(NHOST> 

COMMON /TRNSM7 /TRANS <  NMOST , NUMCHN ) 

COMMON/ INPUT' ALT (NMOST) >  TEMP ( NMOST ) rPREt NMOST) .H20< NMOST) t 
*  THKNS  <  9  > . NTH ICK ,NBASE ,TRN( NMOST  f NUMCHN) ,NEULEV< 10) 

REAL  NBASE 

DATA  RSUBV/4 . 4 1 EE6/RSUDD/2 . 87E6/ 

C 

r:  P  =  PRESSURE  (MD) 

C  T  =  TEMPERATURE  (DEG  K.) 

C  W  =  MIXING  RATIO  ( GM/KGM ) 

C  H  =  HEIGHT  ( KM ) 

C  RHOW  =  WATER  VAPOR  DENSITY  (GH/M**3> 

C  KERNEL  =  WEIGHTING  FUNCTION  (TRANSMITTANCE/KM) 

C  GAMT  =  ATTENUATION  IN  LAYER  I  TO  1  +  1  (DB/KM) 

C  FRGHZ  =  CHANNEL  FREQUENCIES  IN  GHZ 
C 

PRINT  16f (FRGHZ(I)»I=1. NUMCHN)  9  ECHO  CHECK 
PRINT  17.LEV1 ,P(LEV1  )fP< l)  9  ECHO  CHECK 
C 

C*»* *********** ******* **** ****** * ****************************** ******** 

C  COMPUTE  WATER  VAPOR  DENSITY  (RHOW)  FROM  PRES  r TEMP  AND  MIXING  RATIO. 

C* ********************************************************************** 

DO  30  1=1 .LEVI 

RSUBM*  < 1 . / ( 1 . FW< I)*l .E-3))*!W! I >*1 . E -3*RSUBV+RSUBD ) 

PCGS*P ( I ) *  1 . E3  9  DYNES /CM**2 

RHOM=PCGS/(RSUDM*T<I) )  9  MOIST  AIR  DENSITY  (GM/CM**3) 

RH0W<I)=RH0M*U(I)*l.E-3/(l.+W(I)*l.E-3>  9  H20  VAPOR  DENSITY  <GM/CM**3) 

RH0W(I)=RH0W(I)*1.E6  B  H20  VAPOR  DENSITY  (GM/N**3) 

30  CONTINUE 
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c 

C ****************************************************»****,**»*»»*»*»*** 

C  COMPUTE  ATTENUATION  COEFFICIENTS 

C****************************************************************t****** 

PRINT  18 

DO  35  I-LEV1.1.-1 

35  PRINT  19,P(I>.H(I>.T(I).RH0U(I) 

NUMLAY=LEV1-1  8  NUMBER  OF  LAYERS  IN  VERTICAL  PROFILE 

DO  500  J  =  1 » NUMCHN 

8*1.0  e  TRANSMITTANCE  AT  EARTHS  SURFACE. 

A=0  •  0 

DO  215  1=2. LEVI 
TRANS ( 1-1 > J)=B 
U(I-1 . J)*A 

PTORR=P(I >*760./1013.25  B  CONVERT  PRESSURE  TO  TORR  (MM  HG) 

GAM  1 “GAM 02 < FRGHZ < J ) . T ( I ) .PTORR )  0  ATTENU.  AT  PTORR.T(I)  BY  02  (DB/KM) 

GAM2=GAMH20(FRGHZ(J).T(I) . PTORR . RHOU ( I )  >  0  ATTENU.  AT  PTORR » T ( I > , RHO ( I > 

GAMT ( I )=GAM1+GAM2  B  TOTAL  ATTENUATION  BETWEEN  LEVEL  I  AND  1+1  (DB/KM) 

C  FOR  GROUND  UP  USE  DH=H ( I +1 ) -H( I ) 

C  FOR  TOP  DOWN  USE  DH*H< I > -H < 1+1 > 

DH=H (T-l)-H(I)  B  THICKNESS  OF  LAYER  (KM) 

IF(HC-l)  .EQ.  H(  I )  )  DH=H  ( I  - 1  >  -H  (  I  +  1  )  SFOR  PR  IC I P I  TAT  I NG  CASE 
C****t ************************************************************ ****** 

C  LET  ADB=GAMT( 1 )*DH  (FLUX  ATTENUATION  IN  DB) 

C  BUT  DB=10LOG10( I ( 1 >/I(2> )  WHERE  J(1).GT.I(2) 

C  THEREFOR  ADB  =  GAMT < 1 ) *DH= 10LOG 1 0 ( I ( 1 ) / 1 ( 2 ) ) 

C  ADB/10=LOG10( I ( 1 >/I (2) ) .  OR  10** ( ADB/1 0 > =1 ( 1 > /I ( 2 > 

C  ( ADB/10)LN( !0)=LN(I ( 1 ) /I (2) ) .  OR  ( LN ( 10 > / 1 0 ) *ADB=LN ( I ( 1 ) /I ( 2 ) ) 

C  SO  (.23025F>*GAMT(1)*DH=LN(I(1)/I(2))=-LN(I(2)/I(1))=A 
C  SINCE  TRANSMITTANCE.  TNU2= I ( 2 ) / 1 ( 1 ) 

C  THEREFOR  (. 230259 ) *GAMT ( 1 ) *DH=-LN ( TNU2 > =A .  AND  EXP(-A)=TNU2 
t 

C  SIMILARLY  FOR  LEVEL  3.  ADB=GAMT ( 2 > *DH 
C  AND  (.230259>*GAMT(2 ) *DH=-LN ( I ( 3 ) /I ( 2 ) ) 

C  AND  (  .23025<?)*GAMT(2)*DH+A  =  A  YIELDS 

C  -LN( I (3) /I (2) >-LM( I(2)/I( 1 > )=-LN( I (3)/I ( 1 >  >=A 
C  THEN  -A*LN(I(3)'Id>>=LN(TNU3>.  THUS  EXP(-Ai=TNU3 
C 

C  SIMILARLY  FOR  ANY  LEVEL  N.  EXP(-A)=TNUN 

C*********** *****].******** ********************************************** 

A=A(.23025*GAMT(I)*DH  B  -LN(TNUN) 

WATEN=(B-EXP(-A) >/DH  0  (DELTA  TNU ) / ( DELTA  HT) 

KERNEL  < I >  J ) =WATEN  B  WEIGHTING  FUNCTION 

B=EXP ( -A )  0  TRANSMITTANCE  AT  LOWER  LVL  OF  NXT  LAYER. 

215  CONTINUE 

TRANS(LEV1 . J)=B 
U ( LEV  1 »  J ) =A 
500  CONTINUE 
1000  CONTINUE 

16  FORMATdHO. 'CHANNEL  FREQS  <  GHZ )  ' .  ( 10F7.2.2X) ) 

17  FORMATdHO. 'PROFILE  CONTAINS  '.13.'  LEVELS  FROM  '.F7.2.'MB  TO'. 

*  F7.2. 'MB' ) 

18  FORMATdHO. 'ECHO  CHECK' ,T15, ' PRES ( MB ) ' . T35 . 'HEIGHT (KM) ' .T55. 

*  'TEMPdiEG  K ) '  .  T75 .  '  H20  VAPOR  DENS  I T  Y  ( GM/M**3  )  '  .  /  ) 

19  F0RMATdX.T15.r8.3fT35.F7. 3. T55 .FA.2.T75.E9.4) 

RETURN 

END 


FUNCTION  GAM02(FKEw,T.P) 

c«***»**«*. 

C  COMPUTES  OXYGEN  ATTENUATION  USING  MEEKS  AND  LlLLEY  PARAMETERS 

c 

C  FREQ  =  FREQUENCY  (GHZ) 

C  11U  =  FREQUENCY  <hZ> 

C  I  =  TEMPERATURE  (DEGREES  KELVIN) 

CPs  PRESSURE  (TORR) 

C 

REAL  NU 

REAL  NUPc  (4S) #NUMI (45) 

IFdl.EQ.  10396)  GO  TO  2 

C 

C  OXYGEN  ABSORPTION  lines 

c 

NUPL(1)=56,2648E9 
NUPC(3>=58.4466e9 
NUPl ( 5 ) ==9 , 5910E9 
NUPCl7)=t>0.4348£9 
NUPL ( 9 ) sol , 1S08£9 
NUPc ( 11 ) =61 ,  8002E9 
NUPC(13)=62.41l2E9 
NUPci15)s82.9960E9 
NUPL ( 17 ) =63 , S685E9 
NUPc ( 19 ) =64 , I272E9 
NUPC(21)=64.o779£9 
NUPl  ( 23 )  =65 , 22**0e9 
NUPC(25)=65.7626E9 
NUPL127  )  s66 . 297 8E9 
NUPL(29)=66.83l3E9 
NUPL ( 3 1 ) =67 , 3627E9 
NUPL ( 33 ) =67 , 8923E9 
NUPL ( 35 ) s68 , 420SE9 
NUPc(37)=8fl,9478E9 
NUPc ( 3°) =69 ,4741£9 
NUPC ( 41 ) =70 . OOOOE9 
NUPL(43)=70.S244£9 
NUPL(4S)=71,0497E9 
NUM1 ( 1 ) =118, 7505E9 
NUMI(3)=o2,48fa3E9 
NUMI(5)=o0,306l£9 
nUM1(7)SS9.1642e9 
NUM1(9)=S6.3239e9 
NUMI ( 11 )=57,ol25E9 
NUMI(13)=S6.96S2E9 
NUMI(15)=56.3634E9 
nUMH17)=55.7839E9 
NUMI (19)=55,2214£9 
NUMI (21)=S4,6728£9 
NUMI 123)=S4,1294E9 
NUMI (25)=S3,S960E9 
NUMI i 27 ) =53 , 069Se9 
NUMI ( 29 ) =52 , S4  5fi£9 
NUMI|91>s52,02S9E9 
NUMI (33)=5l ,S09lE9 
NUMI(35)=50.9949e9 
NUM1(37)=S0.4830E9 
NUMI (39) =49 ,9730£9 
NUMI (41 )=49,464fit9 
NUMI (43) =48 ,9S82£9 
NUMI (45)=46,4530e9 
PM=267.4l 
11=10396 
2  CONTINUE 

NU=1,E9*FREQ 
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c 

C  COMPUTE  cINEmIOTH 

c 

u=.ki> 

IF(P.LT.PM)  bs •  25* ( 1 .  +  ( AL.06 ( PM/P ) ) /1 . 323 ) 
IF<P.lT.(18.957>)  0s.75 

UNU1  =1.9St+o*P  *<0.2l+0.7B*b)*<300,/T  )**.85 

UNU12=DNul**2 
C 

C  GAM02  =  ATTENUATION  COEFFICIENT  (DB/KM) 

C 

SUMsO.O 

UO  1103  N=l»Hb.2 

FN=N 

SUMsSUM 

l+(  tl./ I  (NUPulN)-NU) **2+0^12) 

2+l./( (NUPL(N)  +  NU) **2  ♦  DNU12 )  > 

3*<FN*(2.0*FN  ♦  3.0)/(FN  ♦  1.0)) 

4+U,/l(NuMI<N)  -  NU)**2  ♦  CNU12) 

5+1.  /((NUMI(N)  +  NU)**2  +  DNU12 ) ) 

6* (Fn  +  1,0 ) * (2 >0*Fn  “  l.OJ/FN 
7+l./(NU**2  +  0NU12J 

b*2.0* (FN**2  +  FN  +1.0)  * (2, 0*FN  +1.0>/(FN**2  +  FN)) 
9+E  *P  ( -2  .  Ot>84*FN*  <  FN+1 . 0 )  /T ) 

1103  CONTINUE 

GAM02  =2.67h2£-9*P/T**3+NU**2*SUM*DNU1 

RETURN 

ENO 


FUNCT  IOim  GAMHZO  <  FREto  .  T  .PTORR  .  RHO  ) 

C 

C  Ml  ~  FRtGUENCY  (HZ) 

C  1  =  TEMPERATURE  (DEGREES  KELVIN) 

C  PTo«R=  PRESSURE  (TORR) 

C  rHO  =  MATER  VAPOR  DENSITY  (G/K3) 

C  OAMH2U  i  ATTENUATION  COEFFICIENT  (DR/KM) 
C 

REAL  NU . nun 
DATA  NUN/22.235E9/ 

P=PTORR 

IF (FKEG.oT .6(1 . )  GO  TO  300 


COMPUTES  MATER  VAPOR  ATTENUATION  USING  BARRETT  AND  CHUNG  FORMULA. 


GARRETT  ANU  CHUNG  LINcMILTH. 

UNU*=  (P.62E+9+P/7h0.)/(T/31B.)»*.b2S+(1.0+.012  +RHO+T/P) 
COMPOTE  *ATER  VAPOR  ATTtNUATlON  COEFFICIENT  GAMH20. 

NUSFkEG* I «E9 

GAMhZO  =+.5faE-23*3.35Elb«RH0  +NU**2/T»*1 .5 
lMlEAPf-GMU./Tj/T+lOi^/UNU  -  NUN ) +*2  ♦  DNU2**2) 

2  ♦  UMI2/ ( (NU  ♦  NUN ) **2  ♦  DNU?»*2))  ♦  7.23RE-24+DNU2) 

RETURN 
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c 


an o  continue 


COMPUTES  WATER  VAPOR  ABSORPTION  FROM  6AUT  LINE  MOOEL  PON  183  6H2 
LIME. 


PMH=lQlJ.25/7r>0.*PTuRR 

CALL  ALPha5 ( FREQ . P Mh . T , RhO • NES • RESNON  > 

&AMhPO=RlSEKESNON 

RETURN 

END 


SUBROUTINE  ALPHAS  (FREQ.P.T.  KHOBAR# RES. RESNON) 

••••*. ***».**.*,*. •.*»*»«••«.*..***••**,*••»*•*, *•*•*••••••*«•** 

C  COMPUTES  ATTENUATION  FR0E1  IG.3  GH?  WATER  VAPOR  LINE  PLUS 
C  NONRESONANT  BACKGROUND 
C  ORIGINAL  PROGRAM  HY  N.  faAUT 

c»*»**** *•***»•*****»»***•***•**•••* A********************************* 

c 

C  FREQ  S  FREQUENCY  (GHZ) 

CP  s  PRESSURE  (MILLIBARS) 

C  1  =  TEMPERATURE  ( DEGRlES  KELVIN) 

C 

SCRTE (X)=SORT(X> 

EXPF  <X)stXP(X) 

CONv-lOud. 

FU4.i-323.7Stl 

F1G«S1H3.31U 

UtLUSU=(3.5hl •1(i.»*(-7)*FREG)**2*T  /18.0 

(.NU0u3=2. 79*  ( P/1013. <'S)*(300./T)**(0.t>19)*(l.  +  0.0066»RHOnAR) 

()NU0G3=SoRTK  ( I iELDSQ+LNU043*»?  ) 

DNUMvWSl . 04*0 .200* HHOBAK* (T/300. )•*((). 1) 

UnUwVG=SuRTE  (i)NUMv  W**2«UEL0SQ) 

UNUNW=U.RU*U.ol9*t  *(300./T) *»0 .62 
ONUNWSSuH TF ( DNUNW **  2 ♦ OELnSQ ) 

ONU1HGS2. h8*(P/10l3. 25)* (300. /T)** (0.70) *(1.+0.0064*Rh08AR) 
DNUlMiSSQRTK  <ONUlt>4**2*DEL0SO) 

Sisi./I <F164-FREQ>*a2*0NU164**2)*1./C (F164*FRE<J>»*2*DNU16U**2) 

SE=1 ./( (F093-FREQ>*»2*DNu093**2)tl./( { F093EFREQ ) **2*ONU093«*2 > 

GAM1 SO . 97m*RH0RAR*FkEQ**2/T** (2.5) *EXPF (“197, 3 /T) *DNU164*S1 
GAM?-0.ft49*HHoRAR»FRLQ**?/T**  (2.5)  *EXPF  ( -454 . /T ) *ONU093*S? 
GAMaSP,SS*FRKoA*2/ Flto4**3* (UNUWVUEONUnW )  *  <  300 ./T ) •*  C 1 ,S) /1 000 . *RHO 
lHAR*«.ln»»T«0.7G0/lul3.2S 
RtSsbAM 1 *CON V 
HESmoNS  ( GAM2EGAE.3  >  *CONV 
RETURN 

End 
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SUBROUTINE  ANTEMP  ( T  mO , T  AllSQ . T , TSFC . EM I S , T  A ,  ISFC  »  NUl  ,  NU2  > 


calculate  clear  column  brightness  values  for  each  channel  (antenna 

TEMPERATURES) . 


PARAMETER  NUMCHNSb.wOMLVsAO.  NMOSTSNUMLV-HO 

DIMENSION  T AU < NMOST • NUMChN ) , T  AUSQ ( NMOST . NUMCHN ) »T( NMOST ) • 

A  EM IS (NUMCHN ) *  T  A (NUMCHN ) 

DIMENSION  T AU 1 ( NUMCHN ) • TaU2 ( NUMCHN ) »F< NUMCHN I 
BIST (1) 

UO  3  J=NU) »NU? 

TAUSsTAU(ISFC.J) 

REFls).-EMIS(J) 

F ( J  >  STAUS*T aUSaREF L 
TAUl(J)sTAUU.J) 

3  T  A (o  >  SO . 

DO  S  IS2, ISFC 
B2=T ( I ) 

DO  4  JSNU1.NU2 

1AU2(J)STAU(1.J> 

TA(u)sTA(j)+0.5*(F,l+b?)*(TAUl ( JJ-TAU2 ( J) ) a ( 1 ,aF ( J) /TAUSG (I # J ) > 

4  TAU1 ( J)sTAU2( J)  ' 
b  B1SH2 

00  ft  JSNU1 »NU2 
HSstsFCaEmIS(J) 

TAUbsTAU(ISFC.J) 
ft  Ta(>j)sTa(J) aRS*T AUS 
RETURN 
END 


SUBROUTINE  BUFF (P*T . A.H.U) 

CA*»A*«a*»»*A««»AAA«AAAAAAAAAA*AAAAAAAAA«AA»*M>AAAAAAAAA««**«***«* 

C  INVERT  PROFILES  FOR  COMPATIBILITY  WITH  REMAINING  ROUTINES. 

C  MAKE  THE  OPTION  FOR  CLEAR.  CLOUDY,  OR  PRECIPITATING  ATMOSPHERES. 

CA*AA«*A«A**A*»AAA******4**AAA**AAAAAAAAAAAAAAAAAAA***AAA*A*AAAAAA 


parameter  LEVsAO.NUmCHNs*.  NMOSTsLEV+iO 

DIMENSION  P(NMOST) ,T (NMOST) .* (NMOST) »H (NMOST ) ,TRN( NmOST, NUMCHN) 
A  ,U( NMOST .NUMCHN) ,TAUABS( NMOST. NUMCHN) 

COMMON/ INPUT /alt ( NMOST ) . TEMP ( NMOST  >  »PRE ( NMOST ) , H20 ( NMOST  > , 

A  THKNS(9) .NThICK .NBASE.TRN.NEmLEV ( 10 ) .TAUABS 
COMMON/T  AUV AL/TR ANS ( NMOST , NUMCHN ) 

COMMON/SFCEMS/EMIS (NUMCHN) 

COMMON /NEALE V/AODHT (10) , NUMNEk .LVNUM 

JSLVNUMA1 

DO  100  lSl.LVNUM 

KSj-1 

PHEIK)SP(I) 

TEMP (K ) ST ( 1 ) 

H?OU)S«(I) 

ALT(K)SH(I) 

DO  SO  LSI .NUMCHN 
TRNIK.L)sTRANS(l.L) 

T  AUhHS ( K  «  L ) SU ( 1 , L ) 

SO  C0N11NUE 
100  CONTINUE 
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READ  IS.200)  ICHOSE 
60  10  (125. 150. 17b. 225) » ICHOSE 
12b  RETUHN 
150  CALL  DRIVER 
RETURN 

175  CALL  PF1LE 
RETURN 

22b  CALL  SlMPLR 
RETURN 

200  FORMAT!  ) 

ENO 


SUBROUTINE  driver 

C  rHlS*Is**H*ViNG*KOUTINL*FOR*CALCUl.ATlON  OF  BRIGHTNESS  TEMPERATURES 
C  WITH  INTERVENING  CLOUD  LATER. 

C  tXtCUTIWE  PROGRAM  For  CLOUO  and  PRECIPITATION— a  solution  of  the 
C  RADIATIVE  TRANSFER  EBUaTION  for  a  PLANE-PARALLEL  CLOUDY  ATMOSPHERE 

c  using  the  oiscrete-okoinate  method,  this  program  will  calculate  the 
c  radiance  < intensity j  above  and  below  non- isothermal  clouds  at  the 
C  hAVENUMBeRS  of  scams  Of  nimbus  6.  CLOUD  TEMPERATURE  gradients 
c  in  the  vertical  are  approximated  by  dividing  the  cloud  depth  into 
C  'N*  4.AYLNS.  each  layer  assigned  a  respective  TEMPERATURE  ACCORINip 
C  10  THE  INPUT  ATMOSPHERIC  profile,  a  system  OF  EQUATIONS  for  the 


C  RADIANCE  IS  UERIVEO  for  The  16  discrete  rays. 

c  •*•••****••••**••*•*••••*•••**••**•••**••**••*•••************»**•*' 

c 

c  *nput  parameters  (Unitsj 

c  ••atmospheric  conditions** 

c  aLI  altitude  OF  ATMOSPHERIC  later  above  the  surface  km 
C  temp  TEMPERATURE  CORRESPONDING  to  *alt*  deg  k 

C  PRE  PRESSURE  CORRESPONDING  TO  'ALT'  MILLIBARS 

C  H20  Water  VAPvR  CONCENTRATION  CORRESPONDING  TO  'ALT*  6M  cN-3 

C  LEV  number  of  atmospheric  profile  Layers  NONE 

c 

c  ••CLOUD  PARAMETERS** 

c  NTHICA  NUMBER  OF  CLOUD  THICKNESS  CASES  PER  RUN  (10  MaX)  NONE 
C  I HKNS  array  UF  LLOUO  THICKNESS  (NAX  of  9  VALUES)  KM 

C  NBASE  ALTITUDE  OF  CLOdD  BASE  (REAL  NUMBER)  KM 

C  NLYK  NUMBER  OF  LATERS  WITHIN  THE  CLOUD  (DEFAULTS 3)  NONE 
C 

C  ••EXTINCTION  PROPERTIES** 

C  PIN1  LEGANORE  POLYNOMIAL  EXPANSION  OF  PHASE  FUNCTION  NONE 

c  pt  single  scattering  albedo  according  to  *nu»  none 

C  NU  WAVE  NUMBbR  (LIMITED  TO  'INT'sTj  CM-1 

C  oEXT  EXTINCTION  COEFFICIENTS  ACCORDING  TO  'Nu*  NONE 

C  LV  LINE  contribution  to  absorption  coefficient  none 
C  CK1  CONTINUUM  contribution  to  absorption  COEPF  taken  none 
C  CR2  From  PAPER  BY  K  J  bIGNELL  (X9T0)  none 

C 
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KfcAL  NU.lV'NbASC 

uOUoLt  PkECISIOn  PiNO.UM.A.Z 

parameter  lev=<*o.  ints6.  mzszs,  ng=I6,  ng2=s.  nhostscev+io 
COMMON/ InPUT/ALT <  NmOST ) , TEMP ( NMOST )»PR£( NMOST ) • H20 ( NMOST t  $ 

•  THKNS (9) • NT HICK .NfaASE # TRN (NMOST , I NT  > .NWLEV ( 10) t 

•  TAUAdS (NMOST .INT) 

COMMON  /ANGLE/  UM(N6).A(NG)«PI 

COMMON/U>MUlT/U0C  *  aTMD*  ATMR , ATMA.SFC . THRUSV ( XnT ) , ThROPT ( I NT J *  ITER 

•  »ECb»IPcHAN(INT) »CLOEM 
COMMON  /ROLY/  PXNO(NG) 

COMMON  /kINOOM/  UP(NG2'lNT),D*(NG2'lNT)'THeTA(NC)'NUCXNT)»LV(XNT)> 

1CK1 < INT) . CK2  < INT  > . 1UVERT ( INT) ’ COVERT ( XNT  > 

COMMON  /MAO AT A/  PlNl(NG.INT) »PT(INT)*BEXT(InT) 

COMMON/NE«LEV/A0OH1 UO) .numnew.lvnum 
COMMON  /FRE0/FGHZ<1NT) 

COMMON  Z(1300) 

COMMON/Fl A6S/ XFL AS  k  2 ) 

DIMENSION  N£«LV<10) 

OATA  SPOlST /2 • 997929E 10/ 

9  FORMAT (  ) 

10  REAL  9 .LAST 

IF  (LAST.NE.O)  SO  TO  too 

C  ••*•••*«••**  INPUT  ClOuO  PARAMETERS  ••*•**.«»***•*«•• 

READ  9.  (IPCHAN(I).lsi.XNT) 

READ  9.  uBC » ATMp . AT MK * ATMA . SPC » ECB » CEDEM » X TER 

READ  9.NTHICK.NBASE.NLTR 

NTHSnTHXcK 

READ  9»(THKNS(X).Isl.NTH) 

READ  9» (FT { I ) » 1st. INT ) 

READ  9.(bEXT(I).Ist.XNT) 

DO  12  Jst.INT 

12  READ  (5,1900)  (PlNt(I,J),Xst.IS) 

1900  FORMAT (SFlS.S) 

READ(S«9) (IFlAS(l) .jst.2) 

lFLASU)s  0  DO  MOT  INCLUDE  SFC  REFLECTION  OF  BRIGHTNESS  CONTRIBUTION 

from  AboVE  scattering  layer  pith  lower  bounoart  value. 

t  INCLUDE  SFC  REFLECTION  OF  BRIGHTNESS  CONTRIBUTION 

from  above  scattering  layer  with  lower  boundary  value. 

(  ASSUME  TOP/DOWN  THROUGHPUT  EttUALS  1.) 

lFLAG(2)s  0  oo  NOT  WRITE  SCATTERING  PARAMETERS  TO  FILE  20. 

.Nt.O  WRITE  SCATTERING  PARAMETERS  TO  FILE  20. 

(IFLAG(2)sRAINFALL  RATe  IN  Mm/HR  100  YILLDS  RFRsO 

15  NU(I?s(FirlZU)*l,E9)/SP0L0T  fl  CHANNEL  WaVE  NUMBERS. 

T^STsNBASE 
DO  50  Ist.LVNUM 

IF  (ABSCAlT (D-AOOHT  (KK)  )«l,E“S)3b»35.30 

35  HEWLV (KK)Sl 
KKSKA+I 

36  IF(AbS(AlT(I)-TeST)-I,E-6)90.90.50 
«*g  NmLEV (K)sl 

KsK«l 

TESTsNB ASE+THKNS ( K-l ) 
if(kk.gt.numnew)  gd  to  go 

50  CONTINUE 
GO  PRINT  2000 
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2000  FORMaTUhO.*THE  LEVELS  ADOEO  TO  THE  PROFILE  ARC  -  »*/.SX, 

•  'PRESS (MB) • rlOx> ’nElGHT(KM) ' »10X» 'TEHP(DEO  K)',10X» 

*  'Mix  Rat (om/kom) *  > 

DO  til  XsltNUMNCH 
JsNEKtVU) 

PRINT  200l.PHE<  J)  .ALT (J)  .TENPC J>  «H20(J) 

2001  FORMAT<5x»F9,2rlOX,F10.2.10X,F11.2,10X,El5.5> 

81  CONTINUE 

KSTOPsO 

C 

C  CHECK  FOR  EAST  DATA  CASE 

C  DETERMINE  EXTINCTION  COEFFICIENTS  ONCE  per  data  case  usino  >kSTOP* 
C  IN  SUBSEQUENT  ROUTINE. 

c  execute  entire  code  once  for  each  input  clouo  thickness  (Thkns> 
c 

DO  90  KTmKsI.NThICK 
as  call  bnoky  <kstop»kthk.nlyr> 
call  RAO  (KTHK.NLYKfKST0P.S6S) 

90  CONTINUE 
100  CONTINUE 

return 

END 
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SUBROUTINE  BNORr(KSTOP,KTHK.NLYR) 


CALCULATE  upward  AND  DoWNWARO  brightness  TEMPERATURES  at  the  cloud 
BOTTOM  ANU  Top  (BOUNDARY  CONDITIONS) 


REAL  LVrNU 

PARAMETER  L£V=G0#  1NT=6,  KZs25,  N6=16.  NG2=0,  NMoSTsLEV^lO 
UOUbLE  PRECISION  Um,A,CCRAO 
REAL  NBASE 

COMMON/LbMULT /UBC » mTMO *  ATMR , aTMA . SFC . THRUSV ( INT ) . THRUPT (INT) *  INTER 

•  . ECo » IPCHAN (INT) * CLOEM 

COMMON/ INPUT/ ALT (NmOST) , TEMP (NMOST) ,PRE< NMOST) .H20< NMOST) » 

*  THKNS{9) »NThICKfNbASEf TRN(NMOST.INT) »NWLEV(10> 

*  . TAUABS (NMOST* INT | 

COMMON  /ANGLE/  UM(NG) »A(NG) .PI 

COMMON  /wINDOW/  UP<NG2,INT>,DW(NG2,INT),THeTA(N6),NU(lNT),LVlINT). 

•  CKi { INT) • CK2  < INT ) , CUVERT  < INT ) « COVERT ( INT) 

COMMUN/Nt*LE»/AOOHT ( 10 ) , NUMNE* *LVNUM 
COMMON  /SFCENS/  EMlSfNB 

common  ZNADaTA/  PIM (N6« INT) ,PT(INT) *BEXT(INT) 

COMMON/FLAGS/ IFLAG (2 ) 

DIMENSION  TaUSO( NMoST, INT) 

DIMENSION  TaUM I NMOST  * INT ) * TAU ( NMOST # INT ) 

o i Mansion  t aumoo ( nmost , i nt  > , em i s ( int ) , eup ( int  > , eow ( int ) 

DIMENSION  CUH0RZ (InT) , CDHORZ l I NT ) , CCRAD ( I NT ,N6 ) 

PI=S.1G1B926&S6 
KSTOPsKSTOPwl 
IF  (KSTOP.GT.l)  go  to  go 
write  (6,160)  N8ASE * ThKNS (KTHK ) 

WRITE  (6,160) 

DO  SO  1=1* INT 

IF  (IPCHAN(I).EO.O)  GO  TO  SO 

WRITE  (6,190)  (NU(i)'PT(1)'B£XT(1),(PXNI(J,X),Jsi,NG>> 

SO  CONTINUE 
GO  CONTINUE 

DETERMINE  GAUSSIAN  QUADRATURE  WEIGHTS  FOR  DISCRETE  ORDINATE  METHOD 
NGsNUMBER  OF  GAUSS  POINTS  (16) 

UMsCOSlNE( THETA) 

AsdUADRATURE  weight  VAlULS 

CALL  GAUSS  (NG*-1. 000,1.000) 

LOOP  FOR  ANGLES 

K0UNTS2 

DO  102  Jsl'NGE 
UU=U*NS2 
UMOsSNGL(UM(U) > 

UM(DD)=-UM(U) 

A(UJ)=A(J) 

UM1=SNGL(UM( J) ) 

UM2sSNGL(UM(uU)) 

THETA(U)sACOS(UM1)»160./PX 

ThET A ( UU ) SACOS ( UM2 ) • ISO , /P I 

CALL  TR ANMW (TAUM, TaUSG * 0 . *  1 » INT.KOUNT , UNO ) 

DdUsLVNUM+1 
DO  70  Ksl'LVNUM 
LsUudoK 
DO  60  M=1*INT 
60  TAU(L»M)sTauM(R»M) 

70  CONTINUE 
N6=NWLEV(1) 
ntopsnwlev(kthk^i) 


c 

C  LOOP  FOR  CHANNELS 
C 

00  SO  Xzl’lNT 
THRUPT(X»Zl.O 

XFUPCHAN(I).EQ.O)  60  TO  *0 
00  7S  MsuT  OP  t  L  VNUM 
IF (TAU(M.I)  ,NE.  0.)  60  TO  74 
TAUMOb(M,X)sO» 

60  TO  75 

74  TAURUS  <M, I ) sTAU(NToP>  X ) /TAU (M* X ) 

75  CONTINUE 
60  CONTINUE 

00  101  XslfXNT 
XFUPCHAn<X).£6.0)  60  TO  *01 

C  COMPUTE  0O»N*ARD  INTENSITY  At  THE  CLOUD  TOP 
C 

£D«(I)sO. 

NT0P1sNToP*1 
00  95  MsnTOPI.LVNUM 
TbARs ( TEMP ( M-l > wT EmP ( M ) ) * . 5 
£0*(l)s£ow(n+TBAR*<TAUM0D(M>l'X)-TAUH00(MtX>) 

95  CONTINUE 
U»  ( J* 1 >s£DW ( X ) *UBC 
CCKaU(I*oU)sUBLE(DwCJ*X)) 

00  97  K=1*LVNUM 
TAULOD(K« I )sTAU(l • i ) /TAU (K> I ) 

97  CONTINUE 

C  COMPUTE  UPNAHO  INTENSXlY  AT  BOTTOM  OF  CLOUO 
C 

bSFCsO • 

OATMsO. 

RATMOsO* 

KATMUsO* 

RCFLsI.-EMISUI 
lF(NB.EO.l)  60  TO  96 
DO  100  KsS'Nb 

T6AR= ( TEMP ( K ) wTEMP ( K-l ) } * .5 
lF(TAU(Nb> I)  .NE.  0.)  60  To  90 
UATMsDATM 
60  TO  96 

90  OATMsUATnwTbaK* ( T AU (A *  I ) -TaU (K-i 1 1 ) ) /T AU INB » X  > *ATMO 

96  RaTMd=RATMB«TBAR* ( T AUMOD (*-1 » I ) -T AUMOD <K *  X ) ) *REFL*T AUMOO ( N8 » 1 ) 

*  «aTM* 

100  CONTINUE 

96  bSFCsTEMP < 1 ) 'EM I S ( I ) «T AUMOD ( N6  *  X ) *SFC 

RATMUSEOW ( I ) *T AUMOO ( Ne . x ) 'REFL'TaUMOO ( Ne . X ) • ATMA*THRUPT ( x ) 
1F(ATMA»£6.0) 

•  RATMUSCDVEKT  < I ) *T AUMOD (NB 1 1 ) 'REFL'T AUMOD  ( NB . X )  *ECB 
XF  (IFLA6<1).E«.0)  RATMUsO, 

£UP  < i ) so aTMwR A  TMb  tRATMU 
UPU#l)StUP(l)*BSFC 
CCRmO 1 1 »  J ) =D6LE ( UP ( J  r I ) ) 

101  CONTINUE 

102  CONTINUE 

C  (.XTKaPOLaTE  INTENSITIES  To  THETAsO,  AND  90  DEO. 

C 

ICCKmDsI 

XEMISsO 

00  110  UaTsI.XNT 
1F(1PCHAn(UaY).E6.0>  60  TO  110 
JOUMsiNT 
XUPDwNsO 
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CALL  CXThPO (THETA , tCRAD.  TVgRT » JAY  r JOUM » X  UPOON | 

CALL  EXTRP9 < THETA  *  cCKAD , THOR  12 #  JAY , JOUM » I UPD*N > 

CUV£kT(JaY)sTV£RT 

CUH0R2<JmY)sTH0R1Z 

1UPU»I«1 

C ALL  CXTHPO ( THETA «  CCRAO  >  R vert • JA Y • JOUM  *  1 UPOBN ) 
CALL  EXTRP9 (T  HET A  * ttR AD , KHOR X2  »  JAY , JOUM , IUPO*N » 
COVERT ( JmY ) sKVERT 


CUH0R2( JaYIsHHORIZ 
Ufi  CONTINUE 

MK1TL  (6,200)  (ThKnS(KThk) , (THETA(J) »Js1»Ng2) ) 

DO  120  1=1* InT 
IPIlPCHAN(l).EO.O)  60  TO  120 

ftRXTE  (6,210)  (NUU)#CUVERT(I)*(UP(J*X),Js1,NG2).CUH0RZ<I)) 

120  CONTINUE 
N621sN62«1 

•RITE  (6,220)  (THETA(J),J=N621,N6) 

DO  ISO  1=1* INT 

IF  (iPCHAN(I).EO.O)  60  TO  130 

*RITt  (6,210)  (NU(i)#CDVERT{I>*(D*(J»X>,Jsl,N62),C0H0RZU>> 

130  CONTINUE 
RETURN 

160  FORMAT  t//,54X* 'CLoUD  STRUCTURE’ »/,22X, * CLOUD  TEMPERATURES  s  (DERX 
1VED  FROM  the  ATMOS.  TEMP.  PROFILE  ♦  CLOUD  THICKNESS) **/.22X, ‘CLOUD 
2  BASE  HEIGHT  S' ,F5.2, 'KM* ,/,28X, 'CLD  THKNS  =',F*.1,'KM  (MAX1MU 
SM  of  10  VALUES) ',//,3SX,« input  ABSORPTION  COEFF  », 

4/,*H»X*  'NO'  ,9A«  'LV  ,8X,  *CKl*  ,6X,  »CK2' ) 

ISO  FORMAT  (/,37A, 'SINGLE  ScaTTERINO  PROPERTIES  OF  BaTeR  DROPLETS* »/# 

111X,  «NU»  ,9X,  'PT'  ,7*,  'BEXT' ,5X*  •». . P 

2INI . . . . . . 

190  FORMAT  (9X*F6.A,2(HX»F6,3) *9X*0(F8,5*1X) ,/ ,39x,6(F6«S, IX) ) 

200  FORMAT (Im1,37X*»  BhIOHTNESS  TEMPERATURE  AT  THE  CLOUD  TOP  AND  •» 

•  'BASE' */,«BX, 'CLOuD  THICKNESS  ',F5.3»'  KM*,//49x*»THE  UPMARU  *. 

•  'BRIGHTNESS  ' *//Aa# • (THETA)  (0,0000) ' »SX,8( •(' .FT,** •)• ,3X) * 

•  • (90.0000) ',//,5X, 'HAVE  NO') 

210  FORMaT(5a,F6.«.F10.3,1X,9(F10.3,2X),F10.3) 

220  FORMAT  (/,A8X» 'THE  OOMNWARO  BRI«HTNESS'.//,AX,. (THETA)  (180.000)', 
•SX,8( • ( • ,F7.3« ' ) • »oX) , • (  90.000) *,//,5X, 'WAVE  NO') 

END 
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SUBROUTINE  GAUSS<N<XL<XUI 


SUBROUT INE  'GAUSS'  DETERMINES  THE  GAUSSIAN  QUADRATURE  VARIABLES 
<»R'='UM'=  COSINE IThETA) ,»W*S'A*s  WEIGHTING  FACTOR)  FOP  THE 
DISCRETE  ORLINATE  METHOD. 


IMPLICIT  DOUBLE  PRECISION  <A-H<0-2) 

COMMON  /ANGLE/  R(lfa)<W(lf.)<Pl 
DIMENSION  2<5o)  <PA(50) 

TOLsl. 00-16 

Pll=3.l415<}2f>53baV7<J3D4Gii 

AAS2.0b400/PH**2 

ABS-h? .00400/  l 3 . 0u4O0*PI 1**4) 

AC=151 In .00400/  ( 15.0U400*PI1**6 ) 

AUS-1 255*474. 0400/  (105. 0n+0U*PU**8) 

PAUUI.OD4UO 

ENSN 

NPl=N4l 

US  1 . 00400- ( 2 . 0D4  0O/PI 11**2 

U=1 .of)4«U/0SGRT<  <LN*0. 50*00 )**24U/4.0U400> 

DO  10  1=1 < N 
SL=1 

A2S4.0U400*SM-l.bu40b 

AESAA/A/ 

AFsaB/A2**3 

AGSAC/A2**s 

AHsaO/A2**7 

10  2< 1 )=0.?oD400*PIl*  < a/4A£*AF4AG4AH) 

DC  nO  K=1  <N 
XsOCoS ( 2 1 K  >  *D ) 

20  PA (2 1SX 

DO  30  NN=3<NP1 
ENNSNN-1 

30  PA(NN)S( <2.0f)*0o*ENW-1.0n400>*X*PA{NN-l)-tENN-I.0D400)*PA1NN-2>)/ 
tENN 

PNPsE  N* <  P A  <  N ) -X  *P A <  NP 1 1  1 / < 1 . 0D400-X*  X ) 

X1=A-PA(nP1)/PNP 
XIJSDaBS  (  X 1  —  X  1 
XOO=XD-TOL 
IF  ( xDU 1  50 <50 <40 
40  XSXI 

GO  TO  20 
50  KlKISX 

no  G(K)=2.0U400*U.bU4uO-X*x>/(EN*PA(N> )**2 
RETURN 
ENO 
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SUBROUTINE  EX  TRPO ( TnE  TA . INTERS .TVERT.UA 1  . JOUM . IUPDWN ) 


This  SUBROUTINE  IS  10  ExTRAPOLTATE  the  INTENSITIES  TO  THETA30'  an.)  90 
RY  APPLIN6  LEAST-SUUaRE  POLYNOMIAL. 


IMPLICIT  double  precision  CA-H.O-ZJ 
REAL  THETA *1 WERT 
DOUBLE  PRECISION  INTERS 
DIMENSION  TBETA(Io) •  INTEmS  (Jf)UM»  16) 

Kl=lw»*IUPO*N 
XXSFLOAT ( lUPDwNwlbO I 
S  K2=K1+1 

TISTNETA(KI) 

T*=THETA<K2> 

T3s1h£TA(K3> 

HISINTENS(JAV.KI) 

82SINTEwS<JAt.K2) 
nJ=INTENS( JAY.K3) 

C3s ) N2-B1 1 ✓ ( C  T2-T 1 ) * ( T2-T3  ) ) ♦ ( R 1-R3  > / C ( T3-T1 ) • ( T2-T3  > ) 

C2= I TR3-R1 >-C3» ( To-Tl >  * « T3+T1 ) > / ( T3-T1 ) 

C1=B1-C2*T1-C3*T1*T1 

TVEn1=C1+C2*XX«C3*Xa*XX 

60  To  SO 

1HE  Follow  I  nc>  Entry  IS  FOR  TbETAs90. 

ENTRY  ExTRPOl THETA. INTERS. TWERT . JAY. JDUM» IUPDMN) 

XX=9b. 

KISto-wRAlUPDRN 
60  TO  S 
50  CONTINUE 
RETURN 
END 
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SUBROUTINE  RAD  <aThK»NLYR»KSTOP,$) 

*♦»**»••*»*  •***•**■***•••♦*•**•••♦•*•• 

FOR  A  GIVEN  ctOUu  D£PTh  ( 'KTHK* )  THIS  ROUTINE  APPLIES  THE  DISCRETE 
ORdINATe  METhOO  TO  THE  RADIATIVE  TRANSFER  EQUATION  FOR  EACH  OF  THE 
SPECIFIED  SAVE  NUMBERS  (*NU»>. 


PARAMETER  LEVsOO,  inTs6,  M2s25,  NGsl6r  N62=8r  NMoSTsLEV+10 
UOUbLE  PRECISION  PlNO.UM. A 
COMMON  /AHGLt/  UM(NG) i A(NG) »PI 

COMHON/LbMULT/U8C#ATM0»ATMR»ATMAiSFC#THRUSvaNT).THRUPT(lNT).lTER 
•  *ECu> IPCHANt INT ) 

COMMON  /POLY/  PINO(NG) 

COMMON  /HADaTA/  PINI(N6iINT)*PTCINT>»BEXT(INT) 

LOOP  Thru  Main  CALCULATIONS  FOR  EACH  save  NUMBER  (INOEXS'UAYM 

00  30  JAYS1,INT 
XF(lPCHAM(UAY).EQ.b)  60  TO  30 
uO  10  Lsl'NG 

PIN0(L)S0BLE(<PlNIiL«UAY))«PT(JAY>) 

10  CONTINUE 

IF  (ASTOP.6T.lt  60  TO  20 
CALL  IGEn  (JAY) 

20  CALL  TGKaD  (OAY.kThK«NLYR> 

30  CONTINUE 

UO  40  1=1 » INT 

1F( IPCHANi 1 ) .EG.O)  60  TO  *0 

IF(AbS(THKUSV(I)-ThRUPT(I)).6T.0.ij  RETURN  4 
40  CONTINUE 
RETURN 
END 
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SUBROUTINE  lGfc'N  (wAY ) 


SlibkOuT  1  nE S  1  IGEN  •  •  »S>  AkT  •  #  *Fl LLH *  ,  •  FlLLC  •  .  •  F I LlD •  *  •  FILLS •  >  •  G*PRD • 
*LEP'  <  •UPHHM*.  AND  ' liPQFB'  DETERMINE  THE  EIGENVALUES  RT(I) 
ACCORDING  TO  THE  WAVE  NUMBER  (LIOU.  1973/ APPENDIX) 


PARAMETER  NG=16.  NG2=8»  lNTSf> 

DOUBLE  PRECISION  PlNO»Z. VALUES 
common  2(i30iu 

COMMON  /EIGEN/  VALUES* I NT.NG2) 
common  /POLY/  PINO(NG) 

NClM=130U 

NXSUO 

NSN6 

N0Im=3*n»<N*1 ) 

Ir  (MOlM.GT.NOIM)  GO  TO  lO 
KU2=(NUlM-N)/i! 

Nl.SN  **2 
nS=3*NN+1 

CALL  START  (/Ml) */(NN+l) «Z(2*NN+l>  »P1K0»Z(1) *Z(NS*2*N) >N>7(2*N+1) • 
12 ( 3»N+1 ) > JAY  > 

GO  To  2u 

10  MRITE  (n/30)  N/NOIM/MDIM 
PRINT  4U 
20  CONTINUE 
RETURN 

30  Format  (315) 

40  format  (*OPROGRAM  Mi  1ST  6F  REDEMENSIONED  AT  MOIm,> 

El/O 
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SuHKOUl lNfc  STmRT  (B.C.P.T.D.S.N. COMPlE  . POL •  uAY ) 


SUHHOUTInE  'START •  IS  To  COMPUTE  EIGENVALUES  ‘VALUES (JAY  *  I ) •  AND 
RETURNS  1  HE  EIGENVALUES  THROUGH  *  I GEN'  TO  'TGRAG'  WHERE  THE 
RADIANCE  CALCULATIONS  ARE  MAllE. 


PAKaaETEH  No2=8>  UG=16,  TNT=6 

UOU6LF.  PkECISIOIm  UM. A.P.C.B. S, D. T .OMU . COMPLE .POL . G .SGRPOT.PROOT . 
1 VALUE .VALUES 

COMa-uN  /ANGLE/  UM(Nt) .A(NG) ,PI 
COMMON  /EIGEN/  VAlUES { INT ,  NG2 ) 

U1 MENS l ON  B(N.N) ,C(n»N) ,P(N.N).T(N>  «D(N> ,S(N) 
dimension  COMPLE(N) .POL(N) ,6(N6) .DMU(30> . VAlUE <NG2) 

NLI=N-1 
Do  10  1=1. N 

10  CALL  LEP  <P(1.I). DM ( 1 ) . Nl 1 ) 

DO  20  1  =  1.  N 
DO  20  J=1 >N 

CALL  F1LLC  (T.P(1.1).P(1.J).A(J).CU.J).N) 

?o  Coni inue 

Call  fillr  (h.c.uw.n.nli i 
CALL  PILLS  IS.B.C.P.N) 

CALL  FIllO  (O.S.N) 

IF  (AHS(T(D-l.O).LE.l.E-lO)  D(N|rO. 

00  SO  1=1. N 

OMU  ( 1 )  =1 .DO/UM  ( 1 ) 

30  CONTINUE 
N2=n/2«1 
1HSN+2 
t(N2)=1.0G 
00  40  1=2. N. 2 
J=(iH-l>/2 

6(J)=0(I> 

40  CoNUNUE 

ChLL  DPR!!*  It. N2.VALUE. CoMPLE. POL. IR.IER) 

00  no  1=1. IK 

SuROuT=USORT ( UAHS ( VALUE  < | ) >> 

PkOOTsSuROOT 
DO  SO  J= 1 »  NL 1 

SO  PK0OT=  I PWOOT  ♦!)  ( J )  JaSGROOT 
PKOOT  SPKOOT  ♦  ()<  N) 

VALUE ( I ) =SOROOT 

IF  (PKOOT. 61.1. 00-02)  WRITE  (6.80)  PROOT 
60  CONUNUE 

UO  70  1=1 >NG2 
VALUES (JAY , I ) = VALUE ( I ) 

70  CONTINUE 
RETURN 

AO  FORMAT  (»0//////  EIGENVALUE  ERROR  ='»Dld.8/> 

ENO 


t 


( 
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SUBROUTINE  LfcP  (Y.X.M) 


•«••••••*•••»  LEGENDkE  polynomial  expansion  ***•***•*•*•*•*•••• 

DOUfiCE  PRECISION  Y.A.S 
DIMENSION  Vttl 

met. 

IF  (lx)  lb.Ib.20 
IU  RETURN 
2b  V (2)SX 

IF  (N-l)  10.10.30 
3b  DO  4b  1=2. N 
G=X»Y(I> 

40  YCWJ  )=G-Y ( l«l)4S*<(x-V ( 1»1 ) ) /FLOAT ( I*1  > 

RETURN 

END 


SUBROUTINE  FILLC  (T.PI.Pj.AJ.C.N) 

COMPUTE  COEFFICIENT  C(I.J) 

SEE  LIOU.  1973  (APPENDIX) 

OOUHCE  PRECISION  T.P1.PU.AJ.C 
DIMENSION  T(1).P1(1I.PU(1) 

C=0. 

DC  lb  L=1 >N 

IU  C=C*T (L)*PI (L)APJ(L) 

C=C*AJ/2. 

RETURN 

ENO 


SUBROUTINE  F ICCH  IB.C.UM.N.NL1) 


C 

C  SEE  Cl Out  1973  (APPENDIX) 

C 

BOUNCE  PRECISION  h.C.UM 

Li  I  MENS  I  ON  B(N.N)  »C  (N.N)  »UM(N> 

00  20  1  =  1. NCI 

IP1=1*1 

00  1(1  J=IP1.N 

h(I.J)=ClI.U)/UM(l> 

1U  h  ( J  •  1  )  =C  ( J 1 1 )  /UN  (01 
20  MI«1)  =  (C(I.I I— 1 .CO )/UM<  1 ) 
H(N.NI=(C(N.N>-1.D0)/UM(N) 
RETURN 
END 


SUBROUTINE  FICCS  (S.b.CtP.N) 

c 

C  S(NI  IS  TRACE  OF  B(N.N) 

C  SEE  C10U.  1973  (APPENDIX) 

C 

OOUhCF  PRECISION  S.B.P.C 
DIMENSION  SU).H(N.N)  tP(N.N).C(N.N) 
DO  l(i  1=1  .N 
S(I)=0. 

S(1)=S(1)+B(1.I> 

00  Id  J=l.N 
10  C  ( I  •  O )  =B  ( I  .  J ) 

NHSN 

IF  (M00(N>2) .EO.O)  NH=N-1 
IF  (N.CE.2)  GO  TO  50 
DO  90  1=2. Nh. 2 
CACC  GMPKD  (N.C.P.N.N.N) 

00  2U  J=1.N 
20  S(I)=S(1>*P(J.JI 

CACC  GMPRO  IB.P.C.N.N.N) 

IPISl'fl 
DO  30  J= 1 . N 

30  S(IP1 >=S(IPl)*C(U.J) 

90  C0N11NUE 

50  IF  (NH.EO.N)  RETURN 

CACC  GMPRO  (R.C.P.N.N.N) 

DO  NO  JS1.N 
GO  S(N)=S(N>+P( J.J) 

RETURN 

END 
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SUBROUTINE  6MPRU  (A.b'R>N«M.l) 

THE  SUBROUTINE  IS  TO  COMPUTE  6ENERAL  MATRIX  PRODUCT. 

DOUBLE  PRECISION  R.A.B 
DIMENSION  AU)ifKI)iK(l) 

IRSO 

IKS*M 

DO  10  KSI.L 

IRSlAAM 

DO  10  JsltN 

1RS1R+1 

JlSU-N 

IHSIK 

R  < IN  ISO 

DO  10  Isl«M 

JlSUl+N 

1RS1H+1 

10  R<IR>SRIIR)*ACJ1)*B(1B> 

RETURN 

END 


SUHHOUT INE  FIELD  (D.S.N) 

oil)  ARE  COEFFICIENTS  OF  THE  EXPANDED  POLYNOMIAL  OF  EIGENVALUES. 
SEE  LIOU.  1973  IAPPENLIXI 

DOUBLE  PRECISION  D.S.X 
DIMENSION  S( 1 ) >0 ( 1 ) 

O(l)i-SU) 

DO  HO  IS2»N 

1L1SI-1 

XsO. 

00  10  JS1.1L1 
10  XSX*D(1-J>*S(J) 

80  DU)S-U*S(1))/FL0AT(I) 

RETURN 

ENO 
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SUBROUTINES  • OPKHM *  AND  •DPnFB'  ARE  FOR  THE  CALCULATION  OF  ROOTS 
(EIGENVALUES )  ACCORDING  TO  POLYNOMIAL  EXPANSIONS 


DIMENSION  C ( 1 ) .«R (1 )  «RC ( t ) .POL ( 1 ) >0 (4) 

DOUBLE  PRECISION  L.HK.RC.P0L.Q.EPS.A.8.H.Q1.G2 

TEST  ON  LEADING  2ERO  COEFFICIENTS 

EPS= 1.0-6 
LlMslOO 
1RS1L41 
10  18=18-1 

IF'  11R-1)  420.420  «2u 
20  IF  (CdRl)  30.10.30 

»08K  UP  2ERU  ROOTS  ANL  NORMALIZE  REMAINING  POLYNOMIAL 

30  IER=U 
JSIR 
L=0 

ASC(IR) 

DO  nO  1=1. IR 
IF  (L)  40.40.70 
40  IF  <CC1>)  60.S0.60 
SO  HR< llsO.DO 

Rcmso.no 
Pol  (wi)so.oo 

JSd-l 
GO  TO  GO 
60  L=1 
ISTsl 
o=0 

70  JSJ41 

cm=cm/A 

PGL(Jl=C(l) 

Call  OVERFL  (Ml 
IF  (N-21  420. AO. GO 
AO  CONTINUE 

STaHT  RaIRSTO*  ITERATION 


01=0. 00 

O2SO.n0 

40  IF  (o-2>  330.100.140 
C  UE6REE  OF  POLYNOMIAL  IS  EQUAL  TO  ONE 
100  AsPOL(l) 

KRUSTIS-A  j 

RC(1ST>=0.D0 
18=18-1 
L2SO.O0 

IF  dR-1)  130.130.110 
110  DO  120  1=2. IR 
Q1SW2 

Q2SPOL (141 ) 

120  POL  1 1 )=A*b2401 
130  P0LdR41lsA4Q2 
GO  TO  340 

C  l*EfakEE  OF  POLYNOMIAL  IS  GREATER  THAN  ONE 
140  DO  220  LSI. 10 
NS  1 

ISO  GdISQl  . 

6(21=62 

CALL  DPwFR  (POL. J.6.LIM. I ) 

IF  (II  160.240. 230 

160  IF  (01)  1A0.170.1GO  i 

170  IF  (02)  lH0.210.lb0  *■ 

1A0  60  TO  (190.200.190.210).  N 
190  61S-fal 
NSN41 
GO  TO  150 
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200  62S-62 
N=N*1 
60  TO  150 

210  uisi.no-»ii 
220  62=1. no-02 

1ER= 3 
Ikslk-J 
RETURN 
230  IERsi 
200  01=0(1) 

02=0(2) 

H=0.l)0 

*=0.00 

1SJ 

250  H=-*l*B-u2*A*P0L(l) 

POL(I)=B 

B=A 

ASH 

1=1-1 

IF  (1-2)  260.260.250 
260  P0L(2)=H 
POLI1ISA 
LSlK-l 

IF  (w-L>  270.P70.i90 
270  00  260  1=J.L 

2AU  POL ( 1  -1 )  =P0L  ( 1-1 )  ♦PC.L  ( 1 )  •02+P0L  (  !♦!»  *61 
200  POL ( L ) =POL I L ) ♦POL ( L*  1  > •QP+01 
POL ( 1R  ) SPOL C IR  ) +02 
HS-. 500*02 

ASHah-01 

H=050RT(0ABS(A) ) 

IF  (A)  300.300.310 
300  RR(1ST)SH 
RCUST)=8 
ISTslST+l 
RR ( 1ST )SH 
RC(1ST)=-B 
60  TO  320 
310  B=H*OSI6N(B.H) 

RR(15T)S61/B 
RC( 1ST  1=0.00 
IST=IST*1 
RRUSTISN 
RC ( 1ST) -0.00 
320  1ST=1ST*1 
JSJ-2 
60  TO  90 
330  1RS1R-1 
3*0  A=0.b0 

DO  360  1=1 .IK 
61=C(I) 

02=P0L( !♦! ) 

POL ( 1 >=62 

IF  tol)  350.360. 350 
350  62=161— 62) /61 
360  62SOARSI62) 

IF  (62— A)  360.3A0.370 
370  AS62 
3A0  CONTINUE 
1=IR*1 
POL (11=1.00 
RR(1)=A 
RCC I >=0.00 
IF  HER)  390.390.610 
390  IF  (A-CPS)  410. *10. *00 
*00  IERs-1 
*10  RETURN 
*20  1ERS2 
1R=0 
RETURN 
ENO 
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SUBROUTINE  UPwFH  (C.IC.C.LIm.IER) 

DIMENSION  Cll)iU(l) 

D OUNCE  PRECISION  A.h* AA.RB.CA.CB.CC.CD* A1 .Bl.Cl .H.HH.01 .02*001 ' 
1002.4001 .Q0tt2.D01»D62. EPS. EPSl.C.O 
lERsu 
JSlC«l 
10  Jsj-i 

IF  lu-1)  400.400.20 
20  IF  ICIJtl  30.10.30 
30  AsClu) 

IF  IA-1.U0)  40.60.40 
40  DO  SO  Is l.J 
C ( I isC( I )/A 
CALL  OVEHFL  CN) 

IF  IN-2)  400. SO. SO 
SO  CONTINUE 

60  IF  lu-31  410.380.70 
70  EPSsl.0-14 
EPSlSl.D-6 
4=0 
LLSO 
01=0(1) 

62=0(2) 

001=0. DO 
002=0.00 
AASC(l) 

BH=L (2) 

CHSOaRSIAA) 

CASUARS(RH) 

IF  (LH-C«)  00. 90. 100 
AO  CCsCA+CH 
CRSCN/CA 
CAS1.00 
60  To  110 
90  CC=C»*CA 
CAS1.00 
CBS 1.00 
GO  TO  110 
100  CCSCaiCA 
CASCA/CR 
CB=1 .00 
110  CDSCCa. IDO 
120  ASO.OO 
T)3A 
A  ISA 
R13A 
JSJ 

6001=91 

6602302 

OOlSHH 

Dfttth 

130  Ms-6148— 62* A*C( I ) 

CALL  OVEHFL  (N) 

IF  IN— 2 )  420.140.140 

140  BOA 
ASH 

If*1:  }-l)  180. iso. 160 
iso  HXO.OO 

160  MS— 61481— 024 A14M 
CALL  OVEHFL  (N) 

IF  IN-2)  420. 170.170 
170  C1S81 
B1SA1 
A1SH 

60  TO  130 

180  MsCaaOABSIA)*CB4DABS(R) 

IF  ILL)  190.190.390 
190  LSL41 
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IF  (DABS ( A >-£PS*DABS(C ( 1 ) ) )  200(200(210 
200  IF  (0AB&<f>)-EPS«DABS(C<2))l  390*390(210 
210  IF  (H-CC)  220*220(230 
220  AASA 
6858 
CCS* 

BO  Iso  1 
662=0? 

230  IF  (L-UIM)  280(280*290 
200  IF  (H-CO)  930(930(230 
250  IF  10(1)1  270(260(270 
260  IF  (0(21)  270(920(270 
270  0(1  ISO. DO 
0(2)50.00 
60  10  70 

280  HHSDOAXI (DABS) A1 ) .DA6S(B1 ) (0ABSCC1 ) I 
IF  (HH)  920(920.290 
290  A1SA1/HH 
H1S«1/HH 
CISCl/Hh 
H5A1*C1-Hl«81 
IF  (8)  300(920(300 
300  ASA/HH 
BSA/HH 

HH5(h*A1-A*B1)/H 
H=(a*C1-H*B1)/H 
61561 FHH 
02SO2AH 

IF  (riABS<HH)-CPS*CARS(Bl)>  310(310(330 
310  IF  ((»ARS(H)“EPS*OABS(82)  )  320.320(330 
320  LL=1 

60  It)  120 

330  IF  (L-l)  120(120(390 

390  IF  (OABS(HH)-EPSl*0ABS(6l>>  350(350(120 
350  IF  (r.AB&(H)-CPSl«0A8S(02))  380(360(120 
380  IF  (DA8S(BB61«HH)-0AHS(B1*0U1))  370(960(990 
370  IF  (f(A66<0Bfa2«N)-(.'A8$(a2«0a2)l  120(990(990 
380  BdlsCtll 
«(2)sC(2) 

0(3)50.00 
0(9)50.00 
RETURN 
390  0(1)501 
•(2)502 
O 13)58 
0(9)50 
RETURN 
900  IER5- I 
RETURN 
910  ICRs-? 

RETURN 
920  ICRs- 3 

60  10  990 
930  IER51 
990  0(1)5661 
B(*)5W62 
6(3)— AA 
6(9)snb 
RETURN 
Ef(0 
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SUBROUTINE  TtoKAD  (JAT .KTHK.NLTR) 

c  •m«*»**«m«mmmmm»mmm»m****m*m»uh«***»*»****mmmhm*m 

C  SUBROUTINE  'TBftAU'  FOrMS  The  MATRIX  OF  LINEAR  EQUATIONS  RESULTING 
C  FKOM  TMt  LAYERED  CLOUD  STRUCTURE  AND  BOUNDARY  CONDITIONS  IMPOSED 
C  bY  THE  SURROUNDING  NON-SCATTER IN6  MOLECULAR  ATMOSPHERE  IN  LOCAL 

c  thermodynamic  equilibrium.  the  principle  calculations  are 
c  performed  in  this  routine  mith  the  resultant  output  of  the  spatial 
c  distribution  of  radiant  intensity  according  to  the  DISCRETE  rays 
C  (*NG2')  above  AND  BELOm  The  CLOUD.  VALUES  FOR  THE  ZENITH  ANGLtS 
C  0  AND  90  UE6  FOR  THE  TRANSMITTED  AND  REFLECTED  RADIANCE  ARE 
c  cxtrapolated  from  a  least. souares  poly,  fit  to  the  calculations. 
c 

c 

PARAMETER  LEVsAO.  aNTs6,  MZsZS.  N6sI6.  NG2s6»  N6lMs**8 

•  iNMwSTslEV^IO  .INT2=2*INT 
REAL  NUfcVtNBASE 

QOUbcE  PRECISION  Aa*01FF»0IFNB*ENTenSiFN0R.FUNC. INTERS. LE»LI»ORTf 
10SADA>PIn0.PULI «HT.SADA) TAUiN* W1.W2«XPL»XXR. YPL> YYR.UM.A. 

2 VALUES 

COMMON/1nPUT/ALT(NmOST) .TEMP(NMOST) .PRE(NMOST) .H20CNM0ST) . 

•  THKNS 1 9  > . NTH ICR  t NBASE , TRN ( NMOST  * INT  > »NE*LEV (10) 

•  • TauABS (NMOST t INT) 

COMMON  /angle/  UM(n6) .A(NG) .PI 

COMMON/LdMULT /UBC  » aTMd , aTMR  * ATMA i SFC  » THRUSv ( InT ) , THRUPT ( INT ) t ITER 

•  *ECb» IPCHAN(INT) tCLUEM 
COMMON  /POLY/  PINO(NG) 

common  /klNOON/  UP(NG2.!NT),D»(NG2'INT),THETA(NG).NU(INT).LV(!NT>. 
1CK1 ( INT ) ,CK2 ( INT ) #cUVERT ( INT >  * COVERT ( INT) 

COMMON  /CLOUo/  TEMfCUO) 
common  /EI6EN/  VALUES ( INT .NG2) 

COMMON/SfCEMS/EMIS ( INT ) ,NB 
COMMON/Ne«LEV/AODH) (10) .NUMNEMtLVNUM 
COMMON/Fc AGS/ 1 FLAG ( 2 ) 

COMMON  /RADATA/  PlNI(NG.lNT) .PT(INT) .BEXT(INT) 

DIMENSION  YPC(NG) . YYR(NG) i1NTENS(3.NG) fPOLI (NQ.NO) 

DIMENSION  LE(NGLM) .MKNGLMtNGLM) .W2(NGLM'NgeM) .ACOSUMiNG) 
DIMENSION  SAdA(NG’NG) »0SADA(NG*NG) »0PTH(9) .W(NG»NG) *LI (NGLM) 
DIMENSION  BINTC (3) . TAUI3) *ENTENS(3.N6) »RT (N6) »SAVE(INT2) 

DETERMINE  extinction  coefficients 

DO  10  ISX.NG2 
11&I*NG2 

RT(I)sVAtUES(JAY,I) 

RT(11)s-rT(I) 

10  CONTINUE 

COMPUTE  LEGENDRE  POLYNOMIAL 

DO  30  Isl'NG 
DO  20  Lsl'NG 
LlsL-1 
XPLSUM(I) 

CALL  LEP  ( YPl#XPL*l1) 

POLKI.L)STPL(L) 

20  CONTINUE 
30  CONTINUE 


uuo  uuo 


UO  40  1=1. NG 
UMI=Sh6L(UM(l)) 

ACOSUM (I)SACOS(UMI) 

40  CONTINUE 

compute  saoa 

00  60  Ksi.NG 
00  60  Lsl.NG 
H=L-1 
XXRsRT (K) 

CALL  USAUA  (YYR.XXK.L1) 
SAOA(K.L)=YYML) 

50  CONTINUE 
60  CONTINUE 

f INO  RIGHT  6  FOR  the  SMALLEST  ROOT 


IF  (N6.LT.10)  GO  To  170 
UO  160  N<»SsN62.N62 
ORTsRT (NGS) 

ISSX 

c  own =0.0 

70  CONTINUE 

count=count*i.o 
RT (NOS) SORT 
00  90  LSI. NR 
LlSL-1 
XXRSRT (N6S) 

CALL  USAuA  (YYR.XXH.L1) 

SAOA(NGS.L)sVYR(L) 

IF  (1S.E0.1)  GO  TO  60 
SAOA  (  NGS .  L )  sOSAO  A  <  NGS  .  L  ) 

60  CONTINUE 
90  CONTINUE 
FUNCsO.GuO 
00  110  Isl.NG 
»(I.M6S)s0.0u0 
00  100  Lsl.NG 

M(I.N6S)sm(I.NGS)+PIN0(L)*P0LI (I.L)*SADA(NGS.L)/(1»0D0+UM(X)» 
IRT(NGS)) 

100  CONTINUE 


FUNCsFUNC*A ( 1 >  *M ( I .NGS ) /2. 000 
110  CONTINUE 

00  120  Lsl.NG 
OSAOA(NGS.L)sO.OOO 
00  120  Isl.NG 

120  OSADA(NGS.L)sOSAOA(NGS.U4AU>*l*OLI(X.L)*V(l'NGS)/2.000 
FNORsOSAgA(NGS.I) 

00  190  Lsl.NG 

OSAOA(NGSfL)sOSAOA(NGS.L)/FNON 
SAOA (NGS.L )sOS ADA (nGS.L) 

190  CONTINUE 

!Ss2 

ORTs- ( 1 ,00-FINO (1 | > /OSAOA « NOS » * I 
OIPNSSl.OO-FUNC 


OIFFsOABS  ( 1 .  DOf  NOR ) 

IF  (OIFF.LE. 1.0-7)  GO  TO  ISO 


IF  (COUNT.9E.90.0)  GO  TO  140 


GO  TO  70 

140  6RITE  (6.410)  OXFF 


GO  TO  400 


ISO  CONTINUE 


uuu  uo  oo 


00  1O0  Lsl.NO 

SADA ( N6S , L j sOSAOA < NOS , L ) 

LL=L.l 

KLsM0U(LL.2j 

N0S2sN6S*2 

IF  Ut.EQ.0j  SADA(u6S2,L>=SADA<NQS,D 
IF  -UL.Eu.lj  S  AD  A  ( N6S2 » L  J  s-SAD A  ( NOS » L 1 
160  CONTINUE 
170  CONTINUE 

00  172  Isl.INT 

J=XdNT 

SAVE ( I  }-PT (I ) 

SAVE(Jj=BEXTU> 

172  CONTINUE 

BSCAsPT (wAYj*BEXT (jAVj 

NB0TsN£BlEV<11 

NTOPsN£WLEV(KTHK*lj 

OT  AU AB=T AUABS (NBOT , JAY ) -TAUABS  (NTOP , JAY } 

SABSsUT  AUAB/THKNS  <  kTHK ) 

BEAT (JAY ) sBEXT (JAY)+BABS 
PT ( JAY ) =BSC A/BEXT (WAY ) 

PRINT  176. JAY  rPT (JmY) .8EXT ( JAY) 

OPTh(KTHk)sTHKNS(KTHK) 

GNSFlOAT (NLYRj 

OPThI=OPTH(KTHK)/On 

00  IBO  Isl.NLYR 

T AU  U ) =DbLE <  FLOAT  d ) *OPTH  X  *BEXT ( JA Y ) ) 

ISO  CONTINUE 

values  of  a  of  characteristic  eq.  for  searchxno  the  eisenvalues. 

00  200  I si. NO 
00  200  Ksl.NO 
•  d.KJsO.ODO 
00  190  Ltsl.NG 

N( I .K ) ( I »k ) ♦PINO (LL) *P0LJ ( I  .LLj  *SADA(K.LL  J/d.ODO+(UM(I j  j*RT (K) ) 
190  CONTINUE 
200  CONTINUE 

MAXIMUM  DIMENSIONS  OF  SOLUTION  MATRIX  s  NSL 
NOLsNLYRwNO 
00  210  1=1. NOL 
00  210  Ksl.NOL 
Mld.Kjs0.000 
210  CONTINUE 

FOR  a<-Uj  —  CLOUO  TOP  BOUNOARV  CONDITION.  USE  W(9.K)  TO  M(1*.K) 

00  220  Isl.N02 
II=I>N02 
00  220  Ksl.NO 
ald.KjSNdI.Kj 
220  CONTINUE 

IF  (NLYR.EO.l)  00  10  2H0 
LlSNLYR-1 
00  230  Lsl.Ll 
00  230  1=1. NO 

iX=UN62*<L-l}*NG 

00  230  Ksl.NO 
KKSK+NO* (L-l j 

aldI.KK}S«(I.Kj*OfcXP(-RT(KJ*TAU(Lj| 

•Idl.NKjS-aldl.KKJ 
230  CONTINUE 
200  CONTINUE 


U2 


FOR  W(*U|  —  CLOuO  BASE  BOUNDARY  CONDITION*  USE  W(1»K>  TO  W(8,K> 

00  250  Isl*N62 
IIsI*NLYR*N62* (NLYh-1 > *N02 
JO  250  Ksl*N6 
KKsK*(NLYR-l)*N6 

W1(1I#KK)SW(I*K)*0£XP{-RT(K>*TAU(NLYR>> 

250  CONTINUE 

jO  260  lsl*N6L 
DO  260  Ksl*N6L 
w2(I'K)s»l(I'K> 

260  CONTINUE 

EVACUATE  THE  *L*  coefficients  OF  the  SOLUTIONS  FOR  inteoral- 
ojfferential  TRANSFER  LOCATION. 

THE  CLOUD  TEMPERATURE  GRADIENT  IS  DETERMINED  FROM  ATMOSPHERIC  PROFILE 
INPUT  OATA. 

HOTsTHKNS(KTHK) 

call  TC  <NLYR*0PTHi.H6T#JAY,KTHK) 

UO  270  lsl#NLYR 
UsNU(JAV) 

BINTC U)sTEMPC( II 
270  CONTINUE 

00  2o0  Isl*N62 

Li ( I ) sOBLE (06 ( I. JAY ) -BINTC (1) ) 

the  array  *le»  is  used  for  emissivity  calculations  in  which  the 

UPmaRJ  ANO  DOWNWARD  INTENSITIES  (UP I OF)  OF  THE  MOLEc  ATMOS 

are  set  to  zero 

LEU)sOBlE(-BINTCU)) 

280  CONTINUE 

IF  (NLYR.EO.l)  80  TO  300 
00  290  Lsl'Ll 
DO  290  Isl'NO 
iJ*I*N62«.|L-l>*N0 
LI<II)=0bLElbINTC(L*l|-8lNTC(L|) 

LE(ll)sLKIl) 

290  CONTINUE 
SOU  CONTINUE 

DO  330  Isl*N62 
IF  (NLYR.6T.I)  00  TO  3i0 
11SI+N62 
60  TO  320 

310  II=1*N62*(2*NLYR-1) 

326  LI ( 1 1 ) sOBLE ( UP ( I . JA Y I -BINTC ( NLYR 1 1 
LE(II) sObLE (-BINTC (NLTR 1 1 
330  CONTINUE 

COMPUTE  UPWARD  ♦  OOWNW*RD  INTENSITIES 

THE  MATRIX  INVERSION  WILL  RETURN  THE  L  -  COEFFICIENTS 

NT0TsN6*NLYR 

CALL  SIMM  (Wl*LI*NTOT,KS) 

CALL  SIMu  (WZ'LE'NTOT'KS) 


C  INTENSITY (l.J)  IS  FOR  THE  ITH  LAYER  -  JTH  ANOLE  »MU» 
00  360  I si* NLYR 
DO  370  Jsl*N62 
JJ— J4N02 

1NTENS ( 1 • J I sO • OODO 
INTENS ( I  *  JJ ) so . OODO 
ENTENS ( I  * J ) sO .0000 
ENTENSdt  JJIsO.OODO 
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00  360  Ksl.NG 
NK-K 

AXsOtXP  t-RT  tK> *TAU (I)> 

IF  U.NE.l)  NK=KfNo*( 1-1) 

IF  (l.GT.l)  00  TO  340 
intenS(X«j)=inTenS(I»j)<*lI(NK)»m(j«k) 

ENTEnS (X » d ) stNTENS ( I » d ) 4LE ( NK  > •* ( d  »K ) 

60  TO  350 

300  INTthS (ltd) SXNTENS ( X • J ) 4lX <NK  >  *W  (d.K ) *OEXP (-RT (K ) *TAU ( X-l >> 
£NTENS(X.d)sENT£NS(l.dJ4Ll(N)OwW(d*K)*DEXP(-RT(K)wTAU(I-l>) 

350  X  NU  »S  ( I ,  J  J )  =  I NTENS  (  X  t  d  J )  ♦(.  X  ( NK )  •  W  ( dd  ,  K )  •  AX 
ENTERS  ( X  #  dd )  SENT  ENS  ( X  t  dd )  4LE  ( NX )  •  W  ( d  J  #  K )  •  AX 
360  CONTINUE 

XNTEnS ( X t J ) s XNTENS ( X » d ) 40BCE < B X NTC ( X ) ) 

INTERS ( X , dd ) sXNTENi ( X , d J ) 4DBLE ( BXNTC ( X ) ) 

EWTEnS (ltd)  stNTENS  ( 1  . d )  +DBLE  ( B INTO  ( X ) ) 

ENTEnS ( X  t  dd ) sENTENS ( 1 1  dd ) 4UBLE (BXNTC ( X ) ) 

3TC  CONTINUE 
360  CONTINUE 
XENX&Sl 
XCCKaOsO 
NOUMsX 
dOUMsS 

C  EXTRAPOLATE  'tNTENS*  (FOR  EMISSIVITY)  VALUES  AT  ThETa  b  00,  AND  90 
XUPUwNsO 

CALL  EXTRPO ( THETA  *  ENTENS • T VERT » NOUM t  dOUN  *  lUPDWN ) 

CALL  EXTRP9 ( THETA • ENTENS »  ThOR  X  2 » NOUN  *  dDUN  *  XUPQMN ) 

IF ( ClUEM . 6T . 0 . 1 )  Eh  XSSsT VERT/BX  NTC ( 1 > 

C  EXTRAPOLATE  'XNTtNS*  VALUES  AT  THETA  s  0  AND  9o  0E6REES  FOR 
C  IRaNSMITTEO  and  reflecteo  RAOXANCE  used  in  determining  the  cloud 

C  TRANSMISSIVITY  and  REFLECTIVITY 

C  THE  EXTRAPOLATED  CLEAR  COLUMN  RADIANCE (UP)  FOR  THETAsO  IS  *CUVERT* 
C 

XEMXSsO 

CALL  EXTRPO(THETA«I NT£NS * T V£RT * NUUM , dDUM 1 1 UPOWN ) 

CALL  EXTRP9 (THETA. XNTENS . Th©RXZ .NDUM . dDUM . XUPOWN ) 

XUPDVNsX 

NOUMsS 

CALL  EXTRPO ( THETA . INT£NS . R VERT »NDUM » dDUM . X UPDWN ) 

CALL  EXTRP9 ( THEY  A . XNT  ENS . RHOR IZ . NOUM . dDUM . XUPDWN ) 

IF (CUVERT (dAY).EO.O)  60  TO  384 
TRANSsTVtRT/CUVERT (dAY ) 

REFLtTsRvERT/CUVERT (dAY) 

384  THRUSV(dAt)sTH«UPT(dAY) 

IF ( ITER,  t.  4,0)  GO  TO  365 
IF(CUVERT(dAY).Nt.O)  GO  TO  383 
THRUPT (dAY )sk.O 

COVERT (dAY )sRVERT 
GO  TO  388 

383  THROFT (dAY ) sRVERT/CDVERT (dAY) 

385  WRITE  (6.420)  dAY.hU(dAY) 

WRITE  (O.430)  OPTHiKJHK).TAU(NLYR) 

WRITE  (8.440)  TVERT.RHORIZ 
UO  390  d-l.No2 
dd=N<j4l>d 

WRITE  (6.450)  THETA(d).XNTENS(X.d).THETA(dd).INTENS(NLYR,dd) 
390  CONTINUE 
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YYYslNTENSU,!) 

WRITE  (6,460)  ThORi2*RVERT 

write  (6,470)  emiss , trans » reflct 

IF ( CLUEM . ST . 0 , 1 )  EMI SSMsENTENS ( 1 , 1 ) /BINTC ( 2 ) 

WRITE  (6,460)  EMISSM,8INTC<2) 

CALL  TEMTOP  (  TEMP  ,  N0AS»E  *  KTHK ,  JAY ,  EMJSSM ,  TWERT ,  Y Y Y  , 
INU(JaY), REFLCT) 

IF(JaY.Ne.INT)  60  TO  400 


1F(1FLA6(2) ,tQ.O)  60  TO  400 

wRITt(20*1000)  L VNUM , IFLAG ( 2 ) , NBASe , THKNS (KTHK ) 
WHITE (20, 1010)  ( (TaUABS ( X , J) • lsl ,LVNUM) , J=l, INT) 
•KITE (20, 1020)  (TEMPlI),I=l,tVNUM) 

•KITE (20, 1020)  (ALT ( 1 )  • Isl ,L VNUM) 

WHITE (20, 1010)  ( (P1NI(1,J) ,1=1,NG) ,JSl,XNT) 

WHITt (20, 1010)  (BEAT ( I ) , 1=1 , INT) 

WHITE (20, 1010 )  <PT(I),Isl, INT) 

400  CONTINUE 

00  405  Isl.JNT 

J=I*1NT 

PT (I)sSAVE(I) 

BEXT(I)SSAVEU) 

405  CONTINUE 
RETURN 


175  FORMAT (1mO,*HEV1SE[,  5IN6LE  SCAT  ALBEDO  AND  BEXT  FOR  CHANNEL  ',12' 

•  /,•  SINGLE  SCAT  ALBEDOS', F5. 3, 2X, «bEXTs' ,F5.3) 

410  FORMAT  (10X,'SSSSS  ITERATION  FAILED  *****  FIFF  =»,E20.14) 

420  FORMAT  (//,' - SPECTRAL  REGION',  12,'-—  NU=',F6.4) 

430  FORMAT  (32X, 'CLOUD  THICKNESSs' ,F6.2, • (KM)  OPTICAL  ThICKNESSs', 
1FG. 5, /,36X,«THETA',6X, 'UPWARD  INTENSITY* ,5X, 'THETA' ,5X, 'DOWNWARD  I 
2NTENSITY ' ,/) 

440  F0RMaT(36X,« 0. 0000', 5x,F13. 5, 7X, '90.0000', 5X,F13.5) 

450  F0RMaT(36X,f«.4,5X,F13,5,5x,F9.4,5x,F13.5) 

460  FORmaT(37X,'90.OO0o',5X,F13.5,6X,'iS0.0OOO',5x,F13.5) 

470  FORMAT  (/,20X,'EM1SS1VITYS  • ,F9,6,5X, 'TRANSMISSIVITYs  ',F9.6,5X,'R 
lEFLECTIWlTYs  ''F9.*,/> 

480  FORMAT  (/,10X, 'MEAN  EMISSIVITY  s  »,F9.6, 

1'  CORRESPONDING  TO  T  s  «,F7.2,'DEG  K») 

1000  FORMAT(2l2,2F10,5) 

1010  FORMAT (350E15.7) 

1020  FORMAT (50F10. 5) 


END 
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SUBROUTINE  TC  (LA.OPTHI.  h&T  rOATi KTHK ) 


1  HIS  SUhkOUT INE  CALCULATES  THE  CLOUD  LATER  TEMPERATURES  FROM  THE 
INPUT  ATMOSPHERIC  TEMPERATURE  PROFILE »  HEIGHT  OF  CLOUD  BASF*  ANO  THE 
MJMHEh  OF  ISOTHERMAL  LAYERS  APPROXIMATING  THE  TEMPERATURE  GRADIENT 


PARAMETER  LEV=40.  InT=6.  M2=25»  NG=lb»  NG2=d,  NMOSTsLEVAIO 
REAc  NRASE 

COMMON/ INPUT/ aLT  (NMOST ). TEMP  INMOST > ,PmE (NMOST)  >H20(NM0ST)  . 

•  TMMmS(R) .NTHICK.l, BASE.  ThMNMOST.  INT) .NEWLEYUO) 

•  , TAUAHS (NMOST , INT  > 

COMMON  /CLOUD/  TEMPO l 10) 

COMMON/LRMULT/UbC . ATMD  * ATMR . ATMA.SFC , THRUSV  ( INT) » THRUPT ( INT ) » I TER 

•  «£LR, IPCHaNI INT J  <CLDEM 
DIMENSION  TLMPNT(IO) 

ASSUME  LINEAR  TEMP  GRADIENT  . . . .MONOTON1CALLY  DEC  M/  HEIGHT 

IuXnOTsnEwLEV ( 1 ) 

10XT0P=NE*LEV (KTHk+1 ) 

TEMPO  < 1 ) =TEMP 1 IUXTOP ) *CLn£M 
TEMPO (LA  >  =  TEMP ( 1DXBOT ) wCLOEM 
UTsTEmPc (LA  l-TEMPC ( 1 ) 

D2=FLOAT(LA)-l. 

LAI— LA— 1 
DO  10  1=2, LAI 

TEMPO ( I ) STENPC ( 1 )  ♦  ( OT/02 ) * ( FLOAT ( I -1 ) ) «CLDEM 
TO  CONTINUE 

IF  (oAY.GT.ll  GO  TO  30 
»RIIE  lb. 40) 

DO  2U  1=1. LA 

TEMPf.T  ( 1 » =TEMPC  ( I  >-273.  lb 
IF  (l.EG.l)  ahITE  <«,. 50)  I.TEMPNT(I) 

IF  (l.Eu.LA)  WRITE  (O.60I  I.TEMPNT(Z) 

IF  U.NE.l.AND.l.uE.LA)  WRITE  (6.70)  l.TEMPNT(I) 

20  CONTINUE 

30  continue 

RETURN 

40  FORMAT  (//.4HX.* CLOUD  LAYER  TEMPERATURES') 

50  FORMAT  (SIX. II, F7. 2.'  (TOP)*) 

HO  FORMAT  (SIX, II, F7. 2.*  (RASE)*) 

70  FORMAT  (SIX. II .F7.2I 
END 


u  o  u  u  u 
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SOHKotJT 1NE  TlmTuP < TEMP »  NhASE •  KTHK .  JAY .EMISSM •  X . Y * II » REFLCT t 


this  suhhoutine  finds  the  brightness  temperature,  at 

1HE  SATELLITE’S  POINT  OF  VIE*  1ft  A  CLOUDY  ATMOSPHERE 


PARmMFTEk  LtV=40.INI=6.  NMOSTsLEVFlO 
HEAL  NBaSE 

COMMON/ INPUT /ALT I  INMOST ) (TEMPI  INMOST ) .PRE INMOST) .H20 INMOST ) » 

•  THftftS  ( 4 )  •  NT  HICK  •  I.BmSE  1  *  TRNINMOST  » INT )  . NEWLLV  (10) 

•  .TAOAHSINM0ST.INI ) 

COMMON/i.t oLc. V  /  AUOHT  110)  •  NUMNE to .  L VNUM 

u i mens ion  t au i nmost . int ). temp i nmost  > . temtop i 2  > . t i sfc i 2 ) *  angmod ( 2 » 

EQUIVALENCE ITLN(1«1)*TAU(1*I)) 

HOLTS I .  J*j046E-lh 

CS2.OR7R2RE10 

CGNVnT=2.»HoLT*C*(U»*2) 

ANGMOD  1 1 ) Si . 0 

AftGMor>(2)Sl.U1110b 

T1SfC(1>SX 

TISFL(2)=Y 

LEVISLVNUM 

IP IoaY.Nl.I)  GO  TO  oU 
HU  1SFlsnE*LEVIKTHh»1) 

ISFlPIsISFC-U 
00  mo  l=i»2 
HISTtMPl ISFC) 

T  AU1  -T  All  1 ISFC  •  JAY )  •ANGMOD (L) 

TEMTOPlLlsO.U 

00  <Ai>  KslSFCPl.LEVl 

h2STEMP(M 

T  AU2=T Ao I K  t UA Y I • ANGMOD <  L I 

TEM 1 oP (L IsTEMTOP I L ) *0 . 5* I H1+R2 ) *  I TAU2-TAU1 > 

TAU1STAU2 
BO  HlShE 

HSsT i SFC IL) *T AU( ISFC « JAY ) 

100  TtMT uP (L ) sTEMTOP (L ) ♦HS 

*kITE(6.RO0)  TEMTOP 1 1 ) .TFMToP 12) 

000  FORMAT  I/.1UX.'  VERY.  R.  TEMP.  s* »FI2 .6. 10X» »R.TEMP. AT  B.S  =• i 
1F12.ci»/> 

RETuKN 

El.Q 
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SUBROUTINE  DSADA  ( Y  »  X  *  W ) 


OENERATION  Of  SAOA  FUNCTION 


DOUBLE  PRECISION  Y.A.PINO.FI.YI1  f 

COMMON  /POLY /  PINO (lb) 

DIMENSION  Yd) 

YtltSl.ODO 
IF  IN)  10.10.20 
10  RE ToNN 

20  Y(2)=-(l.D0-PINO<U )/X 
IF  »N-1>  10.10.30 
3(t  DO  A0  1=2. N 
ASl 

F 1SUHLE ( A) 

Y11=(F1-1.0U0)/FI*Y<1-1) 

Y(l41)=-(2.UUU*(Fl-1.000)4l.0D0-PINO(I)>/(FI»X)*Y<I>“Yn 
40  CONTINUE 
RETURN 
END 


SUBROUTINE  PF1LE 
C 

C  WRITE  ATMOSPHERE  PROFILES  TO  FILE  20 

c 

PARamETEm  LtV=40.  NUMCHN=6.NMOST=LEV410 

COMnoN/1nPUT/H{NMCST ) .T (NMOST) ,P INMOST ) .H20CNMOST) .THK(9I . 

•  NTMCK.NRASt  .TRNINMOST. NUMCHN)  .NENLL/UO) 

COMMON /FRE O/FnU ( NUMCHN ) 

COMmon/nc.»LEW/AOOIiT(10)  .NUMNEW.LVNUM 
COMmON/SFCEMS/EMI S (NUMChN I 

DIMENSION  T R ( NMOST . NUMCHN ) . MT ( NMOST >  » TP ( NMOST ) 

1STA I =NEHTRn ( f>, 'UASo.AX  ?0 . .F/0/TRK/100  .  •) 

PRINT  200. 1ST  AT 
DO  So  1=LVNDM.1.-1 
NSLVI.UM-141 
hi(m:hid 
TP(h)=Ttl) 

I.U  ho  J=1 .7 
L=H-w 

lF(o.Al .«>  L=N-J-j*ll\T(b/U) 

Tk(a.L>=ThN(1.J> 

SO  CONTINUE 

print  auh. (Fniuji .j=i. numCHn) 

DO  100  1=1 .LVNUM 

WHITE (2U.30UI.HTU).TP(I).H2C(!)»(TR(1.U).0=1.NUMCHN> 

PR  I  iv  T  hUO.HT(I).T)  (1) . H20( I ) . ( TR ( I »d) .J=l .NUMCHN)  • 

100  CONTINUE 

200  FORmAT()H1.»STATIS  OF  ASSIGN  FOR  FILE  20.  =  *»T3) 

SO o  FORMAT (F7.3*H>.1*e7. 3. 2X.7F13.fi) 

400  FORMATdhl. 'HEIGHT  TtMP  H20  '.2X.7C  TAU  ( •  .F5 .2 .  •  I  •)>  * 

SOo  FORMAT ( 1 X. F7. 3. F6.1.F7. 3, 2X.7F13.fi) 

RETURN 

EnD 


118 


uuuu 


BLOCK  DATA 


DATA  OF  ANTENNA  GAIN  OF  RADIOMETER*  ANO  CLIMATOLOGICAL  PROFILES 
OF  TEMPEHATURE  *  HUM IT) .AND  HFIGMT. 


PARAMETER  NUMLV=40#  NMOST=NUMLVM0*  NUMCHNS6 
COMmON/PHES/  P(NMCST) 

COMMON/FREG/FNU ( NDMCHN ) 

COMMON/ T ANOM/T (NMOST I • M (NMOST ) 

COMMON/HANDU/H(NMoST > >U(NMOST ) 

COMMON/GAIN/GANT  <  12  7  •  NUMCHN  ) 

Data  (P(l) .Isl.NUMLV)/ 

i IS# (30*  v50*  *60«  * 

70. • Ab. *100. *115. *135. *1 50. *200. *250. *300. *350* ,400.f430. *475. » 

500. • 570. #620. *670. • 700. *7t»«.  >850 .  *920.  ,9S0.  * 1000./ 

DATA  Fnu/22. 235. 31. 65. 52, 85*53. 85. 55. 45*37. 00/ 
DaTa(&ANT(I.1).I=1.127>/  03*0.* 

9h.5*  93.7*  ftb.3*  7a. 9*  63.9*  51.7*  38.5*  28.2*  18.9* 

11.7*  6.12*  1.53*  1.39*  3.21*  3.60*  3.64*  2.51*  1.33* 

u.065*  0.889*  1.43*  1.31*  1.0S*  0.355*  0.124*  0.707*  0.9* 

0.840*  0.556*  0.431*  0*595*  0.875*  0.984*  0.873*  0.703*  0.455* 

0.279*  0.364*  0.564*  0.707*  0.865*  0.733*  0.549*  0.361*  0.153* 

0.125*  0.321*  0.539*  0*565*  0.513*  0.457*  0.309*  0.279*  0.215* 

0.347 *  0.409*  0.520*  0*544*  0.621*  0.594*  0.645*  0.539*  0.477* 

0.422  / 

DATA (GANT (I *2 1 *1=1*127)/  63*0.* 

96.5*  93.7*  85.3*  7a. 9*  63.9*  51.7*  38.5*  28.2*  18.4* 

11.7*  6.12*  1.53*  1.39*  3.21*  3.60*  3.64*  2.51*  1.33* 

0.065*  0.889*  1.43*  1.31*  1.05*  0.355*  0.124*  0.707*  0.9* 

0.846*  0.556*  0.431*  0.695*  0.875*  0.984*  0.873*  0.703*  0.455* 

0.279*  0.364*  0.564*  0*707*  0.A65*  0*733*  0.549*  0.361*  0.153* 

0.125*  0.321*  0.539*  0*565*  O.S13*  0.457*  0.309*  0.279*  0.215* 

u.347*  0.409*  0.520*  0.544*  0.621*  0.594*  0.645*  0.539*  0.u77* 

0.422  / 

DATA (GANT (1*3). 1=1.1271/  63*0.. 

9*,. 5.  93.7.  85.3*  7a. 9.  63.9*  51.7.  38.5*  28.2*  18.9* 

11.7*  6.12*  1.53*  1*39*  3*21*  3.60*  3*64*  2.51*  1.33* 

0.065*  0.889*  1*43*  1.31*  1.05*  0.355*  0*124*  0.707*  0.9* 

0.846*  0.556*  0.431*  0*695*  0.875*  0*984*  0.873*  0.703*  0.455* 

0.279.  0.364*  0.584*  0.707*  0.A65*  0*733*  0.549*  0.361*  0.153* 

0.125*  0.321*  0.539*  0*565*  0.513*  0*457*  0.309*  0.279*  0.215* 

0.347*  0.409*  0.520*  0.544*  0.621*  0*594*  0.645*  0.539*  0.477* 

0.422  / 

OATA (GANT ( I *4) *1=1*127)/  63*0. • 

*  96.5*  93.7.  65.3*  7e*9*  63.9*  51.7*  38.5*  28.2*  18.9* 


11.7. 

6.12* 

1*53* 

1*39* 

3.21* 

3.60* 

3.64* 

2.51* 

1.33* 

0.065* 

0.889* 

1*43* 

1.31* 

1.05* 

0.355* 

0.124* 

0.707, 

0.9* 

0.846* 

0*556* 

0.431* 

0.695* 

0.875* 

0.984* 

0.873* 

0.703* 

0.455* 

0.279* 

0.364* 

0.584* 

n.707. 

0.865* 

0.733* 

0.549* 

0.361. 

0.153* 

0.125. 

0.321* 

0.539. 

0*565* 

0.513* 

0.457* 

0*309* 

0.279* 

0.215* 

0.347* 

0*409* 

0.520* 

0*544* 

0.621* 

0.594* 

0.645* 

0.539* 

0.477* 

0.422  / 
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QATAIGAUTU.S)  ,1=1,1271/  63*0.  , 

•  96.5,  93.7.  85.3.  78.9.  63.9.  61.7.  38.5.  26.?.  16.9. 

•  11.7.  6.12.  1.53.  1.39.  3.?1.  3.60.  3.64.  2. SI.  1.33. 

•  U.065.  0.689.  1.43.  1.31.  1.05.  0.355.  0.124.  0.707.  0.9. 

•  0.846,  0.556.  0.431.  0.695.  0.875.  0.984.  0.873.  0.703.  0.455. 

•  0.279.  0.364.  0.564.  ft. 707.  0.A65.  0.733.  0.549.  0.361.  0.1S3. 

•  0.125.  0.321.  0.539.  n.Sbb.  0.513.  0.457.  0.309.  0.279.  0.215. 

•  0.347.  0.409.  0.520.  0.544.  0.621.  0.594.  0.645.  0.539.  0.477. 

•  0.422  / 

OaTA (6AN  Hl.f.)  >1  =  1.127)/  63*0.  . 

•  96.5.  93.7.  8b. 3»  78.9.  63.9.  51.7.  36.5.  28.2.  16.9. 

•  11.7.  6.12.  1.53.  1.39.  3.21.  3.60.  3.64.  2.51.  1.33. 

•  0.065.  0.8A9.  1.43.  1.31.  1.05.  0.355.  0.124.  0.707.  0.9. 

•  0.84n«  0.556.  u,431.  0.695.  0.875.  0.9A4.  0.873.  0.703.  0.u5S« 

•  0.279.  0.364.  0.564.  0.707.  0.865.  0.733.  0.549.  0.361.  0.|S3» 

•  0.125.  0.321.  0,539.  0.565.  0.513.  0.457.  0.309.  0.279.  0.215. 

•  0.347.  0.409.  0.520.  0.544.  0.621.  0.594.  0.645.  0.539.  0.477. 

•  0.422  / 


DATA  (T(l).lSl.NUMLV)/ 

•  226.2.  244.2.  256.7.  266.2.  262.2.  258.2.  251.2.  246.2.  242.2. 

•  235.2.  231.2.  2i7.2«  223.2.  221.2.  219.2.  214.2.  215.2.  212.2. 

•  211.2.  209.2.  216.3.  216.3.  216.3.  218.4.  2?5.5»  233.4.  240.4. 

•  251.7.  255.4.  2o9.9.  £62.2.  268.4.  272.4.  276.2.  278.3.  283.6. 

•  287.5.  290.7.  291.6.  293.6/ 

OATo  (wm.lSI.NUMLV)/ 

•  20*0..  .66690-3. 

•  o . 6 1 SbE-3 . . 9703E-3 • . 4807E-2 . . 302 IE-1 . . 8973E- 1 » . 1 703E-0 . . 2957E00 . 

•  0 . 366SEU0 . . 5653E00 . . 6897EQ0 . . 1204E4 1 . . 1B0BE+  1 . . 2639E* 1 . . 313BE+ 1 . 

•  0 . 53U6E+ 1  •  .  7582E41 .  .  1  OloE+2  .  .  1 140E+2  .  .  1 348E+?/ 


DATA  (H(l).Isl.NUMLWt/ 


6n.2S« 

61.28, 

54.30, 

46.85, 

45.60, 

43.35, 

40.24, 

38,27. 

37.1b. 

34.14, 

31.79, 

26.93, 

27.10, 

25.53, 

24.20, 

20.89, 

19.71. 

18.73, 

17.48, 

16.44, 

15.55, 

14.54, 

13.88, 

12.08, 

10.59, 

9.358. 

8.268, 

7,465, 

6,930. 

6.186, 

5.796, 

4.773, 

4.102, 

3.472, 

3.11b, 

2. 228. 

1.505, 

.8297, 

.5533, 

.1113/ 
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